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1.0  INTRODUCTION 

We  examine  in  this  report  the  problem  of  wideband  detection  of  a  spread  spectrum 
signal  from  an  interceptor's  viewpoint;  our  aim  is  to  derive,  evaluate  and  compare  various 
detector  structures  whose  purpose  is  to  intercept  a  spread-spectrum  communication 
transmitter.  The  transmitter  under  surveillance  employs  a  variety  of  modulation/ 
spreading/transmission  techniques  which  are  invariably  assumed  to  emit  the  message- 
bearing  signal  in  deep  background  noise  or  interference.  Thus,  the  challenging  task  of  the 
intercepting  detector  is  to  reveal  as  best  as  possible  any  transmission  of  the  unknown 
spread  waveform  in  the  presence  of  an  obscuring  noisy  environment 

Naturally,  the  degree  of  success  achieved  by  the  interceptor  when  detecting  the 
presence  of  the  communicator's  spread  signal  depends  on  the  amount  of  information 
available  to  him  regarding  the  structure  of  that  signal.  At  one  extreme,  the  interceptor's 
most  fortunate  situation  would  be  to  acquire  the  spreading  code  itself  C'crack"  the  code).  If 
this  is  too  demanding,  he  must  settle  for  less,  such  as  approximate  knowledge  of  the 
signal's  carrier  center  frequency,  code  rate,  code  epoch,  spreading  bandwidth,  etc.,  or  a 
subset  thereof.  At  the  other  extreme,  he  might  know  almost  nothing,  ^  in  which  case,  he 
could  resort  to  a  simple  energy  discriminator  device  (i.e.,''radiometei^)  since  it  is,  on  many 

C 

occasions,  a  low-SNR  asymptotically  optimal  detector  [1,2].  In  this  report  our  main 
interest  is  to  investigate  the  possibility  that  the  interceptor  could  do  better  by  optimally 
processing  whatever  information  is  available  to  him  between  the  two  extremes  cited  above  ^ 

The  detectability  of  the  primary  candidates  for  low-probability-of-intercept  (LPI) 
waveforms,  such  as  direct-sequence  (DS),  frequency-hopping  (FH),  time-hopping  (TH) 
and  their  hybrids,  is  customarily  related  to  two  factors:  (1)  performance  level  of  the 
interceptor's  detector  whose  only  function  is  to  monitor  the  communicated  messages  and. 
(2)  amount  of  signal  processing  (i.e.,  complexity)  associated  with  such  performance. 

‘Approximate  spectral  band  location  is  a  minimum  prerequisite  for  any  interceptor. 


In  subsequent  sections,  it  will  be  shown  that  the  detectability  (exploitability)  of  a 
spread  signal  is  greatly  dependent  on  the  interceptor's  knowledge  regarding  the  form  of 
both  the  waveform  and  the  interference.  In  other  words,  one  should  know  not  only  the 
spreading  format  used  (plus,  possibly,  some  other  parameters),  but  it  is  also  very  helpful  to 
be  able  to  identify  the  kind  of  background  noise  or  interference  involved.  Although  most 
of  this  report  deals  with  white  Gaussian  noise  as  the  dominant  form  of  (onmi-present) 
interference,  a  case  is  made  in  section  5.5.  about  the  impact  of  nonwhite  interference 
(random  tones,  in  particular).  It  will  be  shown  in  Sections  3  through  5  that,  when  white 
noise  is  the  only  deterrent,  the  gains  achieved  by  intelligent  processing  of  the  received 
waveform  can  vary  from  modest  to  significant  (always  as  compared  to  the  performance  of  a 
radiometer),  depending  on  the  scenario  at  hand.  However,  when  random  interfering  tones 
are  present  and  a  DS  or  a  hybrid  FH/DS  is  detected,  those  gains  could  become  impressive 
(many  tens  of  decibels,  for  example).  For  details,  see  Section  5.5. 

In  searching  for  those  structures  which  perform  the  aforementioned  tasks,  we  start 
from  the  optimal  solutions,  as  derived  from  applying  the  optimal  likelihood-ratio  rule.  The 
results  assume  a  good  deal  of  knowledge  about  the  signal  parameters  (although  never  the 
spreading  code  itself)  and  could  thus  become  rather  academic  in  most  practical  situations. 
Furthermore,  they  typically  suffer  from  the  common  symptom  of  prohibitive 
implementational  complexity.  Thus,  it  is  important  to  consider  suboptimal  structures  and 
evaluate  their  relative  loss  with  respect  to  the  (unattainable)  optimal  performance  ;this  has 
indeed  been  the  spirit  which  permeated  the  present  tqsic. 

The  report  is  organized  as  follows:  In  the  remaining  part  of  this  section  we  classify  the 
various  receivers  and  briefly  lay  the  theoretical  groundwork  for  evaluating  the  forthcoming 
detector  structures.  In  Section  3,  4  and  5,  we  develop  the  theory  for  DS,  TH  and  FH 
signals,  respectively. 

There  exists  a  wide  variety  of  interceptor  receiver  structures,  a  rough  classification  of 
which  is  attempted  below  and  is  also  depicted  in  Figure  1.  The  various  factors 


incorporated  in  the  classsiHcation  include  the  spread  spectrum  signal  whose  interception 
(detection)  is  sought,  the  optimality  or  suboptimality  of  the  interceptor's  decision 
algorithm,  the  degree  of  knowledge  of  different  system  parameters,  etc.  Wideband 
detectors  can  first  be  classified  as  optimal  or  suboptimal.  We  term  "optimal"  those 
structures  which  result  from  a  straightforward  application  of  the  generalized  likelihood  ratio 
theory  to  this  problem  and  invoke  a  number  of  assumptions  regarding  the  signal  structure. 
For  instance,  optimal  solutions  typically  assume  knowledge  of  certain  parameters  (code  rate 
or  hopping  rate,  SNR,  carrier  frequencies,  etc.)  and  average  over  the  unknown  ones 
(timing  epochs,  time  hopping  or  code  sequences,  etc.).  Since,  however,  optimal 
structures  are  often  hard  to  implement,  one  resorts  to  "suboptimal"  solutions  whereby 
either  one  (or  more)  assumptions  are  removed  or  certain  "atypical"  nonlinearities  such  as 
the  In  cosh  (•)  function,  are  substituted  by  simpler  ones  as,  for  instance,  a  square— law 
device. 

A  second  classification  results  from  the  nature  of  the  spread-spectrum  signal  for  which 
the  intercepting  device  is  designed,  so  it  can  be  a  DS,  FH,  TH  or  hybrid  detector.  Also, 
DS  and  TH  detectors  can  be  coherent  or  noncoherent,  depending  on  whether  or  not  the 
carrier  phase  is  known.  ITf  detectors  arc  assumed  to  be  noncoherent  from  hop  to  hop  due 
to  the  nature  of  the  communicating  channel  and  the  transmitter  wideband  FH  synthesizers. 

The  presence  or  absence  of  timing  (epoch)  information  about  the  spreading  code 
distinguishes  detectors  as  synchronous  or  asynchronous.  It  is  clear  that  timing  is  initially 
unavailable  since  the  intercept  detector  docs  not  even  know  if  the  signal  is  there;  hence, 
asynchronous  structures  naturally  attract  more  practical  interest.  However,  synchronous 
and/or  coherent  detectors  will  also  be  considered  —  not  only  because  they  provide  useful 
upper  bounds  in  performance  -  but  also  because  of  the  conceptual  possibility  of 
improved  detector  structure  aided  by  the  recursive  estimation  of  pertinent  parameters  [9]. 
Furthermore,  by  quantizing  the  continuous-epoch  uncertainty  region  into  an  adequate  finite 
number  of  points,  matched-detector  structures  can  be  devised  that  arc  synchronized  to  each 


of  those  levels  and  whose  outputs  are  combined  to  yield  the  final  decision  variable.  Thus, 
since  the  maximum  code-epoch  uncertainQr  equals  a  chip  or  hop  time,^  one  concludes  that 
good  nonsynchronous  detectors  can  be  constructed  from  combinations  of  a  few 
synchronous  ones.  Mote  discussion  on  this  topic  is  provided  in  later  sections. 

Finally,  we  classify  the  detectors  according  to  whether  they  expect  to  operate  in  the 
presence  of  white  (and  typically  Gaussian,  such  as  thermal  noise)  or  nonwhite  interference. 
The  maximum-likelihood  approach  renders  itself  to  the  first  case  most  easily,  while  its 
formulation  is  cumbersome  for  the  second.  Accordingly,  the  theoiy  is  fairly  mature  for 
AWGN  interference,  while  it  is  still  at  the  exploratory  level  for  scenarios  which  involve 
more  complicated  interference  forms.  It  is,  nonetheless,  a  fact  that  most  of  the  challenging 
practical  problems  reside  in  the  "dense-interference-environment"  area,  where  one  or  more 
wideband  signals  are  to  be  intercepted  amidst  a  large  number  of  irrelevant,  narrowband  and 
wideband  users,  jammers,  etc.  Such  interference  should  be  excised  to  the  best  possible 
degree  before  the  remaining  detection  and  estimation  tasks  on  the  sought  signal  can  be 
performed  successfully.  We  shall  not  deal  extensively  with  the  excision  problem  in  this 
report  (a  special  case  will  be  discussed  in  Section  5;  rather,  we  shall  assume  that  tic 
remaining  interference  after  excision  consists  of  wideband  Gaussian  noise  and  focus  on  the 
performance  of  associated  wideband  detection  algorithms.  In  terms  of  the  block  diagram  in 
Figure  2,  which  shows  the  simplified  flow  of  operations  in  a  standard  LPI  receiver,  the 
theory  herein  pertains  to  the  "wideband  signal  detection"  box,  while  interference  excision, 
noise-level  estimation,  etc.  are  performed  in  the  first  box  C'data  collection  and 
processing"),  possibly  aided  by  a  continuously  updated  "signal  history  file".  The  ultimate 
goal  of  the  LPI  receiver  is  indicated  in  the  last  box  ("signal  analysis  and  modulation 
recognition"),  wherefrom  the  information  is  made  available  to  other  interfaces. 


^This  clearly  assumes  a  purely  random  code  sequence  for  which  a  full  chip  (or  hop)  time- 
shift  corresponds  to  another  realization  of  the  same  stochastic  process. 
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We  shall  typically  refer  to  the  communicator's  signals  under  interception  threat  as 
LPI  signals.  For  the  sake  or  clarity  we  should  deflne  here  what  we  mean  by  such  an  LPI 
signal:  it  is  a  spread-spectrum  waveform  whose  code  is  never  identiflable  (decipherable) 
from  an  interceptor's  viewpoint;  the  latter  is  then  forced  to  wideband  processing  throughout 
(as  opposed  to  the  communicator's  narrowband  processing  after  despreading).  The 
detectability  of  such  an  LPI  signal  is  a  function  of  both  scenario-dependent  factors 
(topology,  distances,  power  budgets,  etc.)  as  well  as  scenario-independent  quality  factors. 
namely  the  structure  of  the  interceptor's  signal-processing  algorithms  (detection  rule).  The 
first  five  sections  address  the  algorithmic  aspects  of  the  overall  detectability/covertness 
issue  and  establish  associated  figures  of  merit  parameterized  by  the  spreading  format, 
existing  SNR  and  model  dimensionality  (time-bandwidth  product).  The  ultimate  quantity 
of  interest,  to  be  extracted  from  the  analysis  of  these  detection  algorithms,  is  the  required 
carrier-to-noise-ratio  CNR  =  (S/No)req  at  the  input  to  the  detector  in  order  to  achieve  a 
certain  performance  level  (see  Section  2).  Then,  any  value  of  (S/Nq)^  can  be  linked  to  the 
aforementioned  scenario-dependent  factors  for  a  global,  LPI  vulnerability  assessment.  It  is 
clear,  however,  that  the  two  aspects  of  detectability  can  be  treated  separately.  Furthermore, 
different  detection  algorithms  can  be  compared  on  the  basis  of  their  associated  (S/Nq),^, 
their  relative  merit  increasing  as  (S/No)req  decreases;  conversely,  from  the  communicator's 
viewpoint,  spreading  waveforms  with  lower  (S/No)req  are  more  vulnerable  to  interception. 
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2.0.  LKEUHOOD  RATIOS  AND  PERFORMANCE  MEASURES 


It  is  well  known  [3]  that,  under  the  common  detecdon-performance  criteria  (most 
notably,  the  Neyman-Pearson  philosophy),  likelihood  ratios  yield  optimal  hypothesis¬ 
testing  solutions,  and  performance  is  measured  by  the  resulting  pair  of  detection  and  false- 
alarm  probabilities  (Pd,  Pfa)-  Furthermore  it  will  be  shown  later  that,  in  many  models  the 

overall  likelihood  ratio  (LR)  is  expressible  as  the  product  of  individual  LR’s  due  to  the 
statistical  independence  of  successive  code  chips  (random  DS  code)  or  hopping  slots  (TH 
or  FH).  When  the  true  model  of  the  problem  does  not  allow  for  such  independence 
assumptions  (as  in  the  case  of  a  random,  but  fixed,  carrier  phase),  the  resulting  LR  is  too 
complicated  for  either  exact  analysis  or  implementation  3.  In  this  case,  the  performance  of 
those  (complicated)  optimal  decision  schemes  is  assessed  by  upper  and  lower  bounding 
arguments  (see  the  following  sections);  in  conclusion,  the  product  LR  can  always  be 
thought  of  as  the  analytical  cornerstone.  Since  the  log-likelihood  ratio  (LLR)  is  not  only 
theoretically  equivalent,  but  also  more  convenient  to  implement  than  an  LR  (it  turns 
products  into  sums),  it  will  serve  as  our  point  of  departure. 

The  number  of  terms  entering  the  summation  of  the  LLR  hops  is  typically  large  —  it 
equals  the  number  of  chips  (DS),  or  time  frames  (TH)  or  frequency  (FH)  observed. 
Therefore,  it  is  customaiy  to  model,  via  a  central-limit-theorem  (CLT)-type  argument,  the 
distribution  of  the  resulting  LLR  as  approximately  Gaussian  under  either  hypothesis.  This 
constitutes  a  basic  assumption  upon  which  performance  is  derived  in  this  report. 
However,  the  fact  that  no  such  LLR  possesses  an  fi2ta£l  Gaussian  distribution  brings  about 
certain  ramifications,  as  explained  below. 


^In  addition,  it  can  be  argued  that  the  supeiposition  of  random  data  on  DS  or  TH,  such  as 
PSK  modulation,  invalidates  the  modeling  assumption  of  a  fixed  value  of  the  signal  phase 
over  a  multi-bit  observation  interval. 


Furthermore,  it  should  be  pointed  out  here  that  the  ensuing  discussion,  which  is 
perfectly  suited  for  DS  waveforms,  is  also  applicable  to  FH  or  TH  only  if  the  ovwall  LR 
includes  a  large  number  of  independent  hops  (FH)  or  frames  (TH),  so  that  the  product-LR 
and  some-LLR  forms  are  valid  and  useful.  As  a  suboptimal  alternative  to  such  optimal  FH 
or  TH  detection  procedures,  it  is  possible  (and  many  times  preferable,  from  an 
implementation  viewpoint)  to  perform  individual  LR  decisions  per  hop  (FH)  or  frame  (TH) 
and  combine  afterwards  these  "tentative"  LR-decisions  over  many  hops  or  frames  (post 
detection  combing).  In  such  cases,  they  one  should  address  the  quantitative  properties  of 
the  LR  directly;  this  will  be  further  elaborated  upon  later. 

Let  us  now  return  to  the  sum-LLR  framework.  To  be  specific,  assume  N 
independent,  identically  distributed  (i.i.d.)  measurements  rj,  j  =  1,2,...,N  and  let  A(r) 

indicate  the  individual  LR  associated  with  any  one  of  them.  The  corresponding  individual 
LLR  is  denoted  by 

l(r)  ■  In  A(r)  (2.1) 

The  overall  LR  and  LLR  A  and  I  respectively,  are 

N 

A  -  n  A(r.)  (2-2a) 

i-i  ' 

and 

N 

^  "  I  “-tv  (2.2b) 


due  to  the  i.i.d.  assumption. 


Let  X(r)  denote  any  "well-behaved"  transformation^  of  r.  Although  it  is  convenient,  in 
this  context,  to  think  of  X<r)  as  an  approximation  to  (or  simplification  of)  t(r)  as  per  (2.1), 


it  need  not  be  so.  Clearly, 


X  =2  A.(rj) 

j=l  ^ 


constitutes  another  (typically  suboptimal)  decision  statistic  for  the  hypothesis  testing  at 
hand.  Furthermore,  the  CLT  approximation  will  be  applicable  to  both  land  X,  as  long  as 
N  is  sufficiently  large.  With  the  resulting  "Gaussian"  performance  or  (2.2)  or  (2.3)  in 
mind,  consider  the  following  model:  Let  X  denote  a  Gaussian  decision  statistic  under  both 
hypotheses  Hji  i  =  0.1.  The  performance  pair  (Pp,  Pp^^),  associated  with  the  test 


X  >  threshold 


tU-:: 

»  I  • 

i  i 

..v.y 
-y. .%  N 


is  simply 


m  „  -  m 
x,0  X, 


where  Q(x)  is  the  Gaussian  integral  function 


r  1 

=  l  .r5=  exp  {-— }  dz 


VTjt 


^For  our  purposes,  this  means  that  the  transformation  possesses  finite  first  and  second 
moments  under  either  hypothesis. 
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with  Q-i  [x]  its  inverse,  and  j,  i  =  0,1  are  the  means  and  variances  of  X  under 

either  hypothesis.  In  many  cases  where  the  inclusion  of  a  weak  signal  does  not  alter  the 
variance  a\  Q  significantly,  i.e.  a\o  *  c\i ,  equation  (2.5)  simplifies  to 


™x.l  ■  ”’x.O  A 


=  d, 


'x,0 


(2.7) 


where  the  deflection  (or  distancel  d,  has  been  appropriately  defined  in  (2.7).  Thus  under 
the  equal-variance  Gaussian  assumption,  d,  quantifies  performance  uniquely.  In  other 
words,  any  desired  performance  level  (Pq,  Pfa)  equivalently  translates  to  a  specific 
"desired"  distance  d,  via  (2.7);  the  converse,  howevw  is  not  true. 

The  above  discussion  applies  equally  well  to  X  representing  a  "Gaussian"  t  or  any 
other  X.  Noting  that,  by  definition. 


(m  -  m  ) 

i.l  1.0 


Nd, 


(2.8  a) 


l.O 


and  similarly 


d^  =  Nd^,  (2.8  b) 

X  X 

where  the  distances  d£  and  d^  pertain  to  the  individual  statistics  of  £and  X,  respectively, 
we  conclude  that  (2.5)  -  (2.8)  provide  a  set  of  approximate  performance  estimates,  reliable 
to  the  degree  that  the  Gaussian  assumption  is  itself  sufficiently  accurate.  Simulation 
results,  to  be  discussed  later,  confirm  a  very  satisfactory  accuracy  for  large  values  of  N. 

Use  of  the  aforementioned  distance  measure  hinges  upon  the  evaluation  of  the 
appropriate  moments.  That  is  no  particular  problem  for  certain  standard  forms  of  X  (r);  in 
fact,  d2  can  be  identified  as  an  output-SNR  measure  for  familar  quadratic  detectors  (see, 

for  instance,  [4]-[8]. 


On  the  other  hand,  the  exact  calculation  of  can  be  a  formidable  task  indeed. 
Fortunately,  this  situation  is  remedied  within  the  Gaussian  framework  due  to  results  by 
Peterson,  Birdsall  and  Fox  [3].  For  the  sake  of  continuity  and  compactness,  the  main 
result  is  summarized  below  without  substantation;  however,  the  analytical  reader  is  urged 
to  go  through  Appendix  A  of  this  report  Therein,  the  tool  of  generating  functions  is 
employed  to  provide  an  independent  proof  of  the  aforementioned  result,  Chemoff  bounds 
on  detector  performance,etc. 

Theorem  (P.B.F.,  1954):  Consider  an  LLR  I  which  is  Gaussian  under  Hq.  Then 
(a)  I  will  be  Gaussian  under 


(b)  -nij„o  =m|^,, 


(c) 


where  d3,is  given  by^ 
G 


d^  =  ln(l+var{A|H })  *  ln^{A2|Ho}  (2.9) 

G  0 

and  A  is  the  LR  associated  with  1 .  Interestingly  the  moments  A  of  are  sometimes  much 
easier  to  obtain  than  those  of  t, ,  hence  the  usefulness  of  (2.9).  The  quantity  d^  will  be 

referred  to  as  the  "Gaussian"  distance. 

Corollary:  With  N  i.i.d.  measurements,  performance  is  uniquely  quantified  by 


d J  .  N  In  (1  +  var{A  |  »  (2. 10) 

It  is  interesting  that  the  Gaussianess  of  £  under  Hq  suffices  to  specify  everything  else  in 
terms  of  the  single  number  (^.  Furthermore,  (2.10)  has  been  the  basis  of  analytical 
evaluations  for  detectors  employing  optimal  nonlinearities,  whenever  N  is  thought  to  be 
sufficiently  large  to  support  an  approximate  Gaussian  assumption.  We  shall  follow  the 


^Here  €[•}  and  var{*}  stand  for  the  mean  variance,  respectively. 


same  approach  here,  quickly  cautioning  however  against  two  possible  sources  of 
inaccuracy  and  discrepancy,  namely  (I)  different  variances 


"to"  "^1  * 


under  two  hypotheses,  which  would  negate  conclusion  (c)  of  the  theorem,  and  (11) 
intermediate  values  of  N,  coupled  with  low  PpA  and/or  high  Pq,  which  would  simply  put 

the  accuracy  of  the  Gaussian  approximation  into  question.  It  is  clear  that,  for  case  (I) 
above,  (2.5)  would  provide  a  more  reliable  answer;  nonetheless,  the  difficulty  in 
analytically  evaluating  mj^,  oj  ^ !  i  =  0.1  remains.  As  for  (H),  it  might  be  desirable  to 
proceed  with  a  numerical  calculation  of  (Pd,  Pfa)  by  computer  via  the  analytical  evaluation 
of  characteristic  functions,  or  resort  to  bounds.  Again,  the  nontrivial  aspect  of  this 
approach  is  the  closed-form  evaluation  of  the  characteristic  functions  for  the  kinds  of 
optimal  nonlinearities  encountered  here. 

In  contrast  to  the  above,  suboptimal  rules  involve  dj  ,  which  is  typically 
straightforward  to  calculate.  It  is  also  quite  useful  for  the  purpose  of  comparing  different 
suboptimal  decision  rules,  applied  to  the  same  hypothesis  testing  model  (see  Section  3.3). 
Finally  we  shall  use  dj  (or  (2.5),  whenever  appropriate)  to  compare  with  d|  =  Nc^,  in 
order  to  assess  the  cost  of  suboptimality,  keeping  always  in  mind  our  aforementioned 
reservations* .  Note  again  that  dj  cannot  be  derived  from  the  PBF  theorem,  since  X  is  niit 
anLLR. 

Finally,  for  the  FH  and  TH  cases  mentioned  above,  where  the  appropriate  tests  come 
from  individual  LR's  and  not  from  LLR's,  it  is  clear  that  distance  measures  based  directly 
on  A  should  be  evaluated.  A  good  such  distance  candidate  is 


(2.11) 


^  var^/2  {AIHq} 

provided,  of  course,  that  A  is  now  the  statistic  which  is  approximately  Gaussian  and 
var  (AlHi)  *  var  {AIHq}.  Indeed,  this  is  the  case  in  certain  scenarios  to  be  explored  later, 
making  the  evaluation  of  dy^  in  (2.1 1)  a  necesssaiy  step. 

Using  properties  of  LR's,  expression  (2. 1 1)  can  be  simplified  to 

d^=  Vvar  {AIHq}'  =  (2-12) 

which  will  be  of  use  in  section  5. 

The  theoretical  groundwork  laid  out  in  this  section  will  now  be  applied  to  a  variety  of 


LPI  scenarios. 


3.0  DIRECT-SEQUENCE  WAVEF  ORMS 


In  this  section,  we  consider  and  analyze  wideband  detectors  for  DS  waveforms.  It 
will  be  assumed  throughout  that  the  receiver  has  knowledge  of  the  carrier  frequency  and 
code  rate  and  that  the  code  is  biphase  modulating  die  carrier.  For  the  carrier  phase,  the  two 
possibilities  of  it  being  known  (coherent  detecdon)  and  unknown  (noncoherent  detection) 
are  examined  separately  in  Section  3.1  and  3.2,  respectively.  A  common  assumption  in 
both  of  these  sections  is  that  the  detectors  are  synchronous,  i.e.,  the  code  chip-timing 
epoch  is  known.  In  most  practical  situations  this  is  rather  unrealistic;  the  synchronous 
results  can  then  be  thought  of  as  upper  bounds  on  the  performance  of  any  asynchronous 
detector.  Furthermore,  by  quantizing  the  coniinous-epoch  uncertainty  region  into  an 
adequate  finite  number  of  points,  matched-detector  structures  can  be  devised  which  are 
synchronized  to  each  of  these  points  and  whose  outputs  are  combined  to  yield  the  final 
decision  variable.  The  amount  of  quantization  required  depends  upon  the  chip  envelope;  in 
this  report  we  shall  restrict  attention  to  the  constant-envelope  case.  Since  the  maximum 
code  epoch-uncertainty  for  a  random  code  equals  a  chip  time,  one  concludes  that  good 
asynchronous  detectors  could  be  constructed  from  combinations  of  a  few  synchronous 
ones.  More  discussion  on  the  cost  of  asynchronism  is  given  in  Section  3.3.  Section  3.4 
discusses  some  gains  derived  from  correlation-combining  of  two  independent  receptions, 
ranks  the  performance  of  the  aforementioned  detectors,  and  compares  them  to  that  of  the 
radiometer. 

In  order  to  introduce  some  notation,  let  the  high-rate  ±l-valued  random  spreading  code 


c(t)  be  represented  by 


where  p(t)  is  a  unit  pulse  of  duration  seconds  and  {qi  is  a  sequence  of  i.i.d. 

random  variable  (rv's)  with  Pr{Cn  =  +1}  =  Pr{Cn  =  -1}  =  0.5.  Note  that  the  random 
sequence  {Cq},  as  modeled,  does  not  distinguish  between  a  pure  random  code  and  a  code 
which  is  additionally  modulated  (in  a  synchronous  fashion)  by  random  binary  data. 
Furthermore,  the  chip  epoch  is  modeled  by  the  rv  e,  uniformly  distributed  over  (0,1).  The 
waveform  observed  by  the  detector  is  therefore  given  by 


>V2S  c(t)  cos(o)q  t  +<>)+  n(t)  (Hj  ) 


r(t)  = 


(0<t>T) 


where  S,  (Uq  and  <|>  are  the  average  signal  power,  carrier  radian  frequency  and  carrier  phase, 

respectively  and  n(t)  is  bandpass  AWGN  with  one-sided  power  spectral  density  (PSD)  of 
No  W/Hz.  The  observation  time  is  T  seconds,  which  we  assume  to  be  an  integer  multiple 
of  the  chip  time,  i.e.,  T  =  NT^ ;  N  a  positive  integer.  Such  an  assumption  is  the  least 
restrictive  since,  in  practice,  N  is  generally  a  large  number  (of  the  order  of  hundreds  or 
higher).  However,  a  somewhat  stronger  restriction  is  embedded  in  (3.2).,  namely  the 
assumption  that  under  hypothesis  H],  the  signal  is  present  during  the  whole  observation 

interval.  That  excludes  the  possibility  of  the  signal  either  starting  or  ending  at  any  random 
time  in  (0,T).  Still,  such  a  formulation  is  important  because  it  provides  meaningful  and 
fairly  simple  comparative  conclusions  which  (a)  would  otherwise  be  obscured  by 
mathematical  complexities  and,  (b)  can  be  argued  to  extend  to  more  general  models. 

In  terms  of  the  above  notation,  a  synchronous  detector  implies  that  e  is  identically 
zero,  while  a  coherent  one  means  that  <t>  is  known.  The  detector  for  which  ^  is  unknown 
but  constant  over  the  observation  interval  duration  is  well  within  the  coherence  time  of  the 
spread-spectrum  channel.  For  computational  purposes,  we  also  consider  the  Hetitious 
chip-noncoherent  detector,  where  phase  is  assumed  to  be  totally  random  (uniform  in 
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(0,2}i)  and  independent  from  chip  to  chip.  Although  this  is  not  realistic  for  a  biphase- 
system,  it  is  discussed  here  because  its  readily  derived  performance  serves  as  a  useful 
lower  bound  to  the  performance  of  the  carrier-noncoherent  receiver.  In  addition,  the 
random  phase  model  becomes  increasingly  realistic  if  the  DS  modulation  has  four  or  more 
phase  states. 


When  both  the  random  phase  <l>  and  random  chip  epoch  e  arc  assumed  known,  the 

detector  is  asked  to  perform  the  following  composite  hypothesis  testing  problem;  Decide 
between  the  alternative  Hq  and  Hi,  where 


J  c(t)  +  n  ft) 


(H<,) 


(0  <  t  < 


Hypothesis  Hi  is  composite  because  it  contains  all  possible  patterns  that  the  code  can 
assume  in  NT^  seconds.  Although  we  consider  only  random  codes  (in  which  case  there  are 

exactly  2^  such  patterns),  some  of  the  steps  below  would  also  be  valid  for  deterministic 
codes.  The  subscript  I  in  njft)  indicates  that  only  the  inphase  component  of  the  noise 
contributes  with  a  flat  two-sided  PSD  of  No/2  W/Hz. 

Well-known  results  about  optimal  detectors  (see,  for  instance  [12])  can  be  applied  here 
to  yield  the  composite  likelihood  ratio  test: 

AfrdO  .  expj^  r(t)  c,(t)  dt| 


expj-NxJ  ,  -  I  "‘a 


where 
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is  the  predetection  (or  chip)  SNR,  CjCt);  i  =  is  the  i*  code  pattern,  Cy  is  the  j®  chip 
of  the  i*  pattern,  and  rj  is  given  by 


rj  i  r(t)  at 

(3.6) 

It  is  shown  in  Appendix  B  that  an  expression  equivalent  to  (3.4)  can  be  derived  in  terns  of 
the  cosh  (•)  function,  where  cosh  (x)  =  e*  +  c*  )/2.  It  is  given  by 

N  1  I  /  2^  \  *ll 

A(r(t))  »  '■jj  *0  (3.7a) 

which  yields  to  optimal  LLR  test 


In  A  (r(t))  *  -Ny  +  J  In  cosh  In  An 

c  iti  J/  %  0  (3.7b) 

as  shown  in  Figure  3,  we  note  that  (3.7)  could  have  been  derived  directly,based  on  the 
i.i.d.  and  AWGN  assumptions,  as  per  (2.2)  (see  also  [5],[6],[8]  for  similar  results).  The 
method  preferred  here  can  be  generalized  to  models  with  statistically  dependent  chips  (c.f. 
Section  3.2)  or  nonrandom  code 

It  can  be  shown  that,  the  conditional  mean  of  the  random  argument  of  the  In  cosh  (•) 
function  in  (3.7b),  assuming  Cj  »  1  under  Hj,  is  equal  to  2Yc>  which  is  also  its  variance 


under  Hq.  Therefore,  for  typical  predetection  SNR  values  below  -10  dB  or  so,  the 
approximation  In  cosh  (x)  «  x2/2  is  well  applicable,  in  which  case  (3.7b)  reduces  to  the 
approximate  (suboptimal)  rule 


N  p  H, 

■  'i  i 


depicted  in  Figure  3.  Since  the  aforementioned  approximation  of  the  In  cosh(x)  function 
becomes  increasingly  tight  as  x->0,  it  follows  that  (3.8)  is  asymptotically  optimal  as  Yg-K) 
Within  the  theoretical  framework  of  Section.2.  the  Dcrformance  of  (3.8)  is  easily 


^  .  . 


Therefore,  Pjj  and  PpA  are  related,  as  per  (2.5)  by 


Pq  =  Q 


(3.10a) 


with 


a  =  ^;  b  =  4  (synchronous  coherent  detector).  (3. 10b) 

Since  bYc  is  typically  much  less  than  one,  (3.10)  can  be  combined  with  (2.7)  to  yield 


dx  “  Yc 


(3.11) 


We  note  that  performance  is  dictated  by  the  product  » which  is  typical  of  schemes 

employing  post-detection  integration.  In  contrast,  detection  of  a  known  waveform  would 
involve  the  factor  VNYc;  the  losses  in  performance,  due  to  the  composite  test  (since  the 

signal  pattern  is  unknown)  for  low  values  of  SNR,  are  evident  by  comparison.  Note  also 
that  rule  (3.8)  could  have  been  analyzed  exactly,  but  this  issue  is  deferred  until  noncoherent 
detectors  are  discussed  in  the  following  section. 

In  terms  of  the  overall  observation  time  T  and  the  null-to-null  spread  spectrum 
bandwidth  Wj  =  2'ri  (Hz),  (3. 1 1)  can  be  rewritten  to  establish  the  required  (S/Nol^q 
in  order  to  achieve  the  performance  level  d^  as 


(3.12) 


The  above  refers  to  the  suboptimal  rule  (3.8).  The  Gaussian  distance  dj^of  (2.7)  for 
the  optimal  rule  (3.7b)  can  also  be  derived.  Let 
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with  Tj  as  per  (3.6).  Under  Hq,  yj  is  a  zero-mean  Gaussian  rv  with  variance  =  2Y(.- 
Using  the  fact  that^{cosh(yj)}  =  cxp{o2  tl}  in  conjunction  with  (3.7a)  results  in 


d|j  =  VN  In  cosh(2yc) 


(3.13) 
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For  small  x.  In  cosh(x) «  jtP-H,  which,  upon  substitution  in  (3.13),  verifies  that  dj-^  dx  of 
(3.1 1)  as  Yc->0.  In  conclusion,  the  approximate  rule  is  well  suited  for  the  low  prediction 
SNR  range  of  interest 

The  performance  (Pj),  PpA  =  10"^)  of  the  optimal  and  suboptimal  decision  rules  (3.7b) 
and  (3.8);  respectively,  as  predicted  by  (3.13)  and  (3.11),  is  shown  in  Figure  4.  The  two 
values  of  N  chosen,  a  high  (1000)  and  a  low  (SO),  are  meant  to  illustrate  the  qualities  and 
limitations  of  the  Gaussian  assumption.  The  PpA  =  lO'^  value,  used  here  for  exemplary 
purposes,  was  fixed  both  in  theory  and  simulation.  Although  this  Pp^,  if  considered  by 

itself,  would  correspond  to  unacceptable  false-alarm  rates  it  is  typically  true  that  multiple 
tests  are  employed  before  a  final  decision  is  reached,  thus  lowering  the  overall  Pp^  to  the 

desired  levels . 

The  simulation  results  for  both  rules,  also  included  in  Figure  4,  confirm  the  expected 
trend  of  the  CLT-Gaussian  approximation  to  improve  as  N  and  Pj)  increase,  for  fixed  Pp^. 

This  improvement  is  manifested  in  the  increased  quality  of  both  performance  estimates,  i.e. 
eqns.  (3.11)  and  (3.13),  each  pertaining  to  its  own  test.^  Interestingly,  the  approximate 
estimate  (3.1 1)  seems  satisfactory  for  both  optimal  and  suboptimal  tests.  This  should  be 

‘^For  a  10  Mchips/sec  code  and  N  *  1000,  the  number  Pp^  =  .01  corresponds  to  an  average 
false-alarm  rate  of  100  FA/second. 

^This  is  certainly  true  for  increasing  N;  for  small  values  of  N,  the  applicability  of  (3.13) 
becomes  more  questionable  as  Pq  increases. 


dj^  performance  prediction  for  LLR  test  (3.7b) 

approximate  Gaussian  prediction  (3.10) 

simulation  results  for  LLR  test  (3.7b) 
simulation  results  suboptimal  test  (3.8) 

(a)  Pq  vs,  for  the  synchronous  coherent  detector 

_p 

optimal  and  suboptimal  tests,  Pr«  =  10  . 


attributed  to  the  low  chip-SNR  Yc.  which  renders  rules  (3.7b)  and  (3.8)  almost  equivalent, 
as  attested  also  by  the  closeness  of  the  respective  simulation  results.  In  different 
environments  however,  (eg.  or  TH),  corresponding  to  higher  hop-SNR's,  the  cost  of 
suboptimality  should  be  more  pronounced  (see  Section  4  and  5). 


Let  us  now  relax  the  assumption  that  the  carrier  phase  of  the  DS  waveform  is  known 
by  the  interceptor,  but  still  retain  the  assumption  of  a  known  code  epoch.  We  shall  state  the 
optimal  decision  rule  for  a  carrier-noncoherent  system,  whose  CXACI  analysis,  however, 
appears  unattainable.  Instead,  a  chip-noncoherent  detector  will  be  derived  and  analyzed, 
whose  performance  can  serve  as  an  asymptotically  tight  (as  Yc“^0)  lower  bound  for  the 

carrier-noncoherent  system.^  Thus,  the  performance  of  the  latter  can  be  nicely  bracketed 
between  that  of  the  former  and  the  upper  bounding  performance  of  the  synchronous 
coherent  detector^®  in  the  previous  Section  3.1. 

Consider  first  the  carrier-noncoherent  system.  Under  the  usual  narrowband 
assumption  (Oq  »  2  tt/T,  the  resulting  decision  rule  is 


exp(-  Ny, 


yn  Ml 


(3.14) 


where  Io(*)  is  the  zeroth-order  modified  Bessel  function  and  Rj  is  the  i®*  correlation 
envelope 


*  V''/  *  . ^ 


(U5a) 


In  (3.15a);  the  two  components  are  given  by 


j?i  ''*j  ^ 


(3,15W 


^Since  a  chip-noncoherent  detector  utilizes  less  statisdcal  information  than  a  carrier- 
noncoherent  detector,  it  is  inferior  to  the  latter  for  all  SNR  values. 

•®In  fact,  it  was  argued  in  [4]  that  the  coherent  upper  bound  is  asymptotically  achieved,  for 
large  distance  d,  by  an  approximate  version  of  the  carrier-noncoherent  detector,  because  the 
(unknown)  constant  phase  can  be  estimated  reliably  at  high  SNR. 


cos  u^t 
sin  uQt 


(31.5c) 


'■Ij 

& 

r(t) 

:Qi 

(j-l)T, 
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The  optimal  rule  (3.14)  -  (3.15)  is  hard  to  mechanize.  Instead,  suppose  that  the  received 
waveform  r(t)  is  represented  under  Hj  by 


m 

r(t)  =  Cj  p(t -jTc)cos|u)Qt  ♦j)  +  n 


(t) 


(3.16) 


where  {(t>j}  is  a  sequence  of  i.i.d.  uniform  phases,  thus  modeling  the  aforementioned  chip- 
noncoherent  system.  It  is  then  straightforward  to  show  the  independence  of  chips,  using  a 
procedure  analogous  to  Section  3.1  or  directly  (see  also  [4],  [8]),  that  the  optimal  decision 
rule  is 


‘  -  -"V  *  I,  ■"  ‘o 

J  * 


liJi'o 


where  rj  is  the  envelope  of  the  j*  chip 


(3.17) 


with  r,  ,  rn  as  per  (3.15b).  Again,  under  the  small-argument  approximation 
'j  i 

in  In(x  »  x2/4,  it  follows  that  (3.17)  reduces  to  the  suboptimal  rule 


Rule  (3.19)  is  formally  identical  to  (3.8),  the  difference  being  that  tj  here 
corresponds  to  a  noncoherent-chip  integration.  Furthermore,  it  is  shown  that  as  Yc— >0 
(namely  for  small  arguments  of  the  Iq  (•)  function),  the  optimal  carrier-noncoherent  rule 

(3.14) ,  which  involves  averaging  over  all  possible  2^  N-tuples  of  the  code,  reduces  to  the 
simple  suboptimal  decision  rule  (3.19)  for  the  chip-noncoherent  system.  It  can  thus  be 
concluded  that  the  performance  associated  with  (3.19)  is  asymptotically  close  to  that  of 

(3.14)  on  the  low-SNR  side;  however,  this  probably  occurs  at  an  unacceptably  poor  overall 
performance  level. 

Although  the  distance  dj^  for  the  exact  rule  (3.17)  can  be  derived  (see  Section  4.0  for  a 
closely  related  case),  we  shall  be  content  with  d;^  of  (3.11)  in  view  of  the  low-SNR 

conclusions  of  the  previous  section.  The  random  envelope  rj  in  (3.18)  possesses  a 
Rayleigh  Pdf  under  Hq 


f  (r  jl  Ho)  - 
and  a  Rician  under  Hj 
fCtj  I  Hi )  - 


(3.20b) 

with  o2  =  NoT<./2.  Then,  the  statistics  of  the  r.v.  QJd^)  in  (3.19)  are  known  exactly  [19]: 
under  Hq,  (X/o^)  is  a  central  chi-square  r.v.  with  2N  degrees  freedom  while,  under  Hj, 
(X/o^)  is  a  noncentral  chi-square  with  2N  degrees  of  freedom  and  noncentrality  parameter 
8  =  2Nyc.  Following  Marcum's  expressions  [15]  for  the  performance  analysis  of  pulsed 


radar,  which  involves  an  analogous  detection  model,  we  conclude  that  the  probabilities  of 
interest  are 


Pfa  ■ 


(3.21a) 


where  I(u,p)  is  the  incomplete  Gamma  function  and 


* 


Pd 


Qn 


(3.21b) 


where  Qn(  a  ,P)  is  the  Nth  order  generalized  Q  function  [15, 16] 

Although  the  above  results  are  exact,  they  arc  cumbersome  to  use  from  a 
computational  viewpoint,  especially  for  large  N.  In  addition,  the  Gaussian  approximation 
yields  satisfactory  results  in  the  range  of  N  of  interest  and  thus  becomes  our  alternative 
simpler  route. 

It  is  shown  in  Appendix  D  that 


(3.22a) 


(3.22b) 


Combining  (3.19)  and  (3.22)  and  employing  Gaussian  assumption  results  in  (3.10a)  where 
now 


a=l;  b  =  2 
Equivalently, 


(Syncrochronous  chip-noncoherent  detector).  (3.23) 


(3.24) 
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which  reveals  a  loss  factor  of  V2,  or  1.5  dB,  with  respect  to  the  coherent  detectors.  This  is 
the  penalty  paid  for  the  lack  of  coherency.  Furthermore,  since  the  performance  of  the 
carrier-noncoherent  system  is  bracketed  between  the  coherent  (upper  bound)  and  the  chip- 
noncoherent  (lower  bound)  ones  that  differ  by  only  1.5  dB,  it  is  concluded  that  any 
implementational  complexity  beyond  that  of  rule  (3.19)  is  not  justified  from  a  practical 


standpoint. 
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Asynchronous  Detectors 


In  the  most  natural  problem  formulation,  the  code  timing  (epoch)  e  is  a  random 
variable  uniformly  distributed  within  a  chip  interval.  We  shall  first  show  in  this  section 
how  a  cascade  of  analog  devices  and  nonlinearities  can  be  configured,  in  principle,  to 
optimally  average  the  likelihood  ratio  over  the  epoch  uncertainty.  The  derivation  follows 
the  basic  steps  illustrated  by  Krasner  [8]  in  a  somewhat  simplified  and  expanded  manner. 
Then,  practical  suboptimal  detectors  of  reduced  complexity  will  be  derived,  analyzed  and 
compared  with  the  optimal  one.  The  main  purpose  here,  besides  obtaining  useful 
quantitative  performance  estimates,  is  to  establish  a  certain  performance-bounding 
philosophy;  although  attention  is  confined  to  constant-envelope  DS  signals,  the  line  of 
thought  is  also  applicable  to  nonconstant  envelope  waveforms,  for  which  the  effect  of 
epoch  mismatch  could  be  much  more  pronounced. 

For  the  sake  of  simplicity,  the  coherent  case  is  analyzed  herein.  Most  of  the  ensuing 
steps,  however,  carry  without  change  to  the  noncoherent  rule.  The  starting  point  is  the 
coherent  likelihood  ratio  (3.7a),  which  is  now  interpreted  as  a  conditional  LR,  given  the 
value  of  the  random  offset  e.  In  other  words,  we  can  write 


» 


AC'"(t)|e3  =  TT  expfyj.  cosh  /  _.^(-rj(e))| 
j=l  I  No  / 


where  the  coherent  integration  (3.6)  is  adjusted  to  reflect  the  knowledge  of  e. 


(3.25a) 


rjU) 


yjVc 


r(T)  dt 


(3.25b) 
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(3.26) 

Complicated  as  it  may  seem,  rule  (3.26)  can  indeed  be  implemented  exactly  by  the 
means  shown  below.  First,  we  define  the  waveform  y(t)  as 


■-■.vV-'- 


y(t)  =  J‘  r{x)  dx 

t-Tc  (3.27) 

This  process  y(t)  is  easily  mechanized  as  the  output,  due  to  an  input  r(t),  of  a  linear  niter 
matched  to  the  chip-envelope  pulse^i  p(t)  =  u(t)  -  u(t  -  T^),  since 


f  r(T)  dx  =  r  r(x)  p(t-,)  d, 
t-T 


(3.28) 


Note  from  (3.27)  that  the  argument  of  the  In  cosh(»)  function  in  (3.26)  is  just  y(jTp  +  e). 

N 

Thus,  the  next  step  involves  creating  the  sum  of  samples  ^  In  cosh  [y(jTg  +  e)]  of  the 

J-l 

new,  transformed  process  In  cosh  [y(t)l.  To  do  that,  consider  that  the  convolution 
of  In  cosh  [y(t)]  with  the  impulse  response. 


1,,  J 


.•« 


N-1  ,  V 

h(t)  -  I  «(t  -  ml  ) 
maO  ' 


(3.29) 
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^^Here,  u(t)  stands  for  the  unit-step  function. 
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namely. 


g(t)  =  h(t)  ©  In  cosh(]y(t)3  »  ft  lncosh(^(t-T)3dT 

•• 

N-1  ] 

*  2  In  cosh  y(t  -  mT  ')\ 

m=0  ^ 


which ,  sampled  at  t  =  kTj+  e  yields 


s(kT^+e)  »  f  In  cosh[y((k-ra)T^+G))J 


j»k-N+l 


IncoshjypT^+ejj 


(3.31) 


where  the  change  of  variables  j  =  k-m  has  been  used.  At  k  =  N,  i.e.  when  g(t)is  sampled  at 

N 

ts=  NTg  +  e,  the  desired  sum  ^  In  cosh(y(iTc  +  e)]  is  formed.  Since  0  S  e  ^  Tj.  the 
sampling  tunes  tj  belongs  to  the  interval  NT^  (N+1)  Therefore,  the  above  sum 
(or  its  exponent,  according  to  (3.26))  can  be  averaged  with  respect  to  e  by  simply 
integrating  g(t)  (or  exp{g(t)})  over  the  interval  NTg<t<  (N+l)Tc.  Again,  such  an 
integration  with  a  uniform  "weight"  function  can  be  directly  interpreted  and  performed  as  a 
convolution  with  the  square  chip-pulse  p(t)  similar  to  (3.28).  (A  possible  analog 
implementation  of  this  convolutor  with  a  square-pulse  impulse  response,  also  called  a 
"moving  average"  filter,  is  shown  in  Figure  5  .  It  should  be  emphasized  here,  however, 
that  the  use  of  the  same  pulse  shape  for  the  two  convolutions  is  coincidental,  owing  to  the 
square-pulse,  biphase-modulated  pseudo-random  code  considered.  Arbitrary  code  pulse- 
shapes  would  require  general,  matched-filter  structures  for  y(t)  in  (3.25)  -  (3.28),  while  the 
averaging  with  respect  to  the  "uniform"  r  v  e  can  always  be  performed  via  the 
aforementioned  square-pulse  convolution.  The  general  case  is  treated  in  depth  by  Krasner 
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The  coherent  optimal  LR  detector  structure  for  the  continuous-epoch  uncertainty  case 
is  shown  in  Figure  6,  explicitly  indicating  the  series  of  linear  filtering  and  nonlinear 
operations  required  to  implement  it  Note  that  the  impulse  response  h(t)  of  (3.29)  and  the 
subsequent  convolution  (3.30)  are  associated  with  the  tapped-delay-line  box  in  Figure  6, 
which  calls  for  (N-1)  delay  elements,  seconds  apart  Figure  7(a).  Such  an 

implementation,  however,  is  not  recursive  and  would  require  more  delay  elements  for 
higher  N.  As  an  alternative,  the  feedback  structure  of  Figure  7(b)  can  be  employed,  which 
only  needs  a  single  delay  element  in  the  feedback  path.  Ideally,  the  impulse  response  of  the 
latter  structure  would  be  an  infinite-length  delta  train;  thus,  N  need  not  be  specified  a  priori. 
However,  stability  considerations,  due  to  component  noise,  parasitics,  etc.  typically  require 
the  insertion  of  a  gain  a  <  1  (dotted  lines  in  Figure  7(b)),  which  practically  limits  the 
"effective"  length  NTg  of  the  resulting  impulse  response. 

As  witnessed  from  Figure  6,  a  number  of  (possibly  costly)  linear  and  nonlinear 
devices  are  involved  in  the  optimal  asynchronous  processor.  One  step  towards  reducing 
the  complexity  is  to  quantize  the  epoch  uncertainty  region  of  seconds  into  a  (small)  finite 

number  of  alternatives,  i.e.,  points.  This  is  equivalent  to  pretending  that  the  epoch  r  v  can 
take  on  only  those  values  and  accordingly  develop  the  optimal  or  near-optimal  detector  for 
the  resulting  finite  hypothesis  problem.  The  quality  (as  well  as  complexity)  of  such  an 
approximation  will  obviously  increase  with  the  number  of  points  considered.  It  will  soon 
be  clear,  however  that,  for  constant-envelope  waveforms,  the  simplest  case  of  only  a  two- 
point  quantization,  which  will  be  the  one  we  analyze  herein,  provides  quite  adequate 
performance  with  respect  to  the  ideal  (synchronous)  detector;  thus,  higher  level 
quantization  and  complexity  might  even  be  undesirable  from  a  practical  standpoint.  In 
addition,  a  suitable  interpretation  of  this  analysis  can  simultaneously  serve  as  a  bounding 
technique  (upper  and  lower)  for  the  optimal,  continuous-epoch  uncertainty  detector  whose 
exact  analysis  is  intractable. 


Let  us  assume  that  the  two  points  comprising  the  epoch  uncertainty  are  Cq  and  ej;2« 
corresponding  to  the  offsets  e  =  0  and  e  =  1/2  in  (3.1);  that  is  to  say,  the  incoming  code 
chips  are  likely  to  arrive  cither  in  perfect  synchronization  (  e  =  0)  or  half  a  chip 
off  ( e=  1/2)  with  respect  to  the  local  chip- timing  reference,  and  they  are  both  equally 
probable  to  occur.  A  variation  of  the  optimal-likelihood  ratio  (3.7a)  for  the  coherent  case, 
adapted  for  the  above  scenario,  would  imply  that  the  decision  rule  should  be 

A  W  rr  ,  /2/s  N  , 

■  2  r"',!  'jleoD)  ♦  ..pa  In  cosh  (— 


where 


axp{-Ny^} 


(3.32) 


*  C-1.T,  4 


(3.33) 


Clearly,  the  chip-noncoherent  detector  would  form  the  envelopes  rj  in  place  of  the  coherent 

integrations  (3.33)  and  would  proceed  in  an  analogous  way  from  equation  (3.17).  We 
should  note  here  that  the  two  NTg-sccond  observation  intervals  for  the  two  corresponding 

exponential  terms  of  (3.32)  are  slightly  off  by  half  a  chip.  For  large  N,  such  "edge  effects" 
are  insignificant  and  will  be  neglected  in  the  following.*^ 

Equation  (3.32)  can  be  simplified  by  employing  the  small-argument  approximation  to 
the  In  cosh  (•)  function  as  in  Section  (3.1),  resulting  in  the  (non-LR)  test 


^^So,  at  this  pont,  the  epoch  uncertainty  is  assumed  to  be  discrete  by  nature,  not  by 
quantization. 

i^For  that  purpose,  small  adjustments  will  be  made  for  convenience  without  explicit 
acknowledgment 
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In  view  of  the  findings  of  Section  (3.1)  it  could  be  argued  that  the  test  (3.34)  albeit 
suboptimal,  should  perform  close  to  the  LR  test  (3.32).  Nevertheless,  it  is  shown  below 
that  the  former  can  be  outperformed  by  another  suboptimal  test,  at  least  under  certain 
circumstances,  which  is  in  fact  a  further  simplification  to  (3.34).  In  other  words,  one 
should  not  automatically  conclude  that,  in  the  realm  of  suboptimal  tests,  a  "simpler"  rule 
necessarily  implies  a  deterioration  in  performance.  Towards  that  end,  let  rji  and  Tj2  indicate 

the  coherent  integrations  during  the  first  and  second  halves  of  the  j***  chip  interval, 
respectively,  i.e. 


(i-HTc 


(j-l/2)T. 


r(T)dT 


(3.35) 


and  define 
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Then,  (3.34)  can  be  written  as 


^  "  (l)[*>‘p{A5}  ♦  -  PxpCZ}  cosh(W)  >  L, 


or,  through  the  transformation  X,  *  In  t. 


X,  -  2  ♦  Incosh(W) 


^  threshold 


(3.39) 


which  is  thus  equivalent  to  (3.34) 

An  approximation  (3.39)  is  obtained  by  simply  omitting  the  second  term, 


Z  ^  Lg  threshold 


(3.40) 


which,  by  substituting  (3.36)  into  (3.37)  expanding  terms  and  absorbing  constants  into  the 
threshold,  reduces  to 
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'■jZ  ^j2  ('jl  ^  ^  threshold 

J  A  Hn 

®  (3.41) 


An  examination  of  the  decision  rule  (3.41)  reveals  the  source  of  inferiority  of  the 
asynchronous  coherent  detector  in  comparison  with  its  synchronous  counterpart,  namely, 
rule  (3.8).  The  latter  rule  can  be  obtained  from  (3.41)  if  the  factor  rj^j  j  in  the  last  term  is 
substituted  by  rji,  because  the  resulting  expression  is  then  the  perfect 
square  (rjj  +  of  (3.8).  Since  rji  and  rj^.^  arc  independent  r.v.’s  (they  belong 

to  different  random  chips)  regardless  of  the  true  epoch  value,  it  follows  after  some  thought 
that  such  a  substitution  would  increase  the  mean  of  the  decision  statistic  X2  under  and 

slightly  reduce  the  variance-the  effect  of  which  is  a  net  improvement  in  performance  for 
synchronous  detector. 

We  now  turn  to  evaluating  and  comparing  (3.39)  and  (3.41).  Let  us  first  look  at  X,2 
and  return  to  Xi  later.  Without  any  loss  in  generality,  we  can  assume  that  the  true  epoch 
is  £  =  0.  The  symmetry  of  the  problem  suggests  that  the  results  to  be  derived  under  such 


conditioning  also  represent  the  unconditional  (average)  performance  in  this  two-possible- 
value  case.  Using  the  familiar  type  of  Gaussian  analysis,  it  can  be  shown  that 

which,  when  compared  with  its  synchronous  counterpart  (3.1 1)  reveals  a  loss  factor  of 
V4/3,  or  0.63  dB,  as  the  SNR  penalty  (in  y^)  associated  with  the  aforementioned 

asynchronous  detector. 

The  above  estimated  loss  can  also  be  viewed  as  a  lower  bound  on  the  performance 
losses  of  the  continuous -epoch-uncertainty  detector  illustrated  in  Figure  6.  This  is  because 
the  two-point  uncertainty  coherent  detector  discussed  so  far  assumes  a  random  epoch  that 
can  take  on  only  two  values;  thus,  it  faces  less  uncertainty  than  any  other  multipoint  or 
continuous  type  of  epoch-randomness.  Therefore,  it  is  bound  to  outperform  any  other 
asynchronous  detector,  including  the  aforementioned  continuous-uncertainty  one. 

It  is  also  possible  to  derive  an  upper  bound  on  the  performance  losses  of  the 
continuous-uncertainty  detector.  This  can  be  done  by  assuming  a  truly  continuous  epoch 
and  viewing  the  two-point  detector  as  a  suboptimal  quantized  implementation  of  the  optimal 
continuous-uncertainty  one.  The  performance  losses  of  the  suboptimal  detector  under  a 
worst-case  choice  of  the  unquantized  epoch  will  then  clearly  upper-bound  the  losses  of  any 
other,  higher  complexity,  multipoint  or  continuous  detector. 

Some  reflection  could  persuade  us  that  the  worst  epoch  for  the  two-point  quantized 
detector  is  the  one  straddling  between  e=0and  e=  1/2,  namely,  e=  1/4  (or,  equivalently, 
e=3/4,  whose  performance,  however,  is  identical).  This  is  because,  under  Hj,  rji  will 
then  have  the  lowest  possible  mean  (namely,  zero)  half  of  the  time,  while  rj2  remains 
unaffected.  Any  other  offset  will  yield  a  higher  overall  average  for  rji ,  thus  improving 
performance.  The  steps  of  this  worst-case  analysis  are  summarized  in  Appendix  E, 
wherein  it  is  shown  that 


2 
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Thus,  the  upper-bounding  loss  is  (4  V^),  or  1.42  dB,  with  respect  to  the  synchronous 
structure.  In  summary,  lack  of  synchronism  in  the  coherent  case  costs  anywhere  from 
0.63  dB  (lower-bound)  to  1.42  dB  (upper  bound)  Tighter  bounds  can  be  obtained  by 
reiterating  the  above  arguments  for  epoch-quantized  detectors  with  more  than  two 
quantization  points.  As  the  number  of  points  increases,  one  asymptotically  converges  to 
the  actual  performance  of  the  continuous-epoch  detector.  In  any  case,  a  rough  figure  of  an 
average  1  dB  loss  due  to  asynchronism  should  not  be  far  from  exact.  Similar  conclusions 
can  be  drawn  for  noncoherent  detectors. 

Finally,  let  us  briefly  return  to  the  test  (3.39).  In  order  to  assess  the  impact  of 
the  In  cosh(W)  term,  consider  the  rv  W  as  defined  in  (3.36),  (3.37),  and  assume  a  worst- 
case  offset  (e  =  1/4).  It  is  clear  from  the  symmetry  of  the  problem  that  the  rv's  Ag  and  A  j/2 

have  identical  first  and  second  moments  under  either  hypothesis. 

Furthermore,  for  large  N,  the  summations  in  (3.36)  will  make  them  approximately 
jointly  Gaussian,  correlated  rv's.  It  is  then  well  known  that  the  transformations  defined  by 
(3.37)  render  Z  and  W  approximately  Gaussian  and  independent  rv’s.  Thus,  Z 
and  In  cosh(W)  are  also  approximately  independent.  The  mean  of  W  is  zero  and  its 
variance  o^|h  is  essentially  the  same  under  either  hypothesis  H^.  Since,  for  zero-mean 
W,  the  expected  value  e{ln  cosh(W)|H^}  is  only  a  function  of  follows  that  the 

second  term  in  (3.39)  does  not  contribute  to  the  difference  of  the  means  of  Xj  under  Hq 
and  Hj  while  it  does  contribute  to  the  variance.  The  net  result  is  that  d^i  < 

detector  is  better  off  by  neglecting  In  cosh(W).  This  conclusion,  is  true  for  the  worst  offset 
e  =  1/4;  although  other  offsets  will  require  somewhat  more  complicated  arguments,  the 
same  supremacy  of  over  Xj  is  likely  to  hold. 

Some  asynchronous  simulation  results  are  presented  in  Figure  8  for  (a)  N  =  SO  and 
(b)  N  =  1000,  along  with  the  synchronous  (suboptimal  detector)  simulation  points  from 
Figure  4  are  used  here  for  reference.  The  simulation  model  was  coherent,  worst-case 
offset  for  both  rules  (3.39)  and  (3.41).  As  seen,  the  agreement  with  the  theory  is  quite 
satisfactory. 


Osimulation  results  for  synchronous  test  (3.8) 

□  simulation  results  for  asynchronous  test  (3.41) 

worst  offset 

^simulation  results  for  asynchronous  rule  (3.39) 

worst  offset 
(a)  N  =  50 

FigurcSa  Simulation  results  for  Pp  vs. 

asynchronous  detectors,  worst-case 
offset,  PpA  =  10”^ »  N=50. 


Osimulation  results  for  synchronous  test  (3.8) 
^simulation  results  for  asynchronous  test  (3.41) 
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Xsimulation  results  for  asynchronous  rule  (3.39) 

worst  offset 
(b)  N  =  1000 


Figure 8b  Simulation  results  for  vs. 

asynchronous  detectors,  worst-case 
offset,  P..  =  10"^.  N*1000. 


Thus  far  we  have  explored  structures  which  result  from  likelihood-ratio 
considerations  and  variations  thereof.  We  shall  now  take  a  look  at  some  other  detector 


structures  widely  used  in  practice.  Those  are  (a)  the  energy  detector  or  radiometer,  (b)  the 
product  or  crosscomelation  detector  from  two  independent  receptions  and  (c)  the  single-lop 
autocorrelation  or  chip-rate  detector. 

Although  all  three  are  basedband  detectors,  with  respect  to  the  carrier  harmonics,  they 
differ  in  that  the  first  two  utilize  the  DC  compont  of  the  output  as  a  decision  statistic  while 
the  third  filters  and  detects  at  the  chip  rate.  Another  grouping  yet  results  from  the  fact  that, 
although  all  three  rely  upon  a  good  measurement  of  the  noise  level  for  precise  threshold 
setting,  the  last  two  avoid  a  total-noise-power  factor  in  the  output  mean  thus  being  more 
robust  to  noise  fluctnations  than  the  first.  These  facts  are  all  borne  out  in  the  ensuing 
analysis.  We  should  clarify  here  that  the  term  "correlation"  above  is  not  to  be  intepreted  in 
the  statistical  sense;  rather,  it  signifies  a  product-and-filtering  (averaging)  mechanism  in  the 
time  domain.  Such  correlation  operations  comprise  standard  suboptimal^^  timing- 
recovery  methods  in  Pulse- Analog  (PAM)  and  Digital  Modulation  schemes  (see[17);  also, 
the  section  bit  synchronizers  in  [18]).  Furthermore,  the  correlation  detectors  presented 
herein  may  be  viewed  as  a  single-point  ("single-lag")  special  case  of  the  general  multi-lag 
correlation  algorithms  to  be  analyzed  in  Section  5.  Nonetheless,  because  of  their  popularity 
and  implementational  simplicity,  they  will  be  treated  here  separately.  Let  us  note  that  the 
analyses  and  conclusions  herein,  although  described  in  terms  of  DS  waveforms,  basically 
carry  over  to  TH  and  FH  waveforms  with  only  notational  changes.*^ 

i‘*Suboptimality  results  from  the  fact  that  the  received  noisy  waveform  is  not  correlated 
against  a  local  replica  of  the  useful  signal  or  an  estimated  version,  diereof,  but  rather 
against  itself  after  a  proper  delay. 

i^This  statement  neglects  second-order  effects  (for  low  input-SNR)  created  from  signal 
self-noise  and  signal-cross-noise  contributions. 


Detection  based  on  an  energy  measurement  is  an  optimal  LR  procedure,  when  both 
signal  and  noise  are  independent  Gaussian  processes  [12];  it  is,  therefore,  suboptimal  in 
principle,  when  the  sought  signal  is  not  Gaussian,  but  rather  a  BPSK  (or  any  other  type  of 
digitally  modulated)  waveform.  The  performance  of  the  radiometric  detector,  shown  in 
Figure  9,  has  been  analyzed.  It  can  be  shown  that 

-  (NTc)(o2  S6ki  ^  NoWbp)  (3.44a) 


and 


varal  Hkl  -  (NTo)(NoWbp  +  2oBpSNo6k,) 


(3.44b) 


where  the  attennation  factor  a|p  is  defined  as 


«BP 


/^jHL(f)  |2sc(f)dr  /2 


Sgtridf 


^2^  jwWbpTc/2  ^sin 


dx 


(3.44c) 


and  is  a  measure  of  the  energy  loss  due  to  filtering.  In  (3.44c),  HL(f)  stands  for  the  low 
pass  equivalent  of  the  BP  filter  HBp(f)  and  Sc(f)  is  the  power  spectral  density  of  the  code 
sequence  c(t).  Again,  assuming  that  Wbp  »  T,.*!  and  that  the  time-bandwidth  product 
Wbp  (NTj.)  =  N  is  large,  (3.44)  can  be  combined  with  the  "Gaussian"  performance  of 
(2.5).  The  result  is  (3.10a),  with 
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The  radio  a  =  a^gp/VWgpTj  can  be  plotted  versus  (W  BpTp),  and  the  optimal  choice  can  be 
shown  to  be  approximately  1,  from  which  it  follows  that  a^,*  *  0.77.  This 

corresponds  to  a  minimum  loss  of  -1.1  dB  with  respect  to  the  chip  noncoherent  detector 
(3.23).  Since  the  above  conclusions  arc  based  strictly  on  energy  considerations  and 
interchip-interference  effects,  due  to  filtering,  which  have  not  been  taken  into  account,  the 
actual  losses  are  somewhat  higher.  In  the  following  comparisons,  we  select  a  nominal  loss 
figure  of  - 1.5  dB.  We  note  that  the  above  losses  pertain  to  a  rectangular  BPF;  other  filter 
ships  will  result  in  different  loss  figures,  typically  of  the  same  order.  For  the  effect  of 
filtering  on  a  baseband  rectangular  pulse  (in  terms  of  output-SNR  losses)  with  respect  to  a 
matched  filter,  the  interested  reader  is  referred  to  [12]. 

We  now  turn  to  examine  whether  benefits  can  be  extracted  from  utilizing  two 
independent  receptions  of  the  same  signal,  resulting  from  two  RF  circuits  with  independent 
thermal  noises.  Although  only  the  synchronous  coherent  case  is  covered  here,  similar 
conclusions  are  expected  from  the  other  models. 

Let 

*'l(t)  -  /S”  c(t}di(^  *  n^ft);  0<tiT 

r2(t)  -  /S’  o(t'a)6|(i4'  njCt);  k«0,l  (3.46) 

represent  the  two  receptions.  Here,  the  two  noise  processes  njit)-,  m  =  1,2  are  Gaussian 

and  independent,  while  the  signal  componeik  is  the  same  in  both  waveforms. 
Furthermore,  the  second  signal  has  a  delay  of  A  seconds  with  respect  to  the  first  If  A 
were  exactly  known,  the  two  receptions  could  be  adjusted  to  align  in  time  and  be  coherently 
combined  with  a  resultant  3  dB  gain  in  SNR.  In  many  cases,  A  is  unknown  and  its 
estimation  is  of  interest  as  it  relates  to  path  differentials  and,  therefrom,  to  direction  finding. 


A  simple  cross  correlation  receiver  which  performs  this  estimation/detection  problem  is 
shown  in  Figure  10a.  The  corresponding  decision  rale  is 


H| 

X  -  ;NTc  (t)  Tj  (t)  dt  >  Xq 

Ho 


(3.47) 


To  simplify,  let  us  assume  that  detection  occurs  only  when  A  s  0;  then,  (3.47)  can  be 
analyzed  exactly  by  evaluating  the  second-order  statistics  of  the  product 
waveform  gi2(t)  =  ri(t)r2(t)  under  either  hypothesis.  Such  a  procedure  will  be  applied 

below  in  analyzing  the  chip-rate  detector.  Here,  for  the  sake  of  brevity,  let  us  just  consider 
the  alternative  structure  of  Figure  (10b)  which  can  be  thought  of  as  an  approximation  to 
Figure  (10a).  Obviously,  the  constraint  of  individual  chip-synchronism  per  channel  has 
been  imposed  on  the  second  structure,  in  addition  to  the  delay-sync  between  the  two 
channels.  The  rale  corresponding  to  Figure  (10b)  is 


where 


(3.48) 


r;"*  -  / 

‘  (J-I)T^ 


dt 


(3.49) 


are  synchronous,  chip-by-chip  integrations. 


Equation  (3.48)  is  analogous  to  (3.8)  once  rj  is  substituted  rj  rj.  .  Clearly,  in  the 

2 

absence  of  thermal  noise,  both  (3.8)  and  (3.48)  produce  the  same  quantity  NST^ ; 
however,  performance  is  different  in  noise.  Indeed,  from  (3.48),  it  easily  follows  that 


^{X|Hk}=N(NoTe)Ye5k, 


(3.50) 
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From  Antenna  2 

Figure  10b: 

An  approximate  version  of  the  crosscorrelation 
detector  in  Figure  lOa. 

var{XlHk}  =  N(NoTe)2  [1/4  +  yAil 


(3.51) 


The  resulting  performance  is  again  described  by  (3.10a),  where  now 

a  =2;  b  =  4  (Two  independent  coherent  receptions).  (3.52) 

Thus,  a  comparison  with  the  corresponding  performance  parameters  (3. 10b)  for  the  (one- 
reception)  synchronous  coherent  detector  reveals  a  gain  factor  of  •^,  or  1.5  dB  in  SNR  for 
the  present  system.  However,  this  gain  is  clearly  attained  at  the  cost  of  higher  complexity. 
Similar  conclusions  are  true  for  noncoherent  reception  detectors . 

Finally,  we  consider  the  chip-rate  detector,  a  noncoherent  BP  version  which  is 
shown  in  Figure  1 1.  This  is  an  one-antenna  system  Oil^e  the  radiometer)  which,  however, 
detects  power  at  the  first  harmonic  of  the  "signal  x  signal"  term.  Clearly,  this  structure 
could  also  detect  power  at  DC  and  serve  as  a  "generalized"  radiometer  where  the  A  =  0 
case  would  correspond  to  the  familiar  square-law  energy  detector  (Figure  9).  This 
possibility  is  explored  below,  where  it  is  shown  that,  aside  from  sensitivity  considerations, 
SNR  performance  increases  monotonically  as  A-»0.  Thus,  the  standard  radiometer  is 
probably  the  most  meritorious  device  to  consider  at  DC;  contrary  to  that,  detection  at  the 
first  harmonic  is  worth  exploring. 

We  note,  however,  that  the  robusmess  gains  to  be  shown  below  are  predicated 
upon  the  knowledge  of  the  sought  signal  rate,  a  requirement  not  present  in  the  radiometric 
analysis.  We  can,  therefore,  conclude  that  the  chip-rate  detector  trades  increased  parametric 
knowledge  with  a  certain  robustness  to  noise  fluctuations,  something  that  an  energy 
detector  is  not  equipped  to  do. 

Under  H|,  the  filtered  input  is  r(t),  given  by 


r 


c 


r(t)  =  v2aBpS  c(t)  cos(<fltot  +  ♦)  +  n(t)  (3.53) 

where  the  BP  nature  of  the  noise  process  is  reflected  in  its  PSD  8^(0  =  (No/2)  |HBp(f)|2 
and  the  standard  modeling  for  the  filtered  signal  via  the  factor  a^p  has  been  assumed. 

The  product  output  y(t)  =  r(t)r(t-A)  has  a  correlation  function 


Ry(t)  =  Rg  X  ^t)  +  Rs  X  n(‘')  +  «  N^"*) 

where  the  "signal  x  signal"  correlation  is 


(3.54) 
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4  eoa  auot 

(  eoa^wQ A  ♦  —————  )Rg(T)  (3.55a) 

In  (3.55a),  Rg(T)  signifies  the  corrlation  of  the  code-product  process  g(t)  =  c(t)  c(t-A). 
Furthermore,  the  "noise  x  noise"  correlation  is 


RnxN^i)  -  ^n(t)n(t-A)n(t-T)n(t-A-T)J 

-  ♦  &).Rn(t  -  t)  (3.55b) 

while  the  "signal  x  noise"  term  RsxnC"*)  includes  all  other  remaining  terms.  For  low  input- 

SNR,  the  signal  x  noise  term  will  be  neglected  in  comparison  to  the  noise  x  noise  term; 
thus,  the  PSD  Sy(f)  of  y(t)  can  be  approximately  written  as 

ii 
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where  Sg(f)  is  the  PSD  of  g(t). 

In  order  to  pursue  any  further  and  derive  the  second-order  statistics  of  the 
narrowband  output  X,  one  has  to  specify  Sg(0>  S^ff)  and  the  output-filter  characteristic. 
For  a  purely  random  code,  Sg(f)  was  derived  to  be 


where  the  offset  A  is  restricted  to  0  ^  |A|  <  T^.  Note,  the  presence  of  the  "signal  self-noise" 
term,  the  third  in  (3.57).  Because  of  the  weak  input-SNR  assumption,  this  term  will  be 
neglected  in  terms  of  its  contribution  to  the  total  noise  level.  However,  for  high  input-SNR 
environments  (as,  for  instance,  in  certain  FH  scenarios),  this  term  will  be  the  major  limiting 
factor  and  deserves  due  attention.  To  prove  further,  let  us  evaluate  the 
contribution  S^j^^  (f)  for  a  standard  "brick-wall"  noise  PSD 
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Sn(f)  - 


No  .  ,  ,  Wbp 

2  '  '  2 


. . . ,  otherwise 


(3.58) 


which  can  be  shown  to  be 

"  (NoWbp)2  s2(»HBpA)eos2wpaa( f )♦ 

“o 

—  (Wfip-  |f|)£t  4  cos(2«pa)  s,(2wa(ifBp-  (fplJ 
for  |f|<WBP  (3,59) 
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We  shall  restrict  attention  to  the  two  baseband  possibilities:  (a)  DC  Altering  by 
HdcCO  and  0>)  chip-rate  filtering  by  the  latter  corresponding  to  the  m  =  ±1  terms  in 

(3.57).  In  both  cases,  the  postdetection  filter  will  be  assumed  to  be  very  narrow  with 
respect  to  Rc  ;  for  instance,  if  HocCf)  corresponds  to  an  NT^  integrator,  this  implies  that 
N»l.  Note  that  Figure  1 1  refers  to  case  (b)  only.  Furthermore,  the  terms  around  2fo  will 
be  filtered  out  from  (3.56)  (they  have  already  been  excluded  from  (3.59)). 

For  a  narrowband  integrator-LPF(case(a)),  the  output  SNR 

<l2  ^  (/{X  |Hi  }  -  #{A  |Ho})2/varlA  |Ho  1 


is  easily  derived  from  (3.56)  -  (3.59)  to  be 
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(3.60) 


since  the  "noise  x  noise"  DC  contribution  cancels  out  in  the  different  -  ^Z/Hq). 

Certain  interesting  observations  can  be  made  from  (3.60),  especially  in  connection  with  the 
radiometric  output-SNR  d^^nj  from  (3.44);  if  we  form  the  ratio 
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then  it  is  immediately  clear  that  the  ratio  tends  to  1  as  A-^0,  as  it  should.  Furthermore,  if 
A  «  (2Wbp)-J  so  that  the  Sg(*)  function  in  the  denominator  of  (3.61)  is  approximately 
one,  the  ratio  is  roughly  equal  to  (1  -  |A/Tcl)2,  independent  of  the  signal  carrier  frequency 
fg.  Thus,  in  this  range  of  A,  the  ratio  is  maximized  by  choosing  A  =  0,  i.e.,  square-law 

energy  detection  yields  the  highest  output  SNR.  The  same  conclusion  holds  far  higher 
values  of  A,  say,  A  ^  the  Sg(»)  function  diminishes  and  the  ratio  tends  to 

2(l-|Arrel)2cos2(a)oA) 

However,  since,  typical  values  for  the  filter  BW  are  W3p>Rp,  it  follows  that  A  k  Tg/2, 
which  implies  that  the  ratio  is  upperbounded  by  (1/2)  cos^cOgA  ^  1/2  or  a  3  dB  loss  of  the 

correlator  versus  the  radiometer.  In  addition,  one  has  to  worry  in  the  latter  case  about  the 
exact  relationship  between  fg  and  A.  If  f^A  #  k/2,  k  an  integer,  servere  loss  of  "signal" 

power  can  occur  due  to  the  cosine  term.  Clearly,  in  an  LPI  environment  where  the  signal 
parameters  are  hardly  at  the  disposal  of  the  interceptor,  such  limitation  is  of  concern.  On 
the  other  hand,  the  mean  "noise  x  noise"  contribution  at  DC  (the  first  term  in  (3.59))  also 
diminishes  with  increasing  values  of  (WgpA),  a  factor  which,  although  not  directly 

reflected  in  the  SNR  computations,  could  affect  the  sensitivity  of  the  system  in  an  actual 
implementation  with  imperfect  AGC's. 

Let  us  now  turn  to  the  chip-rate  filter  (case(b)).  Evaluating  the  S  x  S  and  N  x  N 
narrowband  terms  at  f  =  R^,  we  get  the  following  SNR  expression. 

OBpS^  (1  ♦  cos  2ta>oA)f(Te)  ^a(~T“ 


No  NbpJ  (1-p)Lt^009(2»QA)S,(2»aWBp(1“p) )) 


‘.’.•‘ry-r: 
■••.nVA 


(chip-rate  detector) 


(3.62) 


where  Wj  is  the  BW  of  the  narrrowband  filter  H  (f)  and  the  ratio  p  =  Rg/W gp  is  assumed 
to  be  in  the  range  0  <p<  1.  Note  that  if  Wj  *  (NTj)*^,  then  the  first  bracketed  term 
corresponds  to  the  radiometric  SNR;  thus,  the  fraction  in  the  second  brackets  represents  an 
SNR-ratio  similar  to  (3.61).  It  is  further  shown  that  the  numerator  of  this  ratio  (i.e., 
the  S  X  S  portion)  is  maximized  for  A  =Tg/2,  assuming  that  the  condition  f^A  =  k/2; 

k  =  integer,  is  satisfied.  For  that  "optimal"  choice  of  the  offset,  the  resulting  SNR  is 

"I  -  ('‘Lra.)(-^ 

(chip-rate  detector) 

(a  -  Tq/2;  TqA  -  k/R)  (3.63) 

Note  that  the  SNR  of  the  chip-rate  detector  increases  as  p-»l,  i.e.,  as  Wgp-^Rp-,  in  fact, 
for  Wgp  =>1.25  Rg,  the  two  detectors  (chip-rate  and  radiometric)  are  roughly  equivalent. 
Since  much  smaller  values  of  Wgp  arc  rather  inappropriate  to  consider  in  view  of  Doppler 
uncertainties,  signal  filtering,  etc.,  such  equivalence  can  be  considered  a  practical  fact. 
Again,  the  aforementioned  tradeoffs  between  robustness  to  noise  fluctuations,  knowledge 
of  chip-rate  and  adjustment  of  carrier  frequency  should  be  kept  in  mind  for  a  true 
comparative  assessment  of  the  various  detectors. 

As  mentioned,  a  complete  analysis  for  all  values  of  the  input  SNR  should 
incorporate  the  signal  self-noise  contribution  from  (3.57),  as  well  as  the  signal  x  noise 
terms. 

Furthermore,  the  dependence  on  the  carrier  frequency  can  be  totally  eliminated  by 
using  the  alternative  structure;  the  penalty  there,  of  course,  will  be  a  somewhat  lower 
output  SNR.  Finally,  let  us  note  that  variations  of  the  schemes  herein  are  also  possible. 
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involving  different  intermediate  frequencies  (IF)  for  the  undelayed  and  delayed  versions, 
respectively.  We  shall  not  pursue  the  topic  here  any  further. 


3.5  Comparisons  and  Discussion 


The  perfonnance  of  any  of  the  schemes  discussed,  so  far,  will  always  be  upper- 
bounded  by  that  of  a  perfectly  code-matched  system  (i.e.,  one  where  the  code  is  totally 
known)  and  lower-bounded  by  the  performance  of  the  radiometer.  The  distance  d^g 

corresponding  to  the  (unattainable  from  an  LPI  viewpoint)  upper  bound  is  easily  shown  to 
be 

due  =  (3.64) 

where  a  coherent  system  has  been  assumed.  Comparing  (3.64)  with  (3.1 1)  implies 
that  the  lack  of  code  knowledge  amounts  to  an  SNR  loss  of  between  two  synchronous 

coherent  detectors,  one  equipped  with  the  code  and  the  other  not;  in  dB,  this  means  that 
knowledge  of  the  code  reduces  in  half  the  input  SNR  (dB)  required  to  achieve  a  certain 
performance  level.  This  is  depicted  in  Figure  12  (compare  curves  and  (^)  along  with 
the  performance  of  the  other  alternative  schemes.  We  observe  that  the  difference  between 
@and(^is  roughly  1.5  dB,  which  is  also  the  difference  between ^and(^ as  well  as(^ 
and@.  It  follows  that,  under  fairly  ideal  conditions,  i.c.,  two  independent  synchronous 
coherent  receptions,  one  could  utilize  up  to  4.5  dB  gain  above  the  radiometer;  each  of  the 
added  features  (i.e.,  two  receptions,  synchronism  anc  ■;oherence)  can  be  thought  of  as 
contributing  1.5  dB  to  the  gain.  We  note  that,  in  arriving  at  these  gain  figures,  system 
imperfections  and  noncalibrated  parameters  such  as  Doppler  frequency  offset  have  not  been 
accounted  for. 

It  is,  therefore,  seen  that  certain  gains  are  plausible  with  respect  to  the  radiometer  if 
careful  designs  are  employed,  but  they  are  certainly  not  overwhelming-at  least  in  the  SNR 

17  For  example,  if  the  unknown-code  detector  requires  -15  dB  for  certain  perfonnance  level, 
the  known-co^  system  will  achieve  the  same  level  at  -  30  dB 


1.  Completely  known  waveform 

2.  Two  independent  receptions,  synchronous  coherent  detector 

3.  Synchronous  coherent  energy  detector 

4.  Synchronous  chip-noncoherent  energy  detector 

5.  Radiometer 


Figure  12.  Performance  Comparisons  for  the  Detection  of  DS  Waveforms 


range  of  interest  and  for  constant-envelope  signals.  This  is  because  a  high-rate  DS 
waveform  is  not  veiy  distinguishable  from  the  background  thermal  AWG  noise  in  which  it 
is  detected  and,  as  is  well  known,  the  radiometer  then  becomes  an  asymptotically  optimal 
detector.  This  is  also  in  accordance  with  Krasner's  approximate  result  regarding 
biorthogonal  waveforms  [8],  which  establishes  the  DS  signals  as  uniformly  good  LPI 
choices  almost  independent  of  the  prevailing  system  parameters  (SNR  and  N).  However, 
if  one  removes  the  biorthogonality  of  the  signal  set  and/or  the  AWG  nature  of  the  noise,  the 
results  can  differ  substantially.  An  illustration  of  the  former  can  be  found  in  the  following 
sections  which  pertain  to  time-hopping  and  frequency-hopping  waveforms.  As  for  the 
latter,  it  is  clear  that  the  white-noise  theory  as  developed  herein  can  be  grossly  inadequate 
in  a  dense  environment,  loaded  with  a  mixture  of  wideband  and  narrowband  waveforms. 
As  mentioned  in  the  introduction,  the  presence  of  narrowband  interference  prediction  and 
excision  algorithms  will  then  become  crucial  to  the  success  of  any  detection  scheme.  The 
example  of  section  5.5  is  meant  to  illustrate  the  dramatic  difference  between  a  somewhat 
"intelligent"  algorithm  and  the  "naive"  radiometer  in  such  a  narrowband  environment 
which,  of  course,  comes  as  no  surprise.  We  note  that  a  general  theory  of  LPI  performance 
in  an  unpredictable-interference  scenario  has  not  yet  matured  (at  least  in  the  open  literature) 
and  thus  constitutes  an  area  of  significant  interest 

We  conclude  by  noting  that  the  implicit  assumption  of  "message-synch"  (i.c., 
signal  occupying  exactly  the  observation  interval)  should  and  could  be  removed  in  a 
realistic  evaluation  of  an  LPI  detecting  system..  However,  the  relative  ranking,  in  terms  of 
merit  of  the  various  detectors  discussed  herein  is  not  expected  to  change  under  such 
alteration  of  the  underlying  scenario. 
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4.0  TIME-HOPPING  WAVEFORMS 


We  now  focus  on  the  wideband  detection  of  a  particular  type  of  orthogonal 
signalling,  namely  TH  waveforms,  where  the  code  dictating  the  hopping  is  again  random. 
We  shall  consider  only  synchronous  detectors  (i.e.,  the  timing  or  epoch  of  the  hopping  slot 
will  be  assumed  known);  asynchronism  can  be  treated  with  methods  similar  to  those  used 
in  Section  3.3.  Furthermore,  both  carrier-coherent  and  carrier-noncoherent  systems  will  be 
examined.  In  particular,  the  former  will  assume  a  known  carrier  phase  combined  with  a 
pulse-position-modulation  (PPM)  format  while,  for  the  latter,  any  modulation  that 
randomizes  the  phase  form  frame  to  frame  (e.g.,  BPSK)  is  well  suited.  We  term  such 
cases  "pulse  noncoherent".  Of  course,  one  could  hypothesize  a  noncoherent  PPM  system 
with  an  unknown,  but  constant  carrier  phase,  and  derive  the  optimal  detection  rule. 
However,  the  resulting  receiver  and  pertinent  analysis  are  again  complicated  (for  a  similar 
situation,  refer  to  section  3.2  for  DS);  therefore,  only  pulse-noncoherent  systems  will  be 
discussed  here.  Again,  the  lower-bounding  performance  of  the  radiometer  will  indicate  that 
little  is  lost  by  using  this  simplification. 

For  our  purposes,  the  received  waveforms  (TH  signal  plus  noise)  can  be  written  as 

00 

r(t)  =  /2S  p(t  -  kTp  -  PkTH)cos(u,ot  +  e,^)  +  n(t)  (4.1) 

where  Tp  is  the  frame  length,  Th  is  the  hop  length  (width  of  each  time  slot),  and  p^,  0^  are 

IS 

random  variables  which  are  independent  from  frame  to  frame  and  of  each  other,  denoting 
the  slot  location  and  carrier  phase,  respectively,  during  the  kth  frame.  Here,  pj^  can  take 
on  any  one  of  the  equiprobable  values  p^  =  0,...,  Np  -1,  where  Np  =  T/T^  is  the  total 
number  of  slots  per  frame  and  6^  summarizes  both  the  unmodulated  carrier  phase  <|>  plus 

**  The  independence  of  0^  is  within  the  aforementioned  spirit  of  a  "pulse-noncoherent" 
system. 
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any  superimposed  PSK  modulation.  For  coherent  systems,  6)^  is  assumed  known  for 
every  k.  Finally,  p(t)  is  a  unit  pulse  of  duration  T|j  seconds,  while  n(t)  is  the  usual 
bandpass  AWGN.  A  noiseless  sample  waveform  (realization)  of  the  signal  is  shown  in 

Figure  13.  We  shall  assume  that  the  total  observtation  interval  consists  of  Q  frames, 
ie.,  T  =  QTp  =  QNpTn.  Note  that  the  average  signal  power  in  (4.1)  is  S,v  =  S/Np. 
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4.1  Synchronous  Coherent  Detectors 


Under  the  coherent  assumption,  VF cos  (  ©ot  +  8^)  in  (4,1)  is  substituted  by  one 
Then,  a  straightforward  application  of  the  composite  lilcelihood  ratio  yields 


ni=l  )  k=l  ( 


where 


|lir  kmu 


A  STh 
"H  ■ 


hop  (slot)  and 

/(in-l)Tp+kT„ 

^  «  r(t)  pCt  -  kTfOdt 

(in-l)Tp+(k-l)Tj^ 


is  the  integral  of  r(t)  in  the  kth  slot  of  the  mth  frame.  If  is  a  small  number  so  that 
exp{x}  •  1  +  x,(4.2)  can  be  simplified  to  the  approximate  rule 


/*T 

ACKt)]  -  /  r(t)dt  ^  Xq 


i.e.,  just  the  coherent  integral  of  r(t).  However,  since  Yh  is  typically  above  0  dB  (much 
larger  than  Yc)i  such  approximations  arc  not  as  successful  as  in  the  DS  case  -  a  fact  to 

which  the  comparisons  in  Section  4.3  will  attest 

The  Gaussian  distance  d^g^h  associated  with  the  optimal  coherent  rule  (4.2)  can  be 

derived  based  on  (2.9)  and  the  following  steps: 


I 


where  Ap  is  the  likelihood  ratio  per  frame.  But, 


Since  rjyn  is  a  zero-mean  Gaussian  random  variable  (rv)  under  Hq.  Combining  (4.6)  and 
(4.7)  results  in 


"coh  =  l«[i  -  l)]  c-*) 


The  distance  d  %  for  approximate  rule  (4.6)  is  derived  much  easier  since 

ff{XlHo}  -  0.  r{XlHi}  -  and  var{X|Ho}  “  (NpQTn)No/2. 

Therefore, 


‘^X.coh 


-  WxIHq) 
(varlxlHg))^^^ 


(4.9) 


It  is  clear  from  (4.8)  and  (4.9)  that  d^  ^  asymptotically  approaches  d  ^  .coh  “ 
an  expected  result  Furthermore,  we  can  compare  the  approximate  coherent  rules  ( 3.8) 
for  DS  versus  (4.6)  for  TH  and  note  that  the  first  involves  a  nonlinear  operation  on  the  data 
prior  to  integration,  whUe  the  second  does  not  This  has  a  reflection  on  their  corresponding 
distances,  as  (3.11)  is  directly  proportional  to  Yc  while  (4.9)  is  proportional  to  Yh. 

for  the  very-low-SNR  case,  the  latter  would  outperfom  the  former. 


Starting  from  (4.1)  and  invoking  the  independent  phase  assumption  results  in  the 


decision  rule 
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where  R|^  is  the  kth  slot,  mu,  frame  envelope 
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=  2  y  ^  ^  r(t) 

cos  Upt 

_Q,km 

(m-l)Tp+(k-l)T„ 

sin  u)Qt 

dt 


(4.11a) 


(4.11b) 


The  distance  di_„o„coh  associated  with(4.l0)is  derived  by  steps  identical  to  (4.7)  as 


'^l.noncoh  “  “  Ol  (4.12) 

Since  Io(x)  <  cxp{x}  for  every  x  >  0,  it  follows  from  the  comparison  of  (4.8)  and  (4.12) 
that  d|^  coh  >  noncoh  for  every  Yh,  as  expected.  We  note  further  that,  unlike  the  coherent 
distance  (4.8)  or  its  approximate  counteipart  (4.9).  dj^  noncoh®^  (4.12)  is  proportional  to 
at  low  SNR. 

Instead  of  analyzing  the  noncoherent  rule  which  approximates  (4.10)  for  low  SNR,  let 
us  just  consider  the  performance  of  the  radiometer  which  operates  on  the  time-hopped 


waveform  (4.1)  for  T  seconds.  Following  the  steps  outlined  in  Section  3.4  it  can  be 
be  shown  that 


^rad  *  ^/^rad^^H 

(4.13) 

where  the  constant  =  o4/(WbpTh)  again  measures  the  loss  due  to  filtering.  As  was 

done  in  Section  3.4,  can  be  set  at  its  optimal  value  =  (0.77)2  «  o.6.  We  should 

note  here  that  (4.13)  is  actually  an  optimistic  prediction  (i.e.,  an  upper  bound)  of  the 
radiometric  performance  since  it  is  based  on  the  assumption  that  the  variance  of  the  test 
statistic  under  Hj  is  approximately  the  same  as  under  Hq.  For  medium  to  high  input  SNR, 

this  is  not  true;  although  a  more  meaningful  performance  description  can  easily  be  derived  , 
we  shall  be  content  with  (4.13)  for  comparison  purposes. 


4.3  Comparisons  and  Discussion 


First,  we  are  interested  in  comparing  the  relative  loss  of  the  approximate  coherent 
rule  (4.5)  versus  the  exact  (4.2),  since  the  former  can  be  implemented  easily  (a  simple 
integrator)  in  contrast  to  the  significant  complexity  of  the  latter.  The  comparison  is  in  terms 
of  the  relative  SNR  values  Yh  dB)  required  by  each  in  order  to  achieve  the  same 
performance  level  (d^^)  for  the  same  fixed  values  of  Q  and  Np.  It  is  based  on  (4.8)  and 
(4.9)  and  is  shown  in  Figure  14  (Np  =  10).  As  expected,  the  difference  diminishes  at  low 
SNR,  but  is  rather  pronounced  at  medium  to  high  SNR.  So,  Yfoh,opt  =  0  dB  corresponds 
to  Y^oh,appr  _  3  9  dB  —  a  gap  that  increases  rapidly  as  opt  increases. 

Second  we  look  at  the  gains  brought  about  by  the  coherent  assumption  versus  the 
noncoherent  by  comparing  the  Gaussian  distances  in  (4.8)  and  (4.12) .  The  result,  again 
in  terms  of  the  required  SNR's  for  the  same  performance  level,  is  shown  in  Figure  15. 
(Note  that  this  comparison  is  independent  of  N  pj.  The  difference  here  diminishes  as  SNR 

increases,  while  it  can  be  substantial  at  low  SNR.  The  final  comparison  is  between  the 
noncoherent  optimal  performance  (4. 12)  (or  its  equivalent  log-likelihood  sum)  versus  the 
radiometer  performance  estimate  (4. 13).  It  is  shown  in  Figure  16  for  =  0.6  and 

Np  =  10.  The  irreducible  distance  of  1.1  dB,  as  the  SNR  goes  to  zero,  is  due  to  the/K^gj 
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factor  .  Again,  we  should  keep  in  mind  that  the  actual  SNR  losses  of  the  radiometer  are 
higher  than  those  shown  in  Figure  16  by  an  amount  that  increases  with  SNR  due  to 
previously  discussed  reasons.  Furthermore,  a  common  trend  is  evident  from  Figures  14 
and  16,  namely,  that  the  simplifying  deviations  from  the  optimal  decision  rules  incur 
comparative  losses  that  increase  fairly  rapidly  with  the  available  hop  SNR  for  values  of 
Yh  above  2-3  dB.  Since  this  is  the  dominant  range  of  importance  in  TH  applications, 
optimal  devices  (albeit  complex)  should  attract  due  attention. 

19 

However,  this  margin  will  probably  disappear  if  the  loss  due  to  asynchronism  is 
accounted  for  in  an  asynchronous,  noncoherent,  optimal  detector. 


So  far  a  multitude  of  optimal  and  subopdma]  receiver  structures  for  the  wideband 
detection  of  DS  and  TH  waveforms  in  AWGN  were  derived  and  evaluated  in  tfiis  report 
The  starting  points  were  those  receivers  resulting  from  the  assumed  knowledge  of  many 
pertinent  parameters  (except  the  codes  themselves);  upon  gradual  removal  of  these 
assumptions,  a  family  of  suboptimal  structures  emerge  whose  associated  losses  furthered 
them  from  the  globally  optimal  performance  towards  the  simple  lower-bounding 
performance  of  the  radiometer.  It  is  generally  concluded  that  the  higher  degree  of 
sophistication  and  complexity  for  such  detectors  is  more  justifiable  in  the  TH  case  than  in 
DS,  both  because  of  the  pirevailing  SNR  values  and  the  signal  structure  of  TH  versus  DS. 


5.0  FREQUENCY-HOPPING  SIGNALS 


For  FH  systems,  the  main  challenge  in  implementing  optimal  receivers  is  the  large 
number  of  hopping  frequencies  (as  many  as  available  by  the  FH  system)  that  must  be 
simultaneously  observed  and  measured  per  hop  before  an  optimal  decision  is  reached. 
Contrary  to  that,  the  two  systems  considered  previously  utilize  a  single  carrier  frequency; 
stated  differently,  the  instantaneous  bandwidth^Oof  those  waveforms  effectively  comprises 
the  overall  spread-spectrum  bandwidth.  The  implication  is  that  any  additional  beneEt 
(always  with  respect  to  the  radiometer)  brought  about  by  the  optimal  detectors  in  the  DS 
and  TH  cases  must  be  extracted  from  the  intelligent  time-domain  processing  of  the  received 
waveform,  rather  than  from  increasing  hardware  complexity.  In  FH/LPI,  however,  we 
shall  see  that  detection  performance  improves  as  the  channelized  width  of  the  spectral 
observables  decreases  towards  the  optimal  hopping-rate  width  which,  of  course,  implies 
more  filters. 

The  FH/LPI  problems  are  mathematically  formulated  as  detection  of  a  complex 
sinusoid  of  unknown  frequency  in  additive  white  Gaussian  noise  (AWGN).  This  is  a 
problem  common  to  diverse  fields  such  as  Doppler  radar  and  sonar  [22-26],  spread 
spectrum  communications  [27-33],  robust  detection  in  the  presence  of  signal  perturbations 
("slewing")[34],  unknown  colored  noise  [35],  unknown  signal  spectrum  [36-37], 
unknown  bandwidth  due  to  digital  modulation  [38]  (see  also  [39]  for  some  simulation 
results),  etc.  Unfortunately,  the  wide  diversity  of  the  underlying  models  precludes  the 
extraction  of  firm,  universal  conclusions  about  the  relative  merit  of  the  different  techniques, 
either  among  themselves  or  with  respect  to  the  simplest  hypothesis-discriminating  device, 
namely  the  energy  detector  (radiometer).  Here,  we  shall  address  the  problem  from  a  basic 

^^By  this,  we  mean  the  bandwidth  associated  with  the  single  chip  pulse  of  the  spreading 
code. 


communications  viewpoint,  cautiously  establishing  the  discrete-space  detection  model,  and 
we  shall  subsequently  explore  a  variety  of  possible  approaches  to  its  solution. 


Let  s(t)  be  the  FH  signal  to  be  detected  solely  in  AWGN  within  a  total  spread 
bandwidth  of  Hz.  If  Rh  L  Th'^  is  the  hopping  rate  and  hopping  frequencies  are 

contiguous  and  cquispaced  by  the  minimum  noncoherent  orthogonal  separation  of  f  =  Rh. 
it  follows  that  the  total  number  M  of  possible  hopping  locations  is  M  =  W,  /  =W,Tj£ 

-- typically,  a  very  large  number.  Let  the  total  observation  time  T  be  Ny  hops,  i.e., 

T  =  NtTh. 

We  restrict  attention  here  to  the  broadband,  AWGN  case,  which  is  the  only 
observable  waveform  under  hypothesis  Hq: 


r(t)  =  n(t);  0  ^  t  <  T 

(Ho) 

(5.1a) 

If,  under  the  alternative  hypothesis  Hi, 

r(t)  =  s(t)  +  n(t);  0  t  <  T 

(Hi) 

(5.1b) 

the  signal  s(t)  is  itself  broadband  C'white")  and  Gaussian,  the  energy  detector  is  tantamount 
to  the  optimal  likelihood  ratio  (LR)  test 


f(r(t)|H,) 

"o 

When,  however,  the  signal  of  interest  is  inherently  narrowband,  as  is  the  case  of  interest 
here,  it  possesses  sufficient  structure  so  that  parametric  approaches  warrant  attention, 
whether  optimal  and  theoretically  justifiable  or  not. 

Implementation  of  (5.2),  or  any  variant  thereof,  over  the  total  RF  bandwidth  W^ 

might  result  in  unacceptable  complexity.  One  way  to  alleviate  the  problem  is  shown  in 
Figure  17  wherein  the  total  bandwidth  W,  is  subdivided  into  large  contiguous  segments  of 


B  Hz,  each  being  much  larger  than  the  optimal  bandwidth  of  Rh  Hz.  Appropriate 

processing  produces  a  per-band  decision  which  is  then  fed  into  an  overall  accumulator 
(e.g.,  a  majority-logic  combiner)  for  a  final  decision  (H^  versus  Hq).  Let  us  note  that 

alternate  reduced  configurations  have  also  appeared  in  the  literature  [40,41],  such  as  the 
partial-band  filter-bank  combiner.  There,  only  a  fraction  of  the  total  number  M  of  the 
hopping  slots  is  being  observed,  but  the  observadon  bandwidth  per  slot  is  optimal  (B=Rh). 

The  focus  of  the  following  sections  is  on  the  waveform-processing  (WP)  aspect 

on  a  per-band  basis.  Clearly,  algorithms  that  improve  performance  on  tiiat  level  will  also 
increase  the  overall  system  performance.  We  shall  assume  throughout  this  development 
that  the  hop  time-bandwidth  product 

G  =  B  Tjj  =  B/Rf{ »  1  (5.3) 

is  very  large,  i.e.,  of  the  order  of  hundreds  or  more.  The  factor  G  also  measures  the  order 
of  the  reduction  in  complexity  when  implementing  the  suboptimal  structure  of  Figure  17 
instead  of  the  optimal.  As  a  consequence,  the  input  SNR  per  band  ,  defined  as 


is  typically  very  small  compared  to  unity  (75^, «  1). 

Now,  two  questions  arise  immediately: 

(a)  Into  how  many  bands  should  the  spectrum  be  divided? 

(b)  What  kinds  of  WP  alternatives  are  available  ot  could  be  used? 

With  respect  to  (a),  it  is  intuitively  clear  that  one  should  use  as  many  bands  as  the  number 
of  devices  which  can  be  afforded  since,  by  decreasing  B  and  thus  increasing  from 

(5.4),  more  reliable  decisions  on  a  per-band  bases  can  be  achieved;  an  elaborate  cost- 
versus-payoff  study  is  of  significant  interest  in  this  area,  but  will  not  be  pursued  here.  As 
for  the  second  question,  we  note  that  the  device  which  immediately  comes  to  mind  is  the 
familiar  radiometer,  designed  to  measure  energy  in  the  band  of  B  Hz.  Under  the  Gaussian 
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assumption  for  the  output-decision  statistic  (which  is  well  justified  for  very  large  G  via  a 
central-limit-type  argument),  its  performance  is  easily  derived  to  be  (see  Sections  3.4, 4.3) 


out 


(5.5a) 


SNR„„t  =  SNrS  =  (Nt)  »  (NtG)-»  Y^h 


(5.5b) 


(PdjPfa)  are  the  detection  and  false  alarm  probabilities,  respectively,  and  Q'Kx)  is  the 
inverse  of  the  Gaussian  integral  function,  and 


No 


(5.5c) 


is  the  hop  SNR.  Note,  that  the  second  equality  in  (5.5b)  signifies  the  small-signal 
suppression  effect  of  square-law  detectors  [42,  page  267].  We  emphasize  again  the  (Pj), 
Pfa)  pertain  to  the  per-band  decision,  not  the  overall  scheme. 

The  radiometer  solution  is  implementationally  simple,  but  otherwise  unacceptable  in 
most  cases  from  a  performance  viewpoint;  we  shall  soon  have  a  chance  to  establish  this 
latter  fact  Alternatively,  more  sophisticated  approaches  must  be  sought,  and  that  is  what  the 
following  sections  are  about,  namely  to  establish  and  analyze  families  of  advanced 
detection  rules  (waveform-processing  schemes). 

A  major  aspect  of  this  detection  problem  is  to  choose  the  domain  in  which  the  decision 
algorithm  is  to  operate,  namely  the  spectral  versus  the  correlation  domain.  Of  course,  by 
virtue  of  the  Wiener-Kinchine  theorem  the  two  domains  contain,  in  principle,  the  same 
amount  of  second-order  information.  However,  how  this  information  can  be  best  exploited 
from  a  detection  viewpoint,  when  a  finite  length  data  record  is  available  and  certain 
parameters  in  the  signal  and/or  noise  models  are  not  completely  known  (in  either  a 
statistical  or  a  deterministic  sense),  remains  still  a  challenging  question.  It  has  prompted 
researchers  to  look  beyond  the  traditional  spectral  approach  (which  is  motivated  mostly 


from  likelihood  ratio  considerations)  into  the  correlation  domain  [43, 44, 26, 30-38];  the 
results  reported  exhibited  a  varying  degree  of  success,  depending  on  the  assumptions  made 
and  the  nature  of  the  adopted  decision  rule.  An  added  degree  of  motivation  for  the 
correlation  domain  comes  from  the  recent  advances  in  the  technology  of  real-time,  large 
time-bandwidth-product  autocorrelation  devices  (see  [45]  and  refs,  therein). 


Now  we  shall  explore  single-hop  algorithms  (Nj  =  1);  in  other  words,  the  data  record 
consists  of  single-hop  0  ^  t  <  of  a  perfectly  epoch-synchronized.signal.  The  results 
should  give  an  indication  of  the  relative  merit  of  the  different  algorithms,  their  capabilities, 
complexity  limitations,  etc.  Futhermore,  we  shall  assume  throughout  this  section  a 
noncoherent  environment,  which  is  mostly  dictated  by  practical  considerations  of  the  FH 
modulation. 

Let  us  first  consider  algorithms  in  the  spectral-domain.  In  order  to  motivate  them,  let 
us  assume  that  the  signal  (under  H^)  is  just  a  sinusoid  of  known  power  S  but  unknown 
frequency  fj  and  phase  0;  frequency  fg  lies  in  an  RF  bandwidth  of  B  Hz,  centered  at  the 
known  frequency  fg,  i.e. 

r(t)^cos(ci)st  +  e)  +  n(t)  (Hj)  (5.6) 


where  [fg-fd  ^  B/2  and  Th  »  i.e.  G  =  BTh  »  1.  This  situation  is  depicted  in  Figure  18. 
The  one-sided  noise  power  spectral  density  (PSD)  Nq  W/Hz  is  also  known.  Under 
the  broadband  asssumption  for  the  noise  and  modeling  the  unknown  frequency  fg  as  a 
random  variable,  uniformly  distributed  in  (fc-B/2,  f^.  +  B/2),  leads  to  the  optimal  LR  test 
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(5.7a) 


where  !{>(•)  is  the  zero*  -order  Bessel  function,  stands  for  expectation  and  R(f,)  is  the 
real  envelope  at  frequency  fg. 
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Figure  18  A  pictorial  representation  of  the  RF  spectrum 


around  center  frequency  f_  under  H,  (signal  present) 


R(g=  VR,2(Q  +  RQ2(g 


(5.7b) 


with 


Ri(fs)=  V2  Jor(t)  cos  27tf,t  dt 

(5.7c) 

RQ(fs)  =  r(t)  sin  27tfst  dt 

(5.7d) 

the  inphase  and  quadrature  components,  respectively.  This  optimal  nile(5.7a)  averages  the 
nonlinear  transformations  of  the  envelopes  at  all  frequencies  in  the  band.  This  being  overly 
demanding  to  implement,  one  discretizes  the  frequency  uncertainty  region  in  G  =  B/(T'*) 
candidates  f^;  m  =  0,...,G  - 1,  spaced  T^^  Hz  apart,  and  approximates  (5.7a)  as 
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threshold 


(5.8) 


with  R(fn,)  as  per  (5.7).  The  set  of  2G  =  2BTh  orthogonal  function  {sin  a)mt,cos  ©„,t}G“* 

1D»0 

can  be  approximately  interpreted  as  a  predominant  set  of  eigenfunctions  of  the  broadband 
noise  with  a  rectangular  spectrum  [12].  Note,  that  if  the  candidate  hopping  frequencies 
were  known  exactly  a_priori.  rule  (5.8)  could  be  derived  directly  from  rule  (4.9),  (with 
Q  =  1)  by  exploiting  the  duality  between  TH  and  FH.  In  fact,  both  spreading  techniques 
belong  to  a  larger  class  of  orthogonal  signal  sets. 

From  a  theoretical  standpoint,  any  mechanism  which  provides  reliable  envelope 
esdmates^i  |R„,|  =  R(fn,)  or,  equivalently,  spectral  estimates  IRmP;  m  =  0,...,G-1  can 

21  Since  both  real  and  complex  quantities  will  be  encountered  in  the  sequel,  the  symbol  "n" 
will  denote  a  complex  number,  while  |*|  will  stand  for  the  complex  norm;  Plain  letters  refer 
to  the  real  counterparts. 


be  coupled  with  the  aforementioned  LR-type  test  or  any  other  approximations  thereof. 
Possible  choices  include  the  aforementioned  quadrature  analog  implementation  (5.7b)- 
(5.7d)  a  bank  of  narrowband  filters,  chirp-transform  (SAWD)  devices  [45]  or  a  complex 
FFT  operating  on  the  baseband  complex  envelope  r(t)  of  the  real  bandpass  observation  r(t): 


rtt)  =  r(t)  exp( 

(5.9a) 

ri<  -  r(kB~''  )  ;  k-0 . G-1 

(5.9b) 

,  FFT{rk}°~’ 

(5.9c) 
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-  Pm  P:  m-0,  ...  ,G-1 

(5.9d) 

As  indicated,  this  requires  sampling  the  baseband  complex  signal  7(0  at  the  Nyquist  rate  of 
B  complex  samples/sec,  and  is  based  on  the  fact  that  r(t)  and  r(t)  possess  identical  spectral 
envelopes,  as  per  (5.9d).  The  information  loss  incurred  in  the  sampling  operation  (5.9b)  is 
negligible  for  large  G;  furthermore,  denser  spectral  sampling  can  be  obtained  by  padding 
{r^}  with  an  appropriate  number  of  zeros,  although  such  is  not  required  for  implementing 

(5.8).  Issues  associated  with  the  standard  generation  of  the  "periodogram”  (5.9d)  such  as 
data  windowing,  segmentation  (Welch  method),  etc.  have  been  studied  extensively  in  the 
context  of  spectral  estimation  (see,  for  instance,  [46]  and  refs  therein)  and  will  not  be 
discussed  here  any  further. 

Let  us  now  return  to  the  performance  of  rule  (5.8).  It  is  clear  from  the  above  problem 
formulation  that,  for  G  large,  the  LR  statistic,  A  of  (5.8)  should  be  modeled  as 
Gaussian  (under  a  central-limit-theorem  type  of  argument)  which,  in  turn,  negates  the 


modeling  assumption  of  a  Gaussian  LLR  i .  As  a  result,  the  more  traditional  measure  d^  is 
not  applicable;  instead  a  "Gaussian”  measure  dy^  must  be  derived  and  employed. 

Let  us  define  dy^  via  the  standard  expression  (2.7) 

d'  k  - 

var{A|Ho} 

Where  S{A\  HJ  is  the  expectation  of  A  under  hypothesis  i  =  0,1  and  var  {A|  Hq}  is 

its  variance  under  Hq.  Let  us  further  note  that,  under  Hq,  all  G  observables 

just  noisy  envelopes,  while  under  Hj  there  is  one  "signal-plus-noise"  envelope  and  (G-1) 

noisy  ones  (recall  the  orthogonality  assumption).  Thus,  for  large  G,  it  is  reasonable  to 
assume  that  var{AlHi}  »  var{AlHo},  which  reinforces  the  utility  of  d;^  as  a  reliable 

distance  measure  between  the  two  (approximately  Gaussian)  distributions  of  A  under  either 
hypothesis. 

It  is  straightforward  to  show  that  (see  also  Appendix  H) 
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Nq  I)  U  1q(2T}j1  isignal-plus-noise  (5.10a) 
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where  Yh  is  the  hop-SNR  as  per  (5.5c).  Thus, 


(5.10c) 


-  (G-1)  ♦  e^”lo(2YH) 


(5.11a) 


(5.11b) 


Combining  the  above  and  the  independence  assumption  between  observables  we  conclude 


(5.12) 


Equation  (5.12)  will  be  used  for  evaluating  the  performance  of  decision  rule  (5.8)  under  the 
assumption  that  the  signal  does  indeed  come  on  one  of  the  candidate  frequencies  (no 
frequency  offset);  we  term  this  Case  I.  Clearly,  the  same  rule  can  be  employed,  neglecting 
the  fact  that  there  can  be  an  offset.  Use  of  (5.3)  under  the  assumption  of  a  worst-case 
frequency  offset  Af  =  Rh/2,  i.e.  the  signal  arriving  amidst  two  candidates  is  termed  Case  II. 


Two  more  cases  arise  from  the  possibility  of  adjusting  the  power  factor  S  in  (5.8)  by  a 
coefficient  a  <  1,  which  accounts  for  the  power  loss  on  the  signal  measurement  due  to 
worst-case  offset: 


V  n  H  1  H  4 
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Thus,  Casein  refers  to 
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(5.14) 


when  there  actually  is  an  offset  Af  =  Rh/2,  while  Case  TV  refers  to  the  adjusted  rule(5.13) 
when  there  actually  is  not  an  offset  (Af  =  0).  All  four  cases  are  summarized  in  Table  1, 
along  with  their  defining  features  (first  four  columns).  The  last  column  lists  the  associated 
distances,  in  analogy  to  (5. 12),  To  illustrate  the  steps  involved,  Case  II  is  analyzed  more 
thoroughly  in  Appendix  I;  the  other  two  involve  similar  calculadcms. 

Let  us  note  here  that  the  amplitude  adjustment  in  (5.13)  is  only  a  partial,  ad-hoc 
remedy  towards  the  frequency-misalignment  problem:  an  exact  approach  should  account 
for  the  fact  that,  whenever  Af  =  Rh/2,  two  adjacent  spectral  slots  simultaneously  receive 
significant  signal  power  (albeit  reduced  a )  and  should  proceed  to  construct  the  LR  from 
there.  However,  the  resulting  "optimal  rule"  is  overly  complex  and  not  measurably  better 
than  the  alternatives  examined  herein. 

The  optimal  rule  (5.8)  might  be  overly  demanding  to  implement.  An  alternative  and 
popular  rule  is  the  spectral-maximum  detector,  which  selects  the  maximum  from  the 
envelope  observables  R„  =  1  and  compares  it  to  a  threshold  [25, 26, 34, 39,  etc.] 


Table  1.  Definitions  and  Distance  Results  for  the  4  Cases  Examined 
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(5.15) 


Rule  (5.15)  affords  a  generalized-likelihood-ratio  tests  [12]  interpretation;  If  the 
narrowband  (hop)  SNR  Yh  =  STh/Nq  is  large,  the  spectral  observable  corresponding  to  the 

signal  location  can  be  expected  to  rise  well  above  die  noise  floor  and  be  the  maximum  with 
high  probability;  thus,  it  can  serve  both  as  a  frequency  estimator,  as  well  as  a  detection 
statistic  for  a  now  narrowband  (matched  Alter)  hypothesis  testing.  The  performance 
associated  with  (S.IS)  can  be  determined  exactly  (i.e.,  no  distance  need  be  evaluated)  by 
noting  that  rule  (5.15)  is  exactly  equivalent  to  the  rule 
(ifallRa,<T^  ,  decideHo) 

^  if  at  least  one  Rg,  ^  T,,.,  ,  decide  Hj  ^  (5.16) 

which  involves  the  comparison  of  each  R^,;  m  =  0,...,G-1  with  the  same  threshold 
Let  Pf,  denote  the  slot  false  alarm  probabiliOr 


(5.17a) 


and  denote  the  slot  detection  probability 


(5.17b) 


The  slot  probabilities  (P^,  Pf*)  are  then  related  to  the  DOC  (system)  probabilities  (Pp,  Ppy^) 


Pd  =  1-(1  -  Pd)(l-Pfa)°-‘  (no  offset) 


(5.19a) 


where  it  is  assumed  that  only  one  slot  contains  the  signal  under  ( no  frequency  offset). 
If,  on  the  other  hand,  a  worst-case  offset  is  assumed,  ( S.19a )  is  modiHed  to 

PD=l-(l-Pd;o)^(l-Pf.)°'2  (worst  offset)  (5.19b) 

where  the  notation  P^.^  is  adopted  to  indicate  the  signal-power  loss  by  a  (i.e.,  a  hop  SNR 
of  oyh). 

The  slot  probabilities  of  eq.  (5.17)  are  related  through  the  standard  Marcum's  Q 
function  [15, 16] 
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Equations  (5.18)-(5.21)  can  be  used  to  relate  PpA  and  Pd  directly:  from  (5.18)  it  follows 
that 

Pf.=  l-(1-PFA)''°  (5.22a) 


wherefrom 
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or  its  sampled  version 
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under  Hj,  and  noise  only  under  Hq,  create  a  set  of  spectral  estimates  IR^j  P;  m  =  1,...,G 

and  use  them  in  conjunction  with  any  one  of  the  rules  (5.8),  (5.15)  or  (5.26)  to  obtain  a 
binary  decision  (Hg  vs.  Hj).  These  spectral  estimates  can  be  created  by  any  judicious 

mechanism.  A  particularly  attractive  class  of  such  novel  spectral  estimators  is  based  on 
autoregressive  (AR)  models  [46].  Those  could  be  employed  to  produce  the  desired  spectral 
observables  in  an  efficient  manner,  according  to  the  formula^^ 
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(5.28) 


Here,  M  is  the  order  of  the  adopted  AR  model  and  a^;  m  =  1,...,M  are  the  estimated 
coefficients  derived  from  a  nonlinear  operation  on  the  raw  data  of  (5.9b)  and/or  a 
weighted  (filtered)  version  of  the  linear  prediction  residuals  [47].  Note,  that  for  the  AR 
spectrum,  it  is  the  area  under  a  peak  which  actually  reflects  the  power  of  the  corresponding 
harmonic  component 


lis  is  the  AR  spectrum  normalized  by  the  estimated  variance  of  the  linear  prediction 


errors. 


The  question  of  identifying  good  spectral 


which  the  FFT  and  AR  models  of  (S.9c)  and  (2.28)  are  two  examples,  is  indeed  quite 
open  and  therefore  exciting  in  its  own  merit  This  is  particularly  so,  because  most  of  the 
efforts  in  the  literature  address  the  improvement  of  spectral  estimators  (in  terms  of  bias, 
resolution,  smoothness,  line-splitting,  sidelobe  behavior,  etc.),  which  says  nothing  of  their 
detection  capability.  Here,  however,  we  leave  that  question  aside  and  proceed  to  examine 
decision  rules  which  do  not  explicitly  utilize  spectral  information;  instead,  these  rules 
operate  directly  on  data  from  the  correlation  domain,  thus  constituting  the  field  of 


Let  the  finite-record  sample  correlation  function  ^  be 


y(T)  =  r(t)r(t-x)dt;  0  ^  x  ^  Th  (5.29) 

where  r(t)  is  as  per  (5.1)  or,  more  specifically,  (5.6).  Since  r(t)  is  a  BP  process,  so  is  y(x ). 
Combining  (5.6),  (5.29)  and  rejecting  double-frequency  terms  leads  to  the  complex 
envelope  y(x )  of  y(x)  as 


V(x) 


^  rltlT*  (t-T)dt 


(5.30) 


where  means  "complex  conjugate."  Sampled  every  sec,  (5.30)  yields 

k  yikB-')  -  ^  /Jjg  rlt)7*  {t-kB'’)dt;  lc-0,...,G-1  (5.31a) 

Furthermore,  if  the  integral  in  (5.31)  is  approximated  by  a  finite  sum  in  terms  of  samples 
of  rtO,  we  arrive  at  the  biased  estimates  of  the  correlation  sequence 
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(5.31b) 


which  could  have  also  been  derived  by  direct  quadrature  demodulation  of  r(t)  and  complex 
sampling  of  the  resulting  7(t);  whether  this  approach  is  preferable  to  the  sequence  (5.29)- 

^^From  now  on,  the  term  "sample"  will  be  omitted;  "correlation"  will  always  pertain  to  time 
averaging,  not  ensemble  averaging. 


(5.31)  is  a  matter  of  implemcntational  convenience.  From  a  constant-falsc-alarm-rate 
(CFAR)  viewpoint,  it  is  typically  advantageous  to  use  the  normalized  correlation 
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(5.32) 


It  is  well  known  that,  within  a  constant,  the  sequence  the  periodogram 

{|RniP}G-i^are  Discrete  Fourier  Transform  (DFT)  pairs.  Thus,  given  (5.31),  a  DFT 

operation  would  allow  the  application  of  the  spectral  algorithms  discussed  in  Section  5.2. 
However,  it  is  possible  to  find  meaningful  detection  statistics  in  terms  of  the  y^'s  without 

an  explicit  evaluation  of  the  spectrum.  In  other  words,  we  can  state  the  following  (single¬ 
hop)  correlation-detection  problem;  Based  on  the  sequence  {yk}^'*  of  (5.31)  or 
equivalently,  on  the  set  [yo>tPk)®'*]  (5.32),  formulate  a  detection  statistic  without 

explicit  evaluation  of  the  spectrum 

Simple  examples  of  ad-hoc  detector  choices  proposed  in  the  past  arc  the  "Semicoherent 
Detection  Statistical  Test" 


-V 


threshold. 


(for  some  k) 


in  a  radar  environment  (see  [26]),  the  RF  correlation  statistic. 


(5.33) 
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(5.34) 


of  [32],  with  0  <  X  ^  1  and  {ak}^  an  appropriate  set  of  coefficients,  the  statistic  |  Pjp  of 

k  “  1 

[34],  etc.  Some  more  advanced  choices  will  be  discussed  later  in  this  section.  For  the  time 


98 


being  let  us  concentrate  on  the  analysis  of  (5.34),  whereupon  it  will  be  shown  that 
significant  gains  can  result  with  respect  to  the  radiometer.  Let  us  note  that  the  zero*  term 
lyoP,  which  is  in  fact  the  radiometric  output,  could  have  been  included  in  (3.54)  with  a 

marginal  improvement  in  the  detector’s  discriminatory  power;  on  the  other  hand,  this 
would  make  the  detector  much  more  sensitive  to  the  (possibly  fluctuating)  noise  variance 
and  is,  therefore,  omitted. 

Decision  rule  (5.34)  can  be  rephrased  in  terms  of  the  continuous-argument  correlation 
function  y(x)  of  (5.29)  as 
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(5.35a) 


where 
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by  virtue  of  the  fact  that  y2(x)|Lp  =  ly(x)l2/2.  In  other  words,  the  complex  baseband  rule 
(5.34)  can  be  reconfigured  in  terms  of  power  measurements  (around  x^)  on  the  real-time 
correlation  output  at  an  RF  or  IF  frequency  as  shown  in  Figure  19.  Again,  whether  the 
rule  is  actually  implemented  at  RF  or  baseband  in  discrete  continuous  time,etc.4S  a  matter 
of  practical  convenience. 

The  set  of  coefficients  {aij}°'i  can  be  chosen  according  to  any  particular  philosophy 

and  is  subject  to  optimization.  It  can  be  shown  that  performance  of  the  summation  (5.35) 
is  insensitive  to  the  exact  value  of  a^'s  and  a  number  of  reasonable  choices  would  work,  as 
long  as  the  upper  limit  of  the  summation  is  properly  truncated  at  a  level  S  =  XG,  which  is 
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less  than  G  (i.e.,  X  <  1).  This  is  because  the  last  samples  possess  such  an  increased 
variance  that  their  inclusion  in  the  summation  is  detrimental.  Analytically,  this  is  equivalent 
to  letting  a^  =  0  for  k  =  S  +  .  Thus,  a  proper  choice  of  X  is  rather  crucial,  while 

that  of  not.  For  reasons  that  will  soon  be  clear,  it  is  mathematically  convenient 

to  let 


Th  ‘  ■^k)'  •  k=l....,xG 

°  i  k-XG+1 . G 


(5.36) 


in  the  following  calculations.  We  note,  that  this  set  of  a^'s,  corresponding  to  the  inverse  of 
a  triangular  window  on  the  correlation  samples  y^  of  (5.31),  is  such  that  the  latter  samples 
become  unbiased  estimates  of  the  true  (ensemble)  correlation  function  of  the  underlying 
discrete-time  process^^  [46,48]. 

To  pursue  the  analysis  of  (5.35),  with  y(t)  as  per  (5,29),  let  us  assume  that,  under 
Hj,  the  input  to  the  correlator  r(t)  consists  of  a  sinusoid  V2S  cos(0}ct4G)>  whose 
frequency  coincides  with  the  center  frequency  of  the  observed  band,  plus  bandpass 
Gaussian  noise  with  the  standard  quadrature  expansion, 


r(t)  =  V2S  cos(cOgt  -f  [ni(t)  cos(t05t  +  G )  -  nqCt)  sin  (Ogt  +  G  ] 
0St<TH 


(5.37) 


In  (5.37),  ni(t)  and  nQ(t)  are  the  inphase  and  quadrature  lowpass  noise  processes, 
respectively,  which  are  i.i.d.,  zero-mean,  with  a  flat  PSD  of  Nq/2  W/Hz  and  two-sided 

24Such  unbiased  correlation  estimates  lead  to  power-spectral-density  estimates  which  are 
not  warranted  to  be  non-negative,  an  undesirable  feature  in  spectral  estimation  [46].  Yet, 
they  arc  quite  appropriate  for  our  detection  procedure  (5.35),  which  further  illuminates  the 
aforementioned  distinction  between  the  detection  and  estimation  goals. 


bandwidth  of  B  Hz  each  (total  power  per  component  =  NoB/2).  The  assumptioi  about  the 
signal's  frequency  is  just  mathematically  convenient;  it  will  be  shown  later  that  it  does  not 
alter  the  results  in  any  way.(see  the  baseband  complex  formulation  in  Appendix  K  as  well 
as  the  results  of  Section  5.4. 

Filtering  out  double  frequencies,  the  autocOTielator  output  y(x)  around  the  center 
frequency  is 

y(t)  =  J  r(t)r(t-t)dt  =  yss(^)  +  VsN  (“f)  +  yNs('^)  yNNW 

T 

(5.38) 

with 

yss(^)  =  S(TH-t)costDc^  (5.39a) 

ysN('t)=  ^  (Nf(x)cos©c''  +  NQ('')sino>c'')  (5.39b) 

yNs('')=  VS(Ni(x)cosa)e'c-NQ(x)sin<OeX)  (5.39c) 

yNwCc)  =  (Nii(x)  +  Nqq(x))  cos  ©c^  +  (Niq(x)-Nqi(x))  sin  (5.39d) 


where  the  correlation-domain  noise  processes  have  been  introduced,  i.  e., 
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nrj-j(t-T)dt 

A 


(5.40a) 
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(5.40b) 
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nj{t)  nj{t-T)dt 


(5.40c) 
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)  *  ^  ^  HpCt)  np(t-T)dt 


(5.40d) 


Niq(t) 


l(t)  np(t-T)dt 


(5.40e) 


/""c 


Nqj(t)  *  y  "  »'Q(t)  nj(t-T)dt 


(5.40f) 


Much  of  the  following  analysis  deals  with  the  statistical  characterization  of  the  above  noise 
processes.  For  instance,  it  is  clear  that  the  first  four  processes  in  (5.40a)  and  (5.40b)  are 
zero  mean  and  Gaussian,  while  the  remaining  are  not  Gaussian.  However,  they  will  be 
approximately  treated  as  such  in  the  analysis,  particularly  for  0  <x  «  ,  by  virtue  of  the 

following  central-limit-type  argument:  The  bandwidth  B  of  n](t)  and  nq  (t)  is  much  larger 
than  Th^'  ;  thus,  each  integral  in  (5.40c)  through  (5.40f)  can  be  approximated  by  a  large 

sum  of  noise-product  samples,  each  sample  taken  B'^  seconds  apart  Since  we  shall  be 
interested  in  values  of  x  that  are  multiples  of  it  can  be  shown  that  those  samples  are 

mutually  uncorrelated.  Although  this  does  not  imply  statistical  independence  (which  is  a 
prerequisite  for  a  rigorous  application  of  the  central-limit  theorem),  it  nonetheless  reinforces 
the  argument.  This  will  permit  us  to  calculate  higher  moments  using  Gaussian  identities 
whose  exact  evaluation  would  otherwise  be  unwieldy.  We  furthermore  note  that  simulation 
has  confirmed  the  validity  of  the  approximation  (see  Section  5.4) 

The  first  point  of  interest  is  the  mean  value  of  y(T)  in  (5.38).  We  shall  assume  that  the 

input  BP  filter  in  Figure  19  has  a  perfectly  rectangular  transfer  function  ("brick  wall"), 
which  implies  that  the  autocorrelation  function  Roq(x)  for  both  n](t)  and  nQ(t)  is  given  by 


{nj(t)nj(t-t)|  =  «?{nQ(t)nQ(t-T)|  = 


(5.41) 
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where  S,(x)  4  (sin  x)/x.  Using  the  fact  that 

'W  .  *{»,}  .  ,{1,51  =  .  ^{n„}.0  (3.42*) 

and 

rf{Nu(x)}  =  {Nqq(x)}  =  (Th-x)R„„(x).  (5.42b) 

we  conclude  from  (5.38)  through  ( 5.40 )  that 

<?{y(x)}  =  [S  +  (NqB)  S,(rcBx)](TH-x)  cos  ©<.1  (5.43) 

Certain  interesting  observations  can  be  made  on  (5.43).  First,  although  the  quadrature 
component  (sin  ©gt)  has  been  eliminated,  the  mean^{y(x)}still  includes  an  unknown 

oscillation  that  is  why  noncoherent  (power)  sampling  must  be  performed  at  the  signal¬ 
processing  unit  Second,  we  note  that  the  mean  value  at  zero  offset 
^{y(0)}  =  (S  +  No  B)Th  is  just  the  average  energy  measured  by  the  radiometer  at  the  end  of 
Th  second,  as  it  should  be.  Since  we  have  assumed  that  Yin«i»  unplying  that  S  «  NflB, 

it  follows  that  the  radiometer  output  is  dominated  by  the  noise  contribution.  Third,  if  the 
output  y(x)  is  power  sampled  at  multiples  of  B'^  i.e.,  x^  =  kB'i;k  an  integer,  the 
contribution  of  the  mean  noise  will  be  eliminated  due  to  the  nulls  of  the  S,(x)  function. 
This  fact  wiU  be  exploited  in  the  following. 

Let  us  return  to  (5.38)  for  a  moment.  Since  we  shall  be  interested  in  low-SNR  case, 
we  make  the  simplifying  assumption  that  the  "signal  times  noise"  terms  ysN('t)  yNs(x) 
are  negligible  compared  to  the  "noise  times  noise"  term  yNsC'^)-  Yin  below  -15  dB,  or 
so,  such  a  simplification  is  quite  justified.  Then,  y(x)  can  be  expressed  as 

y(x)  =  C(x)  cos  ©{jX  4  cos  ©pX  +  NQ‘kt)sin©cX  (5.44) 

^Recall  that  the  signal's  carrier  frequency  is  truly  unknown,  randomly  placed  in  the  B-Hz 
band.  It  was  modeled  here  to  coincide  with  the  l^wn  center  frequency  of  the  band  just  for 
analytical  convenience. 


•N  -N  . 


ivpapippi 


where 


C(t) 


(sCTh-0 
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Q^o) 


(5.45a) 


is  the  envelope  of  the  signal  component  and  the  equivalent  noises  are  defmed  by 


Nf  (t)  4  Njj(t)  +  NqqCt) 


and 


^  Njq(t)  -  Nqj(t) 


(5.45b) 


(5.45c) 


As  mentioned,  y(T)  will  be  power  sampled  at  \  =  kB-^;  at  which  points, 


=  '{n^’OkD}  •  0 


(5.46) 


According  to  our  previous  discussion,  Nj^^and  will  be  treated  as  Gaussian  noises. 

^  ^  eq 

Furthermore,  it  is  shown  in  Appendix  J  that  N|(tit)  (“'j)  ^  uncorrelatcd  for  every 

kj  =  1,...,G;  similarly,  the  pairs ,Nj^(;Ti<23)and  N®J^T|^23) 

are  uncoirelated  as  long  as  #  x,  .  Coupling  this  result  with  the  approximate  Gaussian 

12  G  ^ 

assumption,  it  follows  that  the  sequences  of  noise  samples  |n  ^  ^and 

contain  2G  Gaussian  rv's  which  are  zero  mean  and  approximately  independent.  Thus,  to 
complete  their  statistical  description,  the  second  moments  (or  variances,  for  this  case)  are 
required.  Those  are  obtained  in  Appendix  J,  wherein  it  is  shown  that 
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where  vv*^ 

F.(k)  *  f  (1-p')  -  tj^^pHdp' 

5  ^  (5.48a) 

FjfM  *  f^  (l-p')Sa|^G(Cl  -  cOp*  ♦  h)]  Sa[»6(G-tn>‘'H3‘^' 

0  :•■:•■?: 

(5.48b)  :  ::-:v: 

^  4 

and  the  normalized  parameter  is  defined  as 

•  *'  ^*1 
■*' .  \  •  J 

^  . c'::  4 

(5.48c) 

so  that  0  <  <  1.  Some  meditation  on  the  functional  form  of  F2(k)  reveals  that  it  has  a  i-  '.- V-'l 

negligible  contribution  compared  to  Fi(k),  so  it  will  henceforth  be  dropped.  From  (5.47),  '•‘  ‘■’'.  '‘'j 

*•  *1 

we  then  have  that 

'n'"'  *  ^{(''?''OkD)l  - 

•  •** 

®  f  i(k>  •  C"o'Th)^  0  -  tk)^  Fi(‘)  (5,49) 

u_4 

which  is  a  fairly  simple  expression,  but  also  fundamental  for  the  following.  In  deriving  the 

above,  it  was  assumed  that  (a)  timing  epoch  is  known,  e.g.,  under  Hj ,  the  signal  occupies  -I- 

L  < 

exactly  the  observation  period  and,  (b)  the  input  BP  filter  is  a  rectangular  one  with 

•*'  te*'' 

bandwidth  B.  The  first  assumption  can  be  removed  with  some  SNR  loss,  while  the 
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second  is  only  analytically  convenient;  the  results  can  easily  be  modified  to  account  for  any 
actual  filter  shape. 

Let  us  now  return  to  the  noncoherent  sampling  indicated  by  (3.3Sb).  Such  sampling 
can  be  implemented  via  a  square-law  device,  followed  by  a  wideband  lowpass  filter;  its 
bandwidth  should  be  several  times  B,  but  still  much  narrower  than  twice  the  carrier 
frequency.  Again,  neglecting  "signal  x  noise  terms’’^* ,  we  have  from  (5.44)  and  (5.45a) 
that 

“k  -  PT„3'Cl-'0*♦(HrOO)^(I<Q•''c-lO)^^c•l . 

(5.50) 

By  virtue  of  the  approximate  statistical  independence  of  the  second  and  third  noise  terms  in 
(22),  {Wk}^*  j  is  a  sequence  of  approximately  independent  noisy  samples,  upon  which 

the  decision  is  based.  As  mentioned,  this  valuable  property  of  independence  gradually 
diminishes  as  k  ^G-1  because  the  fundamental  assumption  in  assessing  it  (i.e.,  Gaussian 
equivalent  noise)  weakens  toward  the  end;  this  is  also  why  those  latest  samples  should  be 
ignored.  Nonetheless,  such  a  property  is  the  key  factor  in  the  anticipated  superiority  of  the 
correlator  versus  the  radiometer  -  namely,  the  fact  that  additional  data  can  be  exploited  by 
incorporating  all  these  new  samples  in  the  decision.  Their  approximate  independence 
prevents  these  samples  from  being  statistically  useless. 

In  order  to  proceed  with  the  evaluation  of  rule  (3.35),  we  need  the  mean  and  variance 
of  the  W^'s.  From  (5.47)  and  (5.50),  we  have  that  the  difference  of  the  means  under  the 
two  hypotheses,  A^{Wt},  is  given  by 

A<S{Wk}  ^  {WkIH,}  -^{WklHo}  =  (STH)2(l-;k)2  (5.51) 

^A  more  exact  analysis  is  provided  in  the  discrete-time  model  of  Section  5.4  without  any 
significant  deviation  in  the  results. 


Furthermore,  using  the  independence  between  and  (x^)  along  with  the 

Gaussian  assumption  yields  the  variance  as 

’•rKf  ■  {("fC’Of }  •  -  **i(«PCTO)^ 

with  Oj^j2  as  per  (5.49). 

As  seen  in  (5.51),  (5.52)  and  (5.49),  the  quantities  of  interest  for  include  a 
multiplicative  factor  (1- a  reflection  of  the  fact  that  the  integration  time  reduces 
proportionally  to  Xj.  as  x  Insertion  of  the  coefficients  a^  in  (5.36)  purported  to 

remove  this  factor  without  essentially  altering  the  results.  Thus,  if  we  define 
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it  follows  from  (5.51)  and  (5.52)  that 


(5.53) 


A<y{W’k}  =  S2 

(5.54a) 

and 

varjw^'}  •  4CNo8)^  Fj^(k) 

(5.54b) 

while  the  decision  rule  (5.35a )  reads: 

XG 

t  -  I  W^' 
k-l  ‘ 

If  we  now  define  the  output  SNR  for  the  correlator  in  the  familiar  way 


(5.55) 


SHR^orr  ^  ^2 

out  Y 


(e{Y  |Hi  }  -  e{Y  Ho  J 
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then,  using  the  uncorrelatedness  of  each  we  conclude  that 


(5.56) 


.  -  .  • 

^  - 


4  X6  2 
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(5.57a) 


C  .  . 


or,  using  the  fact  that  Yh  “  ^  *  Tin' 


SNRoSr  “  A(x.G) 


where  the  coefficient  A(X,,  G)  is  defined  as 


A(X,G)  S  ^ 


4G  I  F,^(k) 
k«l  ^ 


Equations  (5.5b)  and  (5.7b)  can  be  combined  into  a  ratio 


(5.57  b) 


(5.57c) 
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^  •  A(X.G).t„‘ 


(5.58) 


The  interpretation  of  (5.58)  is  that,  within  the  framework  of  performance  established 
herein,  the  correlator  will  outperform  the  radiometer  if  the  product  A(X,G)  •  Yh^  ‘S  above 

unity  and  vice  versa.  We  note  that  this  product  signifies  some  kind  of  utility  factor  and  is 
separated  into  two  components:  (a)  A(X,G),  depends  on  features  of  the  device  used  (G) 
and  the  signal-processing  algorithm  employed  (X),  and  (b)  Yh  =  ST^/Nq,  depends  on 

scenario  parameters  (signal  and  noise  strength,  hopping  duration). 

Exact  evaluation  of  A(X,G)  requires  numerical  integration  by  computer.  However,  we 
shall  develop  a  lower  bound  and  an  approximation,  which  are  quite  useful  for  a  wide  range 
of  X.  This  minimum  value  (or  lower  bound)  can  be  used  in  conjunction  with  (5.58)  to 

guarantee  a  minimum  gain  for  the  correlator  versus  the  radiometer.  The  argument  is  based 
on  the  observation  that  the  quantity  Fi(k)  in  (5.48a)  is  a  monotonically  increasing  function 
of  the  argument  thus, 

Fj(k)  .  J  (l-e')S»2[;G(l-x)p^dp'  ;  k-l . xG 

(5.59) 

since  =  X  from  (5.48c).  The  inequality  in  (5.59)  can  be  strengthened  by  the  following 


steps: 


(1-p' )  Sa^pG(l-x)p^dp’ 

0  I-  J 

<  f  Sa2[^G(l-x)p^dp'  .  G'^(l-X)'^  Sa^Cirxldx 


<  G  (1-x)  ^  ^  Sa^[irx]dx 


G'^(l-X)‘'(l/2) 


(5.60) 


which,  when  combined  with  (5.57c),  yields 


A(X,G)  >  A^„(X,G)  =  A^  (X)  =  X(l-X)2 


(5.61) 


We  note  that  the  lower  bound  Amjn  depends  only  or. while  the  exact  gain  A(X,G)  also 
depends  on  G.  This  is  just  a  fortunate  coincidence,  and  indicates  that  A„i„  can  be  used  for 
every  G,  as  long  as  G  »  1.  Thus,  we  arrive  at  the  following  inequality 
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(5.62) 


which,  as  mentioned,  should  be  tight  for  small  X . 

In  addition  to  this  general  bound,  we  can  create  a  good  approximate  formula  for  the 
above  ratio  when  G(1  -  X)  »  1,  as  follows;  since  the  steps  used  to  derive  (5.50)  are,  in 
fact,  tight  successive  approximations,  we  conclude  that  Fi(k)  ->  G’*(l-k.G'i)’i(l/2),  which 
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m»G(l-X)lia‘^  J 


leads  to 


The  last  summation  can  be  well  approximated  by  the  integral 
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m=6(l-x) 
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G(l-X)  * 


provided  that  G(1  -  X)  »  1.  Substitution  into  (5.57c)  leads  to 
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(5.63) 


which  differs  from  the  bound  (5.62)  in  the  exponent  of  (1  -  X).  For  low  X,  (5.62)  and 
(5.63)  agree  closely.  In  any  case,  those  two  equations  indicate  that  the  relative  merit  of  the 
correlator  increases  proportionally  to  the  square  of  the  hop  SNR  and  could,  therefore,  reach 
significant  levels,  depending  on  the  application  (coded  or  uncoded  systems,  slow  or  fast 
hopping,  etc.) 

The  RF  algorithm  (5.34)  -  (5.35)  analyzed  above  is  but  one  attempt  to  exploit  the 
information  imbedded  in  the  correlation  coefficient  {y^}®  •  In  fact,  it  can  be  viewed  as  a 

nonparametric  algorithm,  in  that  no  particular  model  for  the  underlying  process  has  been 

assumed  in  formulating  it;  the  algorithm  simply  measures  the  total  "power"  of  the 

correlation  lags.  As  an  attempt  at  a  more  structured  approach,  consider  modeling  the  time 

domain  process  an  autoregressive  (AR)  process  of  order  M,  i.e. 

k  “0 


^  M.in 


(5.64) 


In  (5.64),  {e^}  represents  the  zero  mean  i.i.d.  ("white")  driving  noise  process  with 
estimated  variance  a\  and  **  estimated  AR  parameter  vector; 


r 


# 


the  normalized  AR  spectrum  (by  o| )  resulting  from  (5.64)  is  expressed  in  (5.28).  The 

fact  that  (5.64)  can  spectrally  model  arbitrarily  closely  any  second-order  ergodic  process 
of  finite  variance  as  G-»«»  (and  possibly  as  M— >««»)  serves  as  a  motivation  for  this 
approach.  Furthermore,  since  the  estimation  of  parameters  for  an  AR  model  results  in 
linear  equations,  it  has  a  computational  advantage  over  moving-average  (MA)  or  ARMA 
formulations. 

From  the  detection  viewpoint,  the  postulation  of  (5.64)  poses  the  interesting  problem 
of  determining  decision  statistics  associated  with  it  One  such  approach  was  outlined  in 
Section  5.2;  furthermore,  Kay  [34]  showed  that,  for  a  one-pole  (M  =  1)  modeling,  the 
spectral-maximum  rule  (5.15)  is  equivalent  to  the  threshold  comparison  of  the  estimated 
(single)  pole  magnitude 


threshold 


(5.65) 


which  is  a  normalized  variant  to  rule  (5.33).  This  fact  of  exact  equivalence,  which 
constitutes  an  interesting  conceptual  bridge  between  the  spectral  and  correlation  detection 
domains,  is  hard  to  generalize  to  higher  M.  Unfortunately,  the  parametric  performance  of 
the  one-pole  spectral  model  is  inferior  to  the  nonparametric  FF/spectral-maximum 
combination  for  the  nanowband  signal  case  [34].  Similarly  discouraging  simulation  results 
were  obtained  in  [38,  39]  for  higher  order  M  ^  I.  In  the  case  of  [38],  a  straightforward 
extension  of  the  statistic  (5.65)  was  attempted  in  the  form 

M  ^ 

I  <  tl'^esholcl 

”o  (5.66) 

wherein  it  was  found  that  for  a  small  time-bandwidth  product  performance  somewhat 
worsens  for  increasing  M;  on  the  other  hand,  when  the  signal  spectrum  broadened,  due  to 
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random  phase  modulation,  performance  of  (5.66)  was  superior  to  the  periodogram.  The 
conclusion  from  the  above  should  not  necessarily  be  that  AR  modeling  is  inappropriate; 
rather,  that  a  search  for  more  meaningful  detection  statistics  is  in  order.  For  more  detailed 
information  see  Appendix  K. 

Let  us,  also,  note  that  AR  modeling  of  either  the  colored  alone  [35]  or  ARMA 
modeling  of  a  Gaussian  signal  plus  noise  [36]  has  been  recently  employed  as  a  tool  for 
signal  detection  under  modeling  uncertainties.  In  the  latter  case,  the  classification  of  the 
algorithm  in  the  context  of  this  report  is  not  clear,  since  the  periodogram  is  used  as  an 
intermediate  step  in  a  spectral  matching  procedure  for  the  estimation  of  the  ARMA 
parameters.  In  both  cases,  however,  a  LR  is  formed,  either  in  its  GLRT  form  or  as  a 
whitening  matched  filter  (asymptotically  equal  to  an  SNR  measure).  When  the  signal  is 
Gaussian  and  broadband  (but  not  exactly  white),  Kay  [37]  suggests  the  use  of  a  GRLT, 
based  again  on  estimating  aj^. 

We  conclude  this  part  by  suggesting  an  alternative  modeling  approach  which  has  been 
shown  to  be  promising  in  the  very  broadband  case.  (See  Appendix  K).  The  approach 
consists  of  modeling  the  correlation  sequence  y^  as  an  AR  process  itself 

Vk  *  I  k-M+t . G-i 

(5.67) 

and  creating  decision  statistics  associated  with  this  model.  For  instance,  an  M  =  1  model 
could  give  rise  to  a  decision  rule  similar  to  (5.65) 


R,(0) 


> 

< 


threshold 


(5.68) 
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where  now  {Ry(k)}j|G  pertains  to  the  correlation  of  the  sequence  {y^}.  In  this  coirelation- 
of-the-correlation  domain  (C^D),  the  aforementioned  rule  (5.34)  can  be  interpreted  as 


Y  = 


Ry(0)  ^  threshold 


(5.69) 


with  a^  =  1.  We  note  that  (5.68),  or  any  other  similar  rule  in  the  C^D  domain  involves 
nonlinear  operations  on  the  data  7k.  in  contrast  with  the  linear  operation  of  the  FFT, 
required  to  produce  the  spectrum  =  IR^P.  The  tradeoff  between  the  merit  resulting  from 
such  nonlinear  operations  versus  the  computational  burden  implicit  in  them  is  examined 
further  in  Appendix  K. 


I 


r 


I 
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In  order  to  substantiate  the  value  of  the  different  approaches  discussed  so  far,  we 
compare  the  performance  of  certain  algorithms  in  detecting  a  sinusoid  in  a  wideband 
scenario  Gafgc  G).  The  comparison  here  is  by  no  means  exhaustive;  rather,  its  purpose  is 
to  give  a  preliminary  flavor  of  the  differences  in  performance  between  approaches.  A  more 
comprehensive  comparison  is  offered  in  Appendix  K,  along  with  a  discussion  of  the 
associated  computational  load  of  each  algorithm,  both  for  low  and  high  TB  products.  We 
note  that  the  large-TB  case  fits  well  the  spread-spectrum  frequency  hopping  detection 
model  described  without,  of  course,  being  confined  to  that 

We  first  look  at  the  performance  of  rule  (5.34),  an  order  to  assess  the  accuracy  of 
our  analysis. 

The  algorithm  has  been  simulated  on  a  digital  computer.  In  the  simulation,  the 
autocorrelation  function  y(T)  of  (5.38)  has  been  carried  out  by  its  discrete-time  equivalent  in 
terms  of  noise  samples  of  njft)  and  ngft).  Since  the  input  process  is  bandlimited  to  B,  the 
analog  operation  of  (5.38)  can  be  simulated  by  its  equivalent  discrete-time  system  with 
sampling  intervals  equaling  B'^. 

In  terms  of  discrete-time  samples,  the  radiometer  ouq)ut  in  (5.26)  is  modeled  as 


(5.69) 

where  G  <=  BT^  is  an  integer  and  ni(k/B),  ngfk/B)  are  quadrature  noise  samples.  Equation 

(5.69)  is  easily  seen  to  be  proportional  to  the  time  average  of  the  squared  envelope  of  r(t) 
over  the  hop  time  T^ 


Similarly,  the  discrete-time  simulation  equation  for  the  autocorrelation  detector  can 


be  expressed  as  follows: 


^CORR  * 


1 

Q-i+1 


[(G-IC+1)S  ♦  e(lc)  ♦J  n(k) 


(5.70) 


where  j  =  '/T,  is  complex  magnitude  squared,  and  ^(k),  q(k)  are  contributions  from  the 
SxN  and  NxN  terms  given  in  (5.39), 

(5.71a) 


and 


(5.71b) 


which  arc  zero-mean  noise  processes  for  the  assumed  filter  characteristics  of  Figure  19. 
Note  that  this  discrete-time  model  incorporates  all  possible  terms  without  any 
simplification. 

For  the  discrete-time  models  (5.69)  -  (5.71),  we  have  analytically  evaluated  the 
means  and  variances  of  the  decision  variable  y  under  either  hypothesis  i  =  0,1.  This 

not  only  allows  for  a  step-by-step  comparison  with  the  corresponding  simulation 
outcomes,  but  also  permits  a  more  precise  performance  description  in  terms  of  (Pj},  Pfa)- 

The  theoretical  results  appear  on  Table  2.  The  simulation  results,  along  with  their 
theoretical  (numerical)  counterparts  from  Table  2,  appear  in  Tables  3  and  4  for  the 
radiometer  and  the  autocoirelator,  respectively.  The  last  entry  of  Tables  3  and  4, 
is  calculated  on  the  basis  of  (5.56),  which  evidently  neglects  var{Yl  Hj}.  The  close 


Table  2.  Moments  of  Detector-Output  Variables  Evaluated  from 
the  Discrete-Time  Model  Equations  (5.69)  and  (5.70) 


Moments  of  I 
Detector-Output 
Variables 

Radiometer 

Autocorrelator* 

<?CYlHo) 

2NqB 

CW'  E^*  1"  (tIx) 

<?CY|H£) 

2CS  +  NqBD 

XGS^  *  j2SNQB+(t(oB]?^  T 

varCYlH^ 

1  (: 

1  •  \  (t^)CNo»/ 

varCY|Hi3 

C2NoBJ^  /  S  \ 

G  *  2  sj?) 

+  CNoB)  (’**><?) 

Table  3,  Comparison  of  Radiometer  Simulation  Results  to  Theory  (200  Trials) 


BT„  -  1000 
H 


1  Input  SNR 

-20  dB  1 

-15  dB  I 

-10  dB  1 

Theory 

Simulation 

Theory 

Simulation 

Theory 

Simulation 

«(y1ho) 

4.00 

4.003 

4.00 

4.003 

4.00 

4.003 

A. 04 

4.044 

4.126 

4.132 

4.40 

4.407 

var(Y|Ho] 

16.00 

16.14 

16.00 

16.14 

16.00 

16.14 

»ar(ylHi) 

16.3 

16.20 

17.0 

16.69 

19.2 

18.57 

Radiometer 

SNRo  (dB) 

-10.00 

-9,73 

0.00 

0.14 

10.0 

10.06 

Table  4.  Comparison  of  Autocorrelation  Detector  Simulation  Results 
to  Approximate  Analysis 

BTjj  -  1000;A.  -  0.1 


^in  (dB) 

-20  dB 

-15  dB 

-10  dB 

Theory 

or  Simulation 

Theory 

Simulation 

Theory 

Simulation 

Theory 

Simulation 

0.417 

0.419 

0.417 

0.419 

0.417 

0.42 

0.47 

0.472 

0.84 

0.863 

4.46 

4.57 

1.76x10'^ 

2.61xl0'‘^ 

1.76x10-3 

2.6x10"^ 

1.76x1*0’^ 

2.61x10-3 

3.37x10"^ 

4.97x10"^ 

4.58x10"^ 

5.34x10"^ 

1.322 

1.41 

Autocorrelation 

Detector 

SNRq  (dB) 

1.16 

0.17 

20.06 

18.76 

39.7 

38.2 

numerical  agreement  between  theory  and  simulation,  over  a  wide  range  of  values,  can 
be  observed  from  the  tables. 

It  is  further  interesting  to  note  that  the  difference  in  variance  of  y  under  the  two 

hypotheses  is  virtually  negligible  for  the  radiometer,  while  it  becomes  quite  pronounced  for 
the  correlator  as  Y|q  increases.  The  implication  is  that  SNRo^  is  well  suited  as  a 

performance  measure  for  the  former,  but  not  so  much  for  the  latter.  A  comparison  between 
theory  and  simulation,  in  terms  of  the  SNR^m  dLTerence 

as  derived  in  (5.63),  is  provided  in  Table  5.  Certain  parameters  were  chosen  so  that  the 
theoretical  assumptions  involved  in  the  analysis  would  be  well  justified:  (1)  large  time- 
bandwidth  product  (G  *=  KXK)  or  30  dB),  (2)  small  input  SNR  (Yin  ^  -10  dB),  (3)  small 
X(X  =  0.1)  and,  (3)  large  number  of  trials  (200),  so  as  to  assume  statistical  confidence. 
Since  those  numbers  were  in  accordance  with  the  assumptions  made,  any  deviations  of  the 
simulation  from  the  theory  would  indicate  an  error  in  the  modeling  process  (Gaussianness, 
etc.).  Fortunately,  that  was  not  the  case,  as  evidenced  by  the  closeness  of  the  theoretical 
and  simulation  results  in  Table  5.  It  is,  also,  clear  that  significant  gains  can  be  attained 
over  the  radiometer,  particularly  as  Yin  increases. 

Next,  we  look  at  the  sensitivity  of  the  decision  rule  (5.55)  with  respect  to  the  choice 
of  k.  This  is  illustrated  in  Figure  20,  where  the  detection  probability  is  plotted  versus 

X  (0  <  X  ^  1)  for  two  values  of  Yin  while  G  =  100  is  constant;  fw  each  set  of 
three  curves  associated  with  a  value  of  Yin.  the  corresponding  values  of  are  10**,  10^ 

and  10'^,  from  top  to  bottom,  respectively.  In  all  cases  examined,  the  trend  is  clear; 
performance  is  insensitive  to  change  of  X  over  a  very  wide  range,  as  long  as  X  is  neither 


Table  5.  Theoretical  and  Simulation  Values  for  (ASNR) 


For  Various  Values  of 


G  =  10^,  X  -  0.1 

y  IdB) 

in 

-10 

-15 

ASNR  (dB) 
theoretical 
(lower  bound) 

29 

19 

(dB) 

Simulation 


28 


18.5 


too  small  nor  too  large.  So,  any  X  in  the  range  (0.1,  0.8)  would  be  a  suitable  choice. 
Although  the  curves  illustrated  in  Figure  5  arc  for  a  particular  choice  of  a^'s  from  (5.36); 
an  identical  kind  of  behavior  was,  also,  found  for  other  choices.  This  is  a  rather 
convenient  conclusion  since  one  does  not  have  to  be  concerned  about  optimizing  X 
whenever  the  constants  are  changed.  Furthermore,  it  was  found  as  a  more  general 
conclusion  that  the  overall  performance  (Pd,  Pfa)  insensitive  to  the  particular  choice  of 
a^’s  .  Thus,  one  can  set  a^  =  1,  k  =  1,2,...,XG,  which  would  simplify  the  decision  rule  to 
a  mere  accumulation  of  the  power  samples  and  a  threshold  comparison. 

A  suitable  set  of  design  curves  for  both  the  radiometer  and  the  correlator  are 
provided  in  Figures  21  and  22.  In  these  figures,  the  input  SNR's  (Yin)  required  to  achieve 

certain  detector  operating  characteristics  are  shown  for  both  algorithms.  The  probabilities 
of  false  alarm  (Ppy^)  and  detection  (Pd)  are  computed  directly  from  the  first  two  moments  of 

the  detector-output  variables  under  either  hypothesis  (Tables  3  and  4).  Gaussian 
approximation  is  used,  which  is  justifiable  since  in  either  algorithm  the  detector  output 
consists  of  the  sum  of  a  large  number  of  random  variables.  With  this  assumption,  Pd  is 
given  in  terms  of  PpA  and  the  moments  of  yas 


L  /vartYlH,)  J 


P0=Q 


(5.73) 


From  Figures  21  and  22,  the  relative  gains  of  the  correlator  over  the  radiometer  in  terms  of 
Yin  can  be  assessed.  These  results  are  plotted  in  Figures  23  ,  24  and,  also,  summarized 

in  Table  6.  Note,  that  the  performance  improvement  of  the  correlator  over  the  radiometer  is 
significant  in  the  (Pd,  Pfa)  range  of  interest.  For  example,  for  BT^  =  1000,  Pd  =  0.9,  Pp;^ 
*  10-3,  the  correlator  will  require  6.8  dB  less  Yin  *^an  the  radiometer.  This  performance 
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Table  b.  Relative  Improvement  of  the  Autocorrelation  Technique  Over  the 
Radiometer  Approach  in  terms  of  Required  Input  SNR 
For  Desired  Pp  and  Pp^  Performance  level.  " 


Improvement  of  Autocorrelation  Technique 
Over  Radimcter  (dB) 


0.9  P_  .  0.8  Pp  -  0.7  Pp  -  0.6  Pp  -  0.5 


11 

ft 

vs^:  ■: 

'•  -‘‘Vv,-'' 
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improvement  will  be  8.3  dB  for  the  operating  characteristic  Pd  =  0.9,  PpA  =  lO"^  (which 
requires  a  higher  than  Pd  =  0.9,  PpA  =  • 

Finally,  the  RF  correlation  and  radiometric  results  are  compared  with  the  various 
spectral  algorithms  of  section  (5.2)  in  Figures  25  a,b,c  (BT^  =  100)  and  26  a,b,c 
(BTh=1000)  for  various  Pfa  levels.  In  particular,  the  curves  on  these  figure  designate 
the  following  cases: 

Curve  1 :  Optimal  LR  algorithm  with  Af  =  0  (Case  I  of  Table  1 ) 

Curve  2:  Power  Adjusted  Rule  with  Af = 0  (Case  IV  of  Table  1) 

Curve  3:  Power  Adjusted  Rule  with  Af = Rh/2  (Case  HI  of  Table  1 ) 

Curve  4:  Optimal  LR  Algorithm  with  Af  =  Rh/2  (Case  n  of  Table  1) 

Curve  5:  Spectral-Maximum  Detector  with  Af  =  0  (eqn.  5.16) 

Curve  6:  Spectral-Maximum  Detector  with  Af  =  Rh/2  (eqn.  5.16) 

Curve  7:  The  Autocorrclator  (which  is  basically  not  affected  by  Af) 

Curve  8;  The  Radiometer(which  is  not  affected  by  Af). 

The  performance  of  the  LR  algorithm,  as  illustrated  in  these  figures,  is  optimal  if 
Af  =  0  only.  With  the  worst-case  frequency  offset  (Af  =  Rh/2)  a  degradation  close  to  5  dB 

will  be  suffered.  With  the  power-adjusted  rule  this  degradation  can  be  somewhat 
minimized:  resulting  in  degradations  (from  optimal)  of  approximately  0.7  dB  and  3.5  dB, 
respectively,  for  Af  =  0  and  Rh/2. 

The  performance  of  the  spectral-maximum  detector,  also,  suffers  a  significant 
(approximately  3  dB)  degradation  when  Af  =  Rh/2. 

The  LR  algorithm  exhibits  a  pronounced  threshold  behavior:  when  the  hop  SNR 
is  above  a  threshold  value  (which  depends  on  the  BTh  and  PpA  setting),  the  probability  of 
signal  detection  Pd  approaches  unity  rapidly.  On  the  other  hand,  if  Yh  below  this 
threshold  value,  Pd  diminishes  to  zero. 


The  other  algorithms'  (i.e.,  the  spectral-maximum,  the  autocorrelator,  and  the 
radiometer)  performances  are  more  robust,  they  do  not  exhibit  pronounced  threshold 
behavior  like  the  LR  detector.  In  fact,  for  large  BTj|  (BTh  =  1000),  moderate  Tij  (Yh  <  10  dB) 
and  acceptable  Pfa(^  lO"^),  the  performance  of  the  spectral-maximum  detector  and  the 
autoconelator  both  exceed  that  of  the  LR  detector,  when  worst-case  Af  is  considered. 

Since  the  cases  of  interest  in  a  LPI  detection  environment  actually  correspond  to 
large  BTjj,  moderate  Yh,  and  PpA  <  lO"^.  and  since  Af  cannot  be  expected  to  be  zero,  it 

appears  that  spectral-maximum  detector  and  the  autocorrelator  are  "better"  LPI  detectors 
than  the  LR  detector,  due  to  their  robustness,  though  the  LR  detector  is  "optimal"  in  the 
theoretical  sense,  assuming  perfect  frequency  and  power  aligmnents. 

The  most  interesting  observation  to  be  made  on  these  figures  is  that  in  the  region  of 
interest,  described  above,  the  autocorrelator  appears  to  be  the  best  choice,  since  it  is  not 
affected  by  Af  and,  thus,  performs  better  than  the  spectral-maximum  detector  when  the 
worst-case  Af  =  R^/S  is  considered. 

For  a  numerical  discussion  and  simulation  comparisons  involving  the 
aforementioned  novel  AR  schemes,  the  interested  reader  is  referred  to  Appendix  K. 
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Thus  far,  the  interference  model  has  consistently  been  additive  broadband  Gaussian 
noise.  As  mentioned  in  the  Introduction,  this  is  hardly  expected  to  be  the  case  in  a  dense 
environment  where  the  sought  signal  coexists  with  a  variety  of  narrowband  and  wideband 
waveforms,  some  of  which  vary  quite  unpredictably  with  time.  Interference  excision  can 
work  well  whenever  sufficient  discrimination  exists  between  the  wanted  and  unwanted 
signals.  Such  discrimination  results  from  distinctive  features,  be  those  long-term 
(multihop)  spectral  information,  directional  information,  different  rates,  etc.  In  this 
section,  we  shall  assume  that  everything  possible  (from  an  excision  viewpoint)  has  been 
done,  and  what  is  left  is  some  random,  unpredictable  narrowband  interference.  Our 
purpose  is  to  show  that,  even  in  such  an  environment,  intelligent  processing  can  provide 
good  detection  capability  while  brute-force  radiometric  detection  could  fail  miserably.  For 
illustrative  purposes  we  shall  assume  that  the  thermal  noise  level  is  negligible  compared  to 
the  narrowband  interference  power.  Furthermore,  the  sought  signal  will  have  a  DS 
modulation  on  it,  either  by  itself  or  in  a  hybrid  DS/FH  form.  In  the  latter  case,  the  spectral 
segments  of  Figure  17  must  have  enough  bandwidth  to  accommodate  the  DS  modulation. 
Since  we  are  interested  here  in  single-hop,  per-band  waveform  processing,  it  follows  that 
the  FH  part  of  the  modulation  does  not  explicitly  enter  in  the  mathematical  model,  except 
through  the  center  frequency  of  the  observed  band.  In  that  sense,  both  pure  DS  and  per- 
band  DS/FH  (single-hop)  in  narrowband  interference  result  in  the  same  mathematical 
model.  Let  us  remark  that  the  present  model  (wideband  signal  in  narrowband  noise)  is,  in 
a  sense,  the  inverse  to  the  problem  tackled  in  Section  5.2  -  5.4,  therein,  having  a 
narrowband  signal  (FH)  in  wideband  Gaussian  noise.  In  view  of  the  encouraging  results, 
in  these  previous  sections,  it  is  to  be  expected  that  a  good  algorithm  could  also  exploit  the 
difference  between  signal  and  "noise"  in  the  present  case,  too.  In  fact,  it  will  then  be 
shown  that  for  observation  data  with  a  large  time-bandwidth-product  G,  a  very  simple 
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algorithm  operating  on  the  output  of  a  real-time  autocorrelation  device  can  achieve  almost 
perfect  performance-in  stark  contrast  with  the  poor  performance  of  the  radiometer. 

Hybrid  modulation  schemes  have  become  increasingly  popular  spread-spectrum 
communication  choices,  due  to  the  enhanced  andjam  margin  which  they  offer.  Adding  DS 
modulation  to  FH,  also,  improves  the  anti-intercept  capability  because  the  "noiselike" 
appearance  of  DS  makes  detection  more  difficult.  On  the  other  hand,  multiple  tones 
constitute  a  common  model  of  nonwhite  interference  and  can  emerge  in  a  number  of 
scenarios,  i.e.,  it  can  be  intentional  (jamnung  of  the  band)  or  unintentional  (multiple  users 
in  a  broadcast  environment,  adjacent  radar  sources,  etc.).  It  is,  also,  conceivable  that  the 
tone  interference  has  been  deliberately  inserted  by  the  communicator  in  a  pseudorandom 
manner,  so  as  to  impede  the  interceptor's  task,  while  it  can  be  pseudorandomly  avoided  by 
the  intended  receiver. 

Let  the  FH/DS  hybrid  signal  to  be  detected  be  represented  by  s(t)=-\/2S  c(t)  cos  (Dot, 
where  Wq  is  an  unknown  frequency  within  the  observed  spectral  band,  c(t)  is  the  DS  code 
of  rate  and  S  is  the  signal  power.  The  unknown  interference  consists  of  M  tones 

(M  is  a  random  variable  in  each  observation  interval,  which  is  equal  to  the  hop  time 
Th  =  R'h  ).  with  ©jj  and  <1)^  denoting,  respectively,  the  power,  radian  frequency  and 
phase  of  the  Jith  tone.  The  total  received  signal  in  (0,  T^),  under  hypothesis  Hj  (signal 
present),  is  given  by 


M 

r(t)  =  fTE  c(t)  cos  u)f,t  +  ^  JZIT  cos('u  t  + 

k=l  (5.74) 

Note,  the  absence  of  thermal  noise  in  (5.74),  as  previously  discussed.  The  code  c(t)  can  be 
modeled  as  either  a  random  sequence  of  independent,  identically  distributed  ll's  with 
Pr[c(t)  =  1]  =  0.5,  or  a  PN  code  with  a  full  period  equal  to  T^.  As  the  ratio 
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N  4  Th/Tc  =  Rc/Rh  increases  (N  denotes  the  number  of  code  chips  per  hop),  the 
performance  difference  becomes  insignificant,  a  fact  also  verified  by  simulation. 

Let  B  denote  the  input  observation  bandwidth.  The  presence  of  the  DS  code  implies 
that  B  should  be  at  least  equal  to  Rg  or  higher,  but  definitely  much  larger  than  Rh. 
Equivalently,  the  time-bandwidth  product  G  ^  BTh  »  1.  Furthermore,  for  simplicity, 
we  shall  assume  that  all  tones  have  equal  power  =  I/M,  k  =  1,...,M,  where  I  is  the  total 

interference  power,  and  they  are  equi-spaced  within  the  bandwidth  B.  In  other  words,  the 
frequency  separation  |f^-  f^+jl  between  adjacent  tones  equals  B/(M  -i-  1),  which  is  much 
greater  than  Th,  i.e.,  If^  -  fk+il  »  Th.  None  of  the  above  assumptions  are  critical  in  the 
forthcoming  conclusions;  they  simply  ease  the  analytical  burden. 

The  real-time  autocorrelation  operation  produces  the  output 

y(T)  =  y  ^  r(t)  r(t  -  Tjdt  ;  0  <  t  < 

(5.75) 

Substituting  (5.74)  into  (5.75)  and  rearranging  redefines  y(  T )  as 

(5.76) 


where  yc(x)  is  the  code  partial-correlation  function  (a  random  variable) 


c(t)  c(t  -  T )dt 


(5.77) 


and  nj*(t),n“.k  =  1,...,M  are  approximately  Gaussian  (via  a  central  limit  theorem-type 
argument),  bandpass  noise  processes  (signal  x  interference  terms)  defined  by 


n,^{T)  *  a|^(T)  cosuijT  - 


(5.78a) 


and 


where 


and 


nJ(T)  = 


6‘(T)j  cosx^T  .  I  2^  6j(T)j  s1n»j 


(5.78b) 


aJ(T)  =  y*  ^c(t)  cos  jj;;a)Q  -  u.^)t  - 


(5.79a) 


“  y  ^  c(t)  sinjj;a.Q-u^3t  - 


(5.79b) 


^J(t)  =  y*  ^c(t  -  t)  cos  (wq  -  ^ 


(5.79c) 


gQ 


=  y  ”c(t  -  t)  sin  (uq -u)|^3t  -  ♦Jldt 


(5.79d) 


Before  discussing  the  statistical  characterization  of  the  above  noise  processes  ,  let  us 
examine  the  noiseless  (mean)  part  of  y(T)  in  (5.76),  as  shown  in  Figure  27.  Of  particular 
interest  here  are  the  envelopes  of  the  useful  signal  (small  shaded  triangle)  and  the 


Figure  27.  Mean  of  Output  y(T):  The  Signal  and  Interference  Additive  Components 


interference  (large  triangle),  respectively,  since  the  actual  components,  e.g.,  the  first  two 
terms  in  (5.76),  are  modulated  by  the  unknown  frequencies.  The  structural  difference 
between  the  two  correlations  is  evident:  the  DS  code  superimposed  on  each  hop  creates  a 
narrow  mean  autocorrelation  function,  since  the  expected  value  of  the  function  yjix)  in 
(5.77)  is  zero  for  It)  ^  T^.  Contrary  to  that,  the  interfering  tones  correlate  for  the  whole 
interval  [0,7^].  Clearly,  then,  a  power  (noncoherent)  sample  Xj  =  would  measure 
interference  only^"^  this  sample  could  be  subtracted  from  the  corresponding  one  at  Tq  =  0, 
so  that  under  H^,  contains  the  full  signal  power  plus  interference.  This  subtraction  would 
approximately  cancel  the  interference  contribution  at  Tq  =  0  so  that,  under  Hj,  only  the 
signal  would  emerge  while,  under  Hq,  the  statistic  would  be  almost  zero.  Thus,  the 
adopted  decision  rule  is  (see  Figure  28) 


A  =  y2(o)  -  y^(t  J  ^ 

I  n  ^  u 


LP 


(5.80) 


where  dg  is  a  fixed  threshold.  In  the  absence  of  thermal  noise  dg  can  be  set  at  a  very  small 

(positive)  level  in  order  to  maximize  the  detection  probability 

We  shall  now  investigate  the  performance  of  nile  (5.80).  Let  us  first  address  the 
ft  B 

statistics  of  n^(x)  and  ng(x),  defined  in  (5.78)  and  (5.79).  Let  <?{•}  and  var{»}  indicate  the 
mean  and  variance,  respectively.  It  is  then  easily  shown  that 


j-I.O 


(5.81) 


and  that 


2'^Any  Xi  =  Tg  would,  also,  do  especially  in  the  face  of  some  uncertainty  about  T^. 


^ 


var^Br 


(5.82) 


In  the  following,  we  shall  concentrate  on  offsets  t  =  rnTg;  m  =  0,1,2,...  In  fact,  only 
T  =  0  and  T  =  are  of  interest  (see  (5.80)).  Then,  lengthy  manipulations  can  establish  the 
following  facts: 


(a)  The  variance  of  (5.82)  is  given  by 


Kc"TcD)‘ 


(N-m)T 


.  2 


(5.83) 


where  S,(x)  =  sin  x/x  is  the  sampling  function. 

(b)  The  components  aj^(mTc)  and  P^mT^);  j  =  I,Q;  is  I,Q;  are  mutually  uncoirelated 
which,  by  the  Gaussian  assumption,  renders  them  independent;  then  the  noises  n"(mTc) 
andn^(mTc)  are  also  independent 

(c)  Each  noise  porcess  has  mutually  uncorrelated  inphasc  and  quadrature  compcments 
(therefore,  by  the  Gaussian  assumption,  independent). 

(d)  Noises  corresponding  to  different  frequencies  (e.g.,  nJ(T)  and  r®  (T)  with  k  ^  i  ) 
are  approximately  uncorrelated 

(e)  Processes  nj(t)  and  nf(x),  however,  have  different  properties,  i.e.,  n^(x)  is  a 
highly  correlated  process  (as  a  function  of  x  ),  while  n^(x)  generally  is  not;  in  fact  samples 
of  no(x)  taken  T^  apart  could  be  uncorrelated. 

Of  the  above  conclusions,  (e)  is  probably  the  most  interesting  from  a  performance 
viewpoint:  If  the  noise  n§(x)  had  been  as  highly  correlated  as  n“(x),  the  decision  rule 

(5.44)  performance  would  have  been  perfect,  since  the  same^*  random  noise  sample  would 
be  obtained  at  x  =  0  and  x  =  T^;  thus,  they  would  cancel  out  This  not  being  true,  a  slight 

degradation  in  performance  is  expected,  as  was,  also,  observed  in  the  simulation. 
2*Same  within  a  totally  insignificant  change. 
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Nonetheless,  since  the  mean  part  of  the  interference  (second  term  in  (5.26))  does  cancel 
out,  the  performance  of  this  scheme  is  far  superior  to  that  of  the  radiometer,  which  is 
oblivious  to  that  term.  In  fact,  the  radiometer  output  is  merely  the  value  y(0),  which  is 
dominated  by  the  power  of  the  random  interference  (Figure  27).  Without  further 
assistance,  it  is  impossible  to  determine  if  there  is  any  signal  in  the  total  observed  power 
since  the  interference  contribution  is  random  and,  hence,  unknown. 

As  mentioned  above,  the  high  correlation  of  ot^ft)  implies  that  a^(0)  *  a^(Tg);  j  =  I,Q. 
On  the  other  hand,  the  degree  of  correlation  between  ^^(0)  and  Pk(Tc)  varies  with  the 
frequency  difference  Af^  =  fo  -  fj..  For  the  special  case,  wherein  one  of  the  interfering 
frequencies  fjj^  coincides  with  fo,  it  is  easily  seen  that  P]^0)  »  P^CTp);  j  =  I,Q.  Those 
facts  are  used  in  the  subsequent  analysis. 

We  can  now  turn  to  the  decision  rule  (5.80):  upon  squaring,  taking  the  difference, 
lowpass  filtering  (e.g.,  rejecting  double-frequency  terms),  assuming  that  f^  =  fofor  one 

frequency29  and  using  the  above  conclusions,  it  follows  that 


Tm 


(s)^hV¥  k?:  *  '^h(r)  (th  - 


where 


"For  the  case  of  a  single  interferer,  it  can  be  shown  that  worst-case  signal  detectability 
happens  when  the  interferer's  frequency  coincides  with  the  signal  frequency. 
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/  V*  ^ 

Lr<: **■ 


(5.85) 


-  I  n“(T)  +nJ(T) 


is  the  total  BP  equivalent  noise.  Let  us  note  that,  at  least  for  the  case  of  a  full  PN  code 
period  per  hop,  ycfT^)  is  approximately  zero,  which  considerably  simplifies  (5.84). 

Furthermore,  even  for  the  random-code  model,  a  comparison  between  analysis  and 
simulation  for  the  one-tone  case  (see  next  section)  has  indicated  that  the  impact  of  the  rv 
yc(Tg)  is  unnoticeable;  hence,  setting  yc(Tc)  =  0  seems  to  be  a  reasonable  approximation  for 

the  general  case.  Still,  a  full  analysis  of  (5.84)  without  further  simplifying  assumptions  is 
extremely  complicated.  In  order  to  gain  some  insight  here,  we  shall  focus  on  the  one-tone 
random  interference  (M  =  1),  with  the  reasonable  conjecture  that  the  multi-tone  case  should 
provide  analogous  conclusions. 

Since,  in  this  case,  Pi(Tc)  «  pj(0)  =  ai(0)  *  alfT^)  =  Ay  and  1)2(0)  =1)2(7^),  it 


follows  from  (5.84  )  that 


(5.86) 


where  we  have  defined  the  signal-to-interferencc  ratio  Yj  =  S/I  and  the  normalized  random 


I  ,  ^ 


-  *  . 


variable 


1 .norm 


-  *  1  /  \  ^ 

~  T~  ‘  '**1  S  *'n 


(5.87) 


In  (5.87),  c^;  n  =  1,...,N  are  the  code  chips  and  (>]  is  the  random  phase  of  the  interference. 
Clearly,  under  Hq,  A  =  0  with  probability  1 .  Under  Hj,  it  can  be  shown  that 
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MM 


(?\y 


1  +  2  Y, 


(S.SSa) 


var{^|Hi}  .  2CSTh)>Yj}-1 


(5.88b) 


Based  on  (5.88)  and  the  Gaussian  assumption  about  A,  the  perfonnance  of  this  decision 
scheme  is  predicted  by  the  detection  probability 


pcorr  ^ 

D  ^ 


-  ♦  V,-'  -  v) 


(5.89) 


where  Aq*  =  Ao/(STh)2  is  a  normalized  threshold  and  Q(x)  is  the  Gaussian  integral 
function.  In  the  absence  of  thermal  noise,  Aq*  can  be  set  arbitrarily  close  to  zero;  thus,  it 
always  yields  the  zero  false-alarm  rate  Pfa«0.  In  practice,  Aq*  would  be  set  according  to 
the  thennal  noise  level  and  the  degree  of  uncertainty  about  the  power  S  of  the  detected 

signal  .30 

The  performance  predicted  by  (5.89)  is  indeed  excellent.  We  note  that  increasing  the 
interference  power  (hence,  increasing  yj  *>)  actually  helps  detection  instead  of  deterring  it, 

while  it  has  no  effect  on  false  alarm.  For  instance,  it  can  be  shown  that,  if  N  >  60,  then 
>  99%  independently  of  Yi.  as  long  as  Aq*  ^  0.4,  i.e.,  one  can  tolerate  40% 

uncertainty  about  the  signal  power  and  still  expect  excellent  detection  capabilities, 
regardless  of  the  interference  power.  The  above  conclusion  is  rather  insensitive  to  N  in 
that,  for  N  =  10,  the  corresponding  minimum  P^is  90%. 


30Altcmatively,  this  can  be  expressed  in  terms  of  the  uncertainty  regarding  the 
transmitter/intercept  receiver's  true  distance. 


,4 

.  •  A  .V 


-Mv-'J 


•'-■wV.'-' 
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The  radiometer  performance  is  easier  to  analyze  for  an  arbitrary  M  and  is  based  on  the 
fact  that  its  output  Yj,,!  =  y(0)  can  be  written  as 

Let  Mniax  be  the  maximum  number  of  tones  which  can  be  expected  in  any  hop.  Here,  for 
simplicity,  we  consider  only  the  case  =  1,  where  Prob[one  interfering  tone]  =  Pr[no 
interfering  tone]  =  1/2  for  each  hop.  Assuming  that  the  threshold  y©  is  set  at  yo  =  IiTh,  so 
that  3'Pfa  =  0,  it  can  then  be  shown  that 


Thus,  the  obtainable  performance  decreases  with  decreasing  signal-to-interferencc  ratio  yj 
and  in  the  limit 

-  I  CUa  =  0 

(5.92) 

which  is  certainly  poor  compared  to  Pj®".  It  can  also  be  shown  that  the  radiometer 
performance  is  a  decreasing  function  of  so  that  the  above  results  constitute  an  upper 

bound  for  the  general  case.  Finally,  we  note  that  the  question  of  threshold  setting  is  much 
more  crucial  for  the  radiometer  than  the  correlator,  since  the  performance  of  the  latter  is 
effectively  independent  of  the  jamming  power. 

^^This  is  selected  so  as  to  match  the  zero  PpA  of  the  correlator:  another  choice  of  Yo  would 
lead  to  PpA  =  0.5. 


The  correlator  with  one  tone  interference  was  simulated  by  computer;  the  results 
shown  in  Figure  29  are  based  on  10,000  independent  trials.  Also,  shown  in  this  figure  are 
the  analytical  predictions  (dotted  lines)  whose  agreement  with  the  simulation  is  quite 
striking,  even  feu*  such  low  values  of  N  as  N  =3  and  N  =  10.  Since  detection  probability 
is  monotonically  increasing  with  N  (see  (5.89)),  those  values  represent  worst-case  designs 
which  nonetheless  yield  excellent  performance  results.  It  was  somewhat  surprising  to  find 
that  the  Gaussian  model  provides  such  an  accurate  analytical  prediction,  even  for  N  =  3; 
furthermore,  setting  ycCTj.)  equal  to  zero  proved  to  be  a  well-justified  simplification.  Note, 
that  no  false  alarm  was  observed  (Pp/^  =  0 )  and  that  performance  is  practically  insensitive 
to  the  amount  of  interference  inserted.  In  contrast,  the  radiometer  performance  (as 
evidenced  by  (5.91))  deteriorates  rapidly  with  decreasing  S/I,  as  expected.  Finally,  let  us 
mention  that  the  performance  shown  in  Figure  29  is  for  a  nonoptimized  (arbitrarily  chosen) 
threshold  Aq*  =  S^T^f^/N.  In  the  absence  of  thermal  noise,  further  improvement  can  be 
attained  for  the  correlator  by  decreasing  Ag*  to  a  very  small  (but  positive)  value. 

Although  analysis  and  simulation  are  not  yet  available  for  the  multitone  case,  it  is 
anticipated  that  the  gap  between  the  radiometer  and  the  correlator  performance  will  increase 
as  the  interference-to-noise  ratio  increases,  independently  of  That  gain  is  due  to  the 

relative  insensitivity  of  the  correlator  to  the  interference  nuisance  parameters. 

As  mentioned,  the  theory  of  this  section  does  not  include  thermal  noise.  The 
determination,  analytical  evaluation  and  simulation  of  algcrithms  which  operate 
satisfactorily  in  a  mixed  environment  (i.e.,  both  thermal  noise  and  random  interference), 
will  be  the  topic  of  future  research. 


Truncation  at  0.9999 


Real-Time  Autocorrelatcr  (Nonoptimized  Threshold) 
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APPENDIX  A 

ON  MOMENT-  GENERATING  FUNCTIONS  AND  LIKELIHOOD  RATIOS 


0 


f 


A  lot  of  powerful  results  in  statistical  decision  theory  rely  upon  the  judicious  use  of 
transform-domain  techniques  and,  in  particular,  on  characteristic  functions  (CF)  or 
moment-generating  functions  (MGF).  Given  the  CF  <I>x(0))  of  a  rv  X, 

OxCco)  =  (A.l) 


assuming  it  exists,  threshold-exceedance  probabilities  for  X  can  be  evaluated  via  the 
formula  [10] 


Pr[X>xo  ] 


(A.2) 


where  x^is  the  fixed  threshold  and  Im{»}  denotes  the  imaginary  part  of  the  enclosed 
quantity.  Expressions(A.l)  and(A.2)  are  particularly  amenable  to  numerical  evaluation  by 
appropriate  algorithms  (c.f.[l  1]).  For  the  purpose  of  proving  the  BFT  results  of  Section  2 
and  deriving  certain  Chemoff  bounds,  we  shall  find  it  convenient  to  deal  with  the  moment¬ 
generating  functions 


Let 


Mx(s)4tf{e«}  =  <I>x(-js) 


(A.3) 


denote  the  MGF  of  X,  where  s  is  a  real  variable;  in  fact,  it  is  typically  adequate  to  restrict 
attention  to  the  s  ^  0  range. 

When  the  averaging  in  (A.3)  is  over  the  pdf  (*)  under  hypothesis  Hj;  i  =  0,1, 
it  will  be  explicitly  indicated  as  a  subindex.  Furthermore, 
let 


H,(s)^lnM,(s) 


(A.4) 
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its  derivative.^  Whether  one  works  with  Hj(s)  (see,  for  instance,  [12])  or  g^fs),  as  we 
shall  do  here,  is  a  matter  of  notational  preference.  From  (A.4),  (A.5)  it  follows  that 

Mx(s)  =  exp  {  r  gx(s')ds'  +  Cg}  (A.6) 

^o 

Since,  however  M^CO)  =  1,  the  constant  Cg  =  0  in  (A.6),  so  that 

Mx(s)  =  exp{  r  g  (s')  ds'}  (A.7) 

^o 

which  is  used  extensively  below. 

Functions  [ix(s)  and  g^(s)  possess  some  interesting  properties.  For  instance,  the 
Taylor  series  expansion  of  4x(s)  around  s  =  0  has  the  form 

Ms)=  i  —  (A.8) 

B-O  n! 

where  C^is  the  mth  cumulant  of  X  (sec  ref.[10]).  Furthermore,  it  is  easily  shown  that 
m 


'The  dot  indicates  differentiation  with  respect  to  s. 


gx(0)  =  <?{X} 

gx(0)  =  var{x} 
gx(s)  2: 0 

i.e.  the  function  g^Cs)  is  monotonically  increasing. 

Let  X  be  some  decision  statistic  for  the  binary  test 


Hi 

X  ^  *0 
Hn 


(A.  10) 


(not  necessarily  an  LR  or  LLR),  with  the  corresponding  pdfs  shown  in  Figure 

A.l.  Then  the  performance  quantities 


PpA  -  Pr{X>XQ  Ihq}  .y*  fj^^j^(x)dx 


Pd  -  Pr{X>xo  ^\)  •  f  fx  |h,  Wdx 

•'xq  ’ 


can  be  bounded  as  follows  (Chemoff  bound):  If  the  r.v.  W  is  defined  as 

^  (O.ifx^xo 

(l,ifX>  Xq 

then  ,  for  every  s  ^  0, 


(A.  11  a) 


(A.l  lb) 


(A.  12) 


exp(-s(X  -Xq))  ^  W  ^  exp(s(X  -xq)) 


(A.  13) 


as  it  is  evident  from  Figure  A.l.  Furthermore,  from  (A.  1 1 ),  (A.  1 2), 


PFA=<f{W|Ho} 


(A.14a) 


t  $ 
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Pd=^{W|Hi} 


(A.  14b) 


From  (A.13),  (A.I4.a)  we  conclude  that 


« «  • 

Furthermore,  since  E(s)  =  g(s)  >  0  from  (A.9),  it  follows  that  the  extremum  is  actually  a 
minimum,  yielding  the  tightest  (or  lowest)  upper  bound  on  of  (A.  15), 


A  lower  bound  on  Pp  in  terms  of  of  gxiH^Cs)  can  also  be  derived  following  similar 
steps;  however,  its  most  utility  is  for  the  particular  case  of  interest,  whereby  X  stands  for 
some  LLR.  Then,  all  the  previous  conclusions  (which  of  course  still  hold)  can  be 
enriched  with  some  additional  ones,  owing  to  the  fact  that  the  MGF  Mj^iH^fs)  of  I  under  Hq 


which  links  the  properties  of  (s)  to  those  of  a  LR  K.  We  conclude  that 
Mt  Pi(a)  “  tfCeSi  pi)-  tf<A»  pi  1  - 


Mjt  P^(a)  -  P  (»♦') 


which  implies  that 


This  mteresting  aspect  of 


LLR’s  is,  also,  manifested  in  the  statement  of  the  PBF  theorem  (Section  2)  for  the  special 
Gaussian  case.  Since  M  £|h^(s)  is  just  a  "shifted"  replica  of  M  dhq(s),  the  latter  contains  all 
the  necessary  information  through  (A. 7),  it  suffices  to  know  go(s)  =  gjIn^Cs).  An 
immediate  corollary  of  (A.  19)  is  that 


(A.20) 


which,  in  turn,  yields  from  (A,7) 


f  So(s': 


go(s')  ds*  =  0 


(A.21) 


The  monotonic  increase  of  go(s)  (recall  (A.9))  along  with  the  fact  (A.21)  implies  that  go(s) 

should  look  something  like  the  curve  of  Figure  A.2,  since,  from  (A.9),  (A.  19), 
go(l)  ®  SjjH,  (0)  AlHj},  etc.  In  particular,  any  LLR  I  must  satisfy. 


(A.22) 


as  a  result  of  Figure  A.2.  Equalities  in  (A.22)  apply  only  to  the  trivial  case  where  A  =1 
almost  everywhere.  An  alternative  proof  of  (A.22)  can  be  found  in  ref  [13]  where 
tfltlHi);  i  =  0,  are  termed  directed  divergences  and  the  difference  J  = 

is  called  J-divergence. 

Another  byproduct  of  (A.  19)  is  a  lower  bound  for  in  (A.  1  lb)  in  a  form  dual  to 
(A.17).  Indeed,  from  (A.  14.b), 

Pd  ■  <?{W  |Hi  }>  ~  Xq)  |h^} 

-  e-sxQ  Mx  ki,  (-S)  -  e“«*o  Mx|h^(1-s) 

or,  changing  the  variable  s  into  (1-s)  without  loss.^ 


Pp  >  )  0  •  exp((i-s)xo  ♦  J  go(s*)dsM  (a,23) 

0 

Again,  by  differentiation  of  the  exponent,  the  extremum  (which  is  easily  shown  to  be 
a  maximum)  occurs  at  s  =  s*  where  go(s*)  =  xq,  resulting  in  the  tightest  (highest,  lower 
bound) 


or 


Pd>  exp{(t-s*)  go(3») 


♦ 


3 

go  (  ®'  )**3'  ) 


(A.24a) 


i  Since  we  are  restricting  attention  to  the  0  ^  s^  1  range  and  0  S  1  -  s  ^  1,  the  substitution 
is  legitimate 


PdI^Pd^LB  “  «*0(Pfa)uB 


(A.24b) 


with  (Pfa)  ub  as  per  (A.  17). 


Equipped  with  this  brief  expose  on  the  useful  properties  of  the  LLR  moment¬ 
generating  functions,  we  can  now  prove  the  BPF  theorem  in  a  rather  compact  fashion. 
Indeed,  assume  that  i  is  Gaussian  under  Hg,  with  ^{i|H  }  =  mg  and  var{i|Ho}  =  its 

second-order  statistics.  Since  the  MGF  of  a  Gaussian  r.v.  is  well-known  to  be. 


”  -  exp  (mo  a  ♦  -  } 

2 


(A.25) 


it  follows  that  the  associated  go(s)  is  simply 


go(s)  (In  M  jth^(s))  -  mo  ♦  o^a 


(A.26) 


In  other  words,  go(s)  is  a  straight  line  as  shown  in  Figure  A.2  with  a  constant  slope 


80 (s)  -  o2 


From  (A.20),  (A.25)  we  have  that 


var{i  fig)  ■  var{l  1 


(A.27) 


^  1  IHo  ^  ^  "  explmo  ♦ 


}  .  1 


2  (A.28) 

i.e.  the  mean  and  variance  of  I  under  Hq  cannot  be  arbitrary  within  the  Gaussian  model, 
but  rather  related  as  per  (A.28).  Inserting  (A.28)  into  (A.25)  results  in 


0^ 

“  expl -  {s-i)sl 

^  (A.29) 


We  can  now  use  (A.  19)  to  assess  that 


o2  o2 

llH<  “  «*P  {“■“  a  ♦  s^l 
^  2  2 

(A.30) 


But  this  is  exactly  the  MGF  of  a  Gaussian  i  under  Hj,  with  var(l  IHj}  =  a|  (in  accordance 
with  (A.27))  and 


o2 

S{%  iHil  -  -  - 

^  (A.31) 


as  shown  in  Figure  A.2.  We  conclude  from  (A.27),  (A.31)  that  the  variance  var{i  IHq}  = 
describes  fully  the  Gaussian  pdfs  of  t  under  either  hypothesis.  In  order  to  complete  the 
statement  of  the  PBF  theorem,  let  us  note  that,  for  any  LLR  t , 


”  ^{a2  ho)  -  1  ♦  varlA  J 

0 


so  that  the  "Gaussian"  distance  d^,  defined  as 


-  ln(i  ♦  var{A  ) )  “  1”  M  (2) 

However,  for  the  particular  Gaussian  model,  (A.29)  yields 

M  IMo  ^  “  ®*P  ^ 
which,  when  compared  with  (A.32)  results  in 


(A.32) 


(A.33) 


0^  «  >  ln(t  ♦  var{A  ^Iq  ^  ^ 


(A.34) 


Note  that  the  LR  A=  e*-  in  (A.34)  is  obviously  not  Gaussian. 

The  content  of  this  Appendix  could  provide  the  basis  for  exploring  the  propoiies  of 
LLR's  beyond  what  is  presented  here. 


I 


APPENDIX  B 


PROOF  OF  EQUATION  (3.7) 


Let  us  start  from  (3.4),  i.c.. 


(r(t))  = 


I" 

i=l  '^0 


,N 


N  )  «1 


(B.l) 


with  as  per  (3.5),  Cy  denoting  the  j*  chip  of  the  i*  code  pattern  and  rj  as  per  (3.6). 
Consider  the  N-dimensional  linear  space  consisting  of  all  binary  patterns  £;  ;  i=l,...,2N  of 

length  N.  Now  observe  that  this  space  can  be  divided  into  two  disjoint  complementary 
subsets,  each  of  cardinality  2^-1,  by  the  rule  that,  for  each  possible  vector  belonging  to 
one  of  the  subsets,  belongs  to  the  other  subset.  Exactly  which  vectors  Cj  are  included  in 
which  subset  is  immaterial  as  long  as  the  above  rule  holds.  We  can  now  perform  the 
summation  in  (B.l)  over  one  of  the  two  complementary  subsets  instead  of  the  whole  linear 
space,  with  the  equivalent  rule. 


1»1 

ie<g 


(B.2) 


where  i  =  (r^,...,rfj),  ^  is  one  of  the  aforementioned  disjoint  subsets  and  any  resulting 
scale  factors  are  absorbed  into  the  threshold  A’q.  In  (B.2),  cosh(x)  is  the  hyperbolic 
cosine  function 


cosh(x) 


Since  cosh(x)  is  an  even  function  of  x,  the  summation  in  (B.2)  can  be  extended  over  the 
whole  linear  space 


2" 

I 

i«l 


(B.4) 


In  order  to  proceed,  we  first  need  to  introduce  the  superscripts  i^and  which  will 
indicate  the  length  of  the  corresponding  vectors,  and  then  prove  the  following; 


T^mma  1.  It  is  true  that; 


(B.5) 


Proof:  The  proof  is  by  induction.  First,  let  N  *  1.  Then, 

2 

^  cosh  =  cosh^r^]  +cosh  *  2cosh^rjj 


(B.6) 

since  cosh  (x)  is  even  function.  Thus,  (B.5)  is  satisfied  for  N  =  1.  Next,  assume  that  it  is 
true  for  N  =  k,  i.e.. 


I  cosh  •  Z'^fl  cosh(rj) 

i-1  "  ■J'l  (B.7) 


and  prove  it  for  N  -  k  +  1 .  Indeed, 


M 

1 


r 


,k+l 


cosh  j^oshlr^’-si''  ♦  r^,i)  ♦  cosh 


(B.8) 


But, 


cosh(a+b)  +cosh(a-b)  * 


,a^b^^-a-b^e^-b^^-a4b 


=  e 


,  0-  (i^)  • 


(,s  ,  ,-0)  (sl^) 


=  2  cosh(a)cosh(b) 


(B.9) 


which,  upon  substitution  in  (B.8),  yields 


cosh(/'^*".o'‘^*»)  -  £  2  cosh(r<'‘'.c<'‘>)  cosh(r,,i) 

.  2  cosh  cosh{r^'‘'.c5'‘') 

.  2  cosh  (Vi)  2'‘  ^  cosh  (fj)  -  Z"*' 


(B.IO) 

where  the  troth  of  (B.8)  has  been  used.  Equation  (B.IO)  completes  the  proof  of  Umma  1. 
If  Lemma  1  is  applied  to  (B.4),  it  renders  it  equivalent  to 


N 

"fj"  cosh 

j-1 


A’, 


(B.ll) 


which  is  (3.7a). 

< 
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APPENDIX  C 


PRCMDF  OF  EQUATIONS  (3.9)  -  (3.10) 


Let  us  rewrite  the  decision  rule  (3.8)  as 


X 


»1 

»0 


An 


(C.i) 


where  the  independent,  identically  distributed  random  variables  (iid  rv's)  yj  =  assume 
the  form 


\  nj2  under  Hq 

|(VST<.  Cj  +  np2  under  Hj 


(C.2) 


In  (C.2),  nj  represents  the  Gaussian  rv 


"J 


j  ^  nj(t)  dt 
(M)T, 


(C.3) 


whose  statistics  are  independent  of  j. 

A  precise  analysis  should  account  for  the  exact  distribution  of  yj  under  either 
hypothesis.  So,  for  instance,  under  Hq  yj  is  chi-squared  rv  with  N  degrees  of  freedom. 
For  large  N  (of  the  order  of  hundreds  or  more,  that  is  the  practical  case),  very  satisfactory 
approximate  results  can  be  obtained  which  circumvent  the  difficulties  of  an  exact,  but 
enormously  complicated,  analysis .  This  is  done  by  invoking  a  central-limit-type  ai;gument; 
The  decision  statistic  X  being  the  sum  of  a  large  number  of  iid  rv's  yj,  is  approximately 

Gaussian  distributed,  with  mean 


<C.4a) 
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and  variance 


varjxlH^J  »  N  varjy.jH^j  ;  k=0.1 

(C.4b) 

under  either  hypothesis.  Thus,  having  the  second-order  statistics  of  the  individual  rv's  yj 
and,  subsequently,  of  X  from  (C.l)  enables  us  to  determine  performance  immediately. 

Let 

m  =  VST,Cj5ki  (C.5) 


so  that 


yj  =  (m+nj)2  =  m^  +  nj^  +  2mnj 


(C.6) 


Furthermore,  the  Gaussian  rv  nj  has  zero  mean  and  variance 

2  _  ^t.2) 


Thus,  from  (C.6), 


Vc 


■H 


2  2 
m  + 

n 


2  5  '*0  ’’’c 

*kl  ^  “1^ 


(C.7) 


i«0\)  •  ("oTc)(^*V«) 


(C.8) 


Finally, 


var{yj}  =  var  {m2+nj2+2mnj}  =  var{nj2+2mnj}  =  var{nj2}  +  4m2  var{nj}  (C.9) 

since  the  random  variables  nj^  and  nj  are  uncorrelated  (^{nj^}  =  0).  Thus, 


f 

r* 

[.■ 

I  ^  var{yj}  =  (3o„4  -  a/)+4m^a2  «  2o„2[o„2  +  2m2l  =  (NoTc)2[l/2  +  2yAi1  (C.IO) 

Equation  (3.9)results  from  combining  (C.8)  and  (CIO)  with  (C.4). 

i- 


APPENDIX  D 


PERFORMANCE  OF  RULE  (3.19) 

A  Gaussian-approximation-based  analysis  of  the  low-SNR,  chip-noncoherent  detector 
(3.19)  can  be  performed  following  guidelines  similar  to  those  in  Appendix  C,  as  follows: 
Let  the  bandpass  AWGN  n(t)  be  represented  by 

n(t)  =  ^2  I  nj  (t)  cos  ojQt  -  nq  (t)  sin  ©ot  |  (D.  1 ) 

where  nj  (t)  and  nq  (t)  are  baseband,  independent,  Gaussian  processes  with  a  flat  one-sided 
PSD  of  Nq  W/Hz  each.  Conditioned  on  the  unknown  phase  (l)j  and  the  ±l-chip  Cj  of  the 
interval  (j-1)  ^  t  >  jT^  the  inphase  and  quadrature  variables  rj  and  rq  of  (3.15c)  are 

Gaussian  with  means 


k=0,l 


and  common  variance 


(D.2) 


(D.3) 

Thus,  averaging  the  phase  ^  <{ij,  it  follows  that  each  of  the  independent  envelope  rv’s  rj  of 
(3. 1 8)  has,  under  H] ,  a  Rician  distribution 


^Again  it  is  easily  seen  that  the  value  of  Cj  becomes  irrelevant  due  to  the  squaring. 


where  is  the  predetection  SNR  while,  under  Hq  (i.e.,  signal  absent,  Yc  “  ®)» 
reduces  to  a  Rayleigh  distribution 


(D.5) 

Therefore,  the  statistics  of  the  decision  rv  X  in  (3.19),  being  the  sum  of  N  iid  rv's,  can  be 
precisely  found  (see  (3.21  a,b)  in  text). 

Although  those  results  are  exact,  they  are  cumbersome  to  use  from  a  computational 
viewpoint,  especially  for  large  N  (N  >  100).  Besides,  the  Gaussian  approximation  about  X 
becomes  suffficiently  tight  in  that  region  and,  thus,  emerges  as  an  attractive  and  simple 
tool. 

Let  us  start  from  the  fact  that 

"ij  ‘  >*1;  *"ij 

(D.6a) 

and 

rg  =  sin* 

''  p.6b) 


In  (D.9),  the  Gaussian  rv's  n| .  and  ng .  are  defined  by 


J  nj(t)  dt 
(j-l)T, 


(DJa) 


’q.  =  J  "gCt)  dt 

^  (j-i)Tj, 


and  are,  therefore,  mutually  independent,  each  with  zero  mean  and  variance 


(D.7b) 


Then,  the  squared  envelope  rj  is 


■  'U  *  ’■q.^  =  *  „  2  *  n,  2  *  2/5  Cj  ♦nqj  six 

J  J  J  J 


It  immediately  follows  that 


‘ki  *  2'’.'  •  sTc'\i*Vc  •  (VcK>*-'c‘u) 


(D.9) 


which  is  (3.22d) 


In  order  to  calculate  the  variance  of  rj  in  (D.8),  we  note  that  the  first  term  is  a  constant 

(thus,  it  can  be  neglected),  while  the  second,  third  and  fourth  terms  are  pairwise 
uncorrelated.  This  is  because:  (a)  n,  2  and  ng  2  are  functions  of  independent  rv’s  and,  (b) 
the  fourth  term  includes  both  cj  and  <|>j ,  which  are  independent  of  the  noise  and  zero  mean. 


As  a  conclusion. 


where  noises  n^  and  nj+j  j,  are  mutally  independent,  as  are  Cj  and  Cj+i. 


Thus, 


4j2''J*1.i('  4t)  *kl  4j^j*l}*-^{"j2Vl.l}  ■  ° 

In  an  analogous  way. 


r ^ 
Jl 


«lcl  ^  "jl 


(E.5) 


which,  combined  with  (E.3)  yields 


ST  2 

ST  2 

< 

4ji''j4 

-  —  \1  ^  <^{"jl  "j2f  = 

—  % 
4 

(E.6) 


by  virtue  of  independence  between  nji  and  nj2.  Finally,  substituting  the  above  into  (3.41) 

|Kk  •.  -Ol-  "If  4jl1  *  4jl'-j2lj 


("o^c)  ■  ■  5T 


(>  *  \  \l)  *  -f-  \l 

(4^)[  *  (4)  Vc  «k.] 


(E.7) 


Before  we  proceed  with  var  {X2I  Ho),  let  us  for  a  minute  turn  to  the  worst-case 

mean  of  (3.41)  as  explained  in  the  report  However,  the  picture  now  is  changed  to 
Figure  E.2,  under  the  assumption  of  a  worst  epoch  e  =  1/4.  Then,  we  can  write  that 


■■^2  ■  ^(f)  ‘j  *kl  *  "j2 


(E.8) 
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and 


Note  that  rj2is  identical  to  (E.3);  thus,  e  =  1/4}  is  the  same  as  above. 


Futhermore 


(E.10) 


where  the  conditioning  events  Ajgnjg  and  Ajjff  refer  to  Cj.i  and  Cj  having  the  same  or  a 
different  sign.  Clearly, 


^rt^same)  ®  ®  1/2 


(E.11) 


a  fact  used  in  (E.  10).  Let  us  now  note  that,  under  A^iff,  the  mean  of  rjj  is  zero  while,  under 
AjanKJ  is  equal  to  VS  Tj,/2.  Thus, 

•  (7)(-V)(l  *  ic  ‘kl)  '  ]i  *  7  \  *ki] 


(E.12) 


In  exactly  the  same  way, 


*\r.,  r 


"...i/*}  -  ..j))*.. '.{"jl  "ji} 

■  •  r^teh, 


(E.13) 


which  is  identical  to  S  {rj2’j+i,i  I  H|,;  e  *  1/4), 
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Put  together,  the  above  result  in 


«{X2|Hk;E=l/4}  =  N 


(E.14) 


In  both  cases,  var  is  the  same.  A  straight-forward  calculation  can  show  it 

to  be 


var{X2(Ho}  =  N  (^q\)2 


(E.15) 


Combining  (E.7)  and  (E.14)  with  (E.15)  results  in  (3.42)  and  (3.43),  respectively. 


~r  ^  r  A 


APPENDIX  F 


I  J 


ANALYSIS  OF  THE  RADIOMETRIC  DETECTOR 

Here  we  briefly  develop  the  LPI  performance  of  the  radiometer,  as  shown  in 
Figure  9,  when  the  input  r(t)  is  as  per  (3.2).  The  output  rgp(t)  of  the  rectangular  BP 


filter  is 


rgp(t)  =  ¥'7$c{t)  cos(wQt  +  (j()  +  f>3p(t) 


where  n^pCt)  is  filtered  AGN  with  the  typical  BP  rep-esentation 


ngp(t)  =  /2  nj(t)  cos  uq!  +  *1^(1)  sin 


V*  - 


In  (F.2),  njCt)  and  nQ(t)  are  independent  zero-mean  jointly  Gaussian  low-pass  noise 
processes,  each  with  a  flat  PSD  of  Nq  W/Hz  (one-sided)  and  bandwidth  of  Wgp/2  Hz  (one¬ 
sided).  The  code  signal,  filtered  by  the  BPF,  is  designated  c(t)  in  (F.l).  The  zonal  LPF 
(shown  for  analytical  purposes  only)  rejects  the  2fo  components  of  the  square-law  output 
z(t)  and  has  an  output  given  by 

2Lp(t)  =  s(c(t))^  +  nj^(t)  +  np^(t)  +  2^c(t)  n^Ct) 


In  order  to  maintain  analytical  tractability,  we  approximate  and  model  c(t)  as  an 
attenuated,  but  undistorted,  version  of  c(t),  namely. 


c(t)  ^  ac(t) 


t  4 


A  .'-V- 
‘AV.’ 


J 
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where  the  "attenuation  factor"  measures  the  (normalized)  power  reduction  due  to 
filtering 


Z  5  X  |H|.(f)|^  S^(f)«lf 

*CD 


■d) 

0 


(F.5) 

since,  for  a  ±1 -valued  code,  the  denominator  of  (F.5)  equals  c2(t)  =  1,  In  the  previous 
expression,  Sc(f)  is  the  PSD  of  the  random  code  signal  c(t)  and  HL(f)  is  the  low-pass 
equivalent  of  the  BPF  HBp(f),  i.e.. 


HL(f) 


‘  |f|  <«Bp/2 

0  ,  otherwise 


(F.6) 


By  suitably  adjusting  Hl  (f).  (F-5)  can  be  generalized  to  filters  other  than  the  perfectly 
rectangular  one  considered  above. 

From  the  above  definitions,  it  follows  that  (c(t))2  = 

Thus, 


Zj^p(t)  =  a^S  +  nj^(t)  +  nq(t)  +  Z/S’ocCt)  nj(t) 

(F.7) 

under  Hj  Under  the  alternative  Hq,  is  given  by  the  above  expression  with  S  =  0.  It 
follows  that 


-  o^S  6.,,  +  N,,  k-0,1 


with  as  per  (3.9c) .  Thus,  the  expected  value  of  E{X|H^}  of  test  statistic  X  at  the 
output  of  the  integrator  in  Figure  9  can  be  calculated  under  either  hypotheses  k  =  0,1, 
as 

■  f  k-O.l 


The  next  step  is  to  obtain  the  second-order  statistics  of  Zy>(t)  and  X.  The  typical 
approach,  based  on  the  assumption  WgpNTg  »  1,  is  to  model  ZLp(t)  as  a  very  wide-band 
("delta-correlated")  process  with  respect  to  the  approximate  bandwidth  (NTg)*i  of  the 
integrator.  Thus,  one  needs  to  evaluate  only  the  PSD  of  ZLp(t)  at  the  origin  f  =  0, 

distinguishing  between  the  average  (constant)  contribution  and  the  random  contribution. 
Using  the  fact  that  ni(t)  and  nQ(t)  are  Gaussian,  the  autocorrelation  function  of  Zu>(t)  can  be 
obtained  directly,  from  which  its  PSD  under  Hj  is  given  by 

2 

hip(f)  *  {ah  +  NpWgp)  «(f)  4  4Sni(f)  ®  SnQ(f)  +  4a^  SSj.(f)  @  Sni(f) 

(F.IO) 

where  ©means  "convolved  with."  The  coefficient  of  the  Dirac  delta  function  6(f)  agrees 
with  the  first-order  statistics  of  Zjj)  (t)  obtained  directly  in  (F.7).  The  two-sided  PSD  of 
zy)(t)  at  f  =  0  measures  its  random  contribution  and  is  obtained  by  evaluating  the  last  two 


terms  S7  (f)  in  (F.9)  with  f  =  0,  resulting  in 
LP 


S,  (f) 

LP 


=  Wgp  4  2o^  SN 


0  ;  ic=o,i 


(F.ll) 


The  delta-correlated  random  part  of  z^(t)  ccmtiibutes  to  the  variance  of  X  as 
v.r{A|H,}  .  *  Za*  SNo  ^l) 


•  to  '  -  - 

-j 


-'T.-n 


(F.12) 
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The  mean  and  variance  of  X  in  (F.8)  and  ^.1 1)  axftice  to  characterize  performance  under 
the  Gaussian  assumption  which  hinges  on  the  fact  that  WgpNTg  «  N»  1.  The  result  is 

(3.10a)  with 


a 


(F.13) 


Clearly,  the  choice  of  Wgp  has  an  impact  on  a  through  (F.12)  and  in  (F.5).  The 


quantity 


a 


0 _ _ 

V^BP  '•’c 


(PM) 

has  been  plotted  in  Figure  F.l  as  a  function  of  the  product  W^pTg  The  maximum 
ttniax  “  ^  ^  attained  at  WgpTc  «  1. 
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APPENDIX  G 


DERIVATION  OF  EQUATION  (3.50)-(3.52) 


Let 


A= 


N  H, 

I  I  '0 


j=l 


be  the  decision  rule  (3.48),  where 


yj  =  r(i)  rj(2) 


and  rjW  ;  i  =  1,2  is  given  by  (3.49).  Then, 


r/0  =  ^/STcCj6l.l+n/i) 


(G.l) 


(G.2) 


(G..3) 


where  nj(i)  and  nj(2)  are  zero  mean,  mutually  independent  Gaussian  rv's  with  common 
variance  0^2  =  NoTc/2.  From  (G.2)  and  (G.3),  it  follows  that 


(2) 


(G.4) 


Thus, 


'{'j}  ■  -  (KoT^)y, 


(G.5) 


which  justifies  (3.50).  Futhermore 


APPENDIX  H 


DERIVATION  OF  (4.13) 


LctR  =  in  (4.10)  and  (4.11).  Then,  following  steps  identical  to  (4.7)  yields 


(4.10) 


1  + 


2^  _\ll 


where  the  density  of  R  under  is  the  Rician 


/(R1Wi)-7j  exp{-Yw}exp{ — )/o 

i  ff. 


1  »} 


2^ 


with 


•  r 


and  Yh  =  ST^/Ng.  We  note  further  that 


(H.1) 


(H.2a) 


(H.2b) 


iJs 

A(ii)-«p{-r„}/J— — 


is  an  LR.  Thus,  from  the  known  properties  of  LR's 

2^ 


— 


Ho 


rnglm 


and 


A[R)f{R\H^)dR. 


(H.3) 


(H.4) 


(H.5 ) 
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r\3^  "A  "ji "A p^KU Piju *ia 


Substituting  (H.2),  (H.3)  into  (H.5)  and  letting  x  *  R/o„  yields 


<f{A*|«,}-exp(-2Tr,)j^xexp|- Yj/o*[y2Y;^x]  dx. 


(H.6) 


From  [21,  pp.  718,  6.633.2]. 


pe-*'’'j,{ax)J,(fix)xdx 


(H.7) 

Letting  p2  =  1/2,  a  =  P  =  >/2yh.  P  =  0  (here  j  =  and  using  the  fact  that  Io(x)  =  JoC^t) 
(H.6)  and  (H.7)  results  in 


/j| 


-«xp{2r,r)/,(2y^) 


(H.8) 


Combining  (H.l),  (H.4)  and  (H.8)  provides  (4.12).  As  a  byproduct,  we  note  from 
(H.3),  (H.5),  and  (H.8)  that 


■«*P{T»)/o(2yw) 


(H.9) 


•  *1^  • 


■» 


lit 

I .  - 


*v^-- 


S-‘w: 


t- 


•V.'.v 


» •  *  •  I 
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APPENDIX  I 


EVALUATION  OF  FOR  CASE  H 


We  consider  here  the  Case  n,  where  no  power  adjustment  talres  place  in  die  decision 
rule  (i.e.,  (S.8)  is  used),  but  the  actual  frequency  offset  of  the  incoming  signal  with  respect 
to  the  observed  spectral  points  is  the  worst,  namely  Af  =  H.  The  case  is  depicted  in 

Figure  1. 1.  Note  that,  in  principle,  a  residual  signal  power  can  be  measured  in  frequencies 
beyond  the  two  adjacent  ones  (f^,  f^^i),  but  that  amount  is  insignificant  and  can  be 

neglected.  On  the  other  hand,  each  of  the  adjacent  spectral  slots  observes  a  tone  of  power 
aS  (  a  as  per  (S.13))in  noise,  so  that  the  the  probability  density  function  (pdf)  of  the 
envelope  R  in  each  adjacent  slot  is  the  Rician  pdf 


where  Yh  is  the  hop  SNR  and  =  NoTh/2  is  the  noise  variance  in  each  quadrature 


component  of  the  bandpass  additive  Gaussian  noise.  Since,  in  the  absence  of  the  signal, 
the  mean  and  variance  of  Io(2VS  R/Ng)  remain  as  per  (5.  lla,b),  we  only  need  to  evaluate 


( 
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The  latter  integral  can  be  evaluated  by  prc^>er  modification  of  the  expression  in  [21,  p.718], 
yielding 


(  /  2/S  I  signal  \)  “atHf  (I^oIyS 

—  Rj  present^!  -  e  ^loCZ/a  Yh)®  J 

-  e  Yh^ 

a.4) 

Fianlly,  noting  from  Figure  (1. 1)  that  the  expression  (1.4)  pertains  to  2  slots  (under  Hi), 
while  the  remaining  (G-2)  are  only  noise,  yields 


,r{A|Hi;Case  11}  -  2e  Io(2/o  Yh)  ♦  (0-2)« 


instead  of  (5.1  la).  Putting  (1.5),  (5.1b),  (5.10c)  and  d2  together  results  in 


d2 

A 


(Io(2/a  Yh)-1)2 
Iq^^Yh) "1 


as  per  Table  1  in  Section  5.2. 

A  final  note  is  in  regard  to  the  notation  in  (I.l)-(1.4):  we  condition  on  signal 
present",  as  opposed  to  "Hi",  because  the  above  expectations  pertain  to  the  slots  which 
truly  contain  a  signal  component,  while  "Hi"  simply  means  "signal  in  some  slots  in  the 


bandwidth  B". 


APPENDIX  J 


STATISTICAL  CHARACTERIZATION  OF  NOISE  PROCESSES 
Ni'iDANDNqIt^ 

We  define  Nj^t)  and  Njj^x)  as  per  (5.45)  and  (5.40).  It  then  follows  from  (5.41), 

eq 

(5.42)  and  (5.43)  that  Ni(Xk),  NQ(Xk);  \  =  KB**  are  zero-mean  random  variables,  i.e., 
(5.46).  Furthermore,  by  virtue  of  a  central-limit-type  argument  (summations  of  a  large 
number  of  rv's),  they  are  assumed  to  be  approximately  Gaussian. 

We  shall  first  show  that  P^^x^)  and  NQ(Xj)  are  uncorrelated  for  every  k  and  j;  in 

which  case,  they  are  also  approximately  independent  (due  to  the  Gaussian  assumption). 
Indeed,  from  (5.45a). 

D}  ”iq^P} 

(J.i) 

Let  us  consider  the  first  term  on  the  right-hand  side  of  (J.I);  a  similar  line  of  argument 
applies  to  the  next  three  terms.  By  definition, 

'{"iiCtkD 

■‘k  ^k 

(L2) 

However,  nQ(t)  is  independent  of  ni(t),  which  implies  that  the  integrand  of  (J.2)  contains 
the  factor  ^{nQ(t2-Xj)},  which  is  zero.  Thus, 


Along  with  (3.1)  and  the  above  remarks,  (J.3)  establishes  the  uncorrelatedness  of 
eq  eq 
Nj(Xk)  and  NqCTj) 

eq  eq 

Let  us  now  consider  Nj  and  Nj(x^^)  for  x^  #  Xic^.  We  have  that 


•I-  « 


{"iiC'ki]''qQ0k2D}  *  *NqC'ki)"uC'i<23} 

(J.4) 


But, 


■  /V "  ^{"iCti>iCti  -  ‘rki>iC^2>iC^2  -  •nt2D} 


■'k  ^^k 

Ki  K2 


^ki  "k2 


(J.5) 


We  shall  approximate  the  integral  (3.5)  as  the  double  sum 


(J.6) 

where  Ax  =  Xj  =  B*i.  We  note  that  nj(k)  and  nj(k)  and  niOc-kj)  are  independent  rv’s  since 
ki  +  0,  as  are  n^fm)  and  ni(m-k2).  Furthermore, due  to  the  independence  of  the  rv's 
involved,  the  expectation  in  (3.6)  will  be  zero  whenever  all  four  sampling  times  k,  k-kj, 
m  ,  m-k2  are  pairwise  different,  as  shown  in  Figure  3. 1(a).  If  we  assume  that  ki<k2,  we 
then  observe  that  the  three  remaining  cases  depicted  in  Figure  3.1  are  (b)  k-kj  =  m-k2;  in 
which  case,  k#m,  or  (c)k  =  m;in  which  case,  k-ki#m-k2,  or  (d)  k=m-k2:  in  which  case, 
k-kj+m.  In  all  those  cases,  the  common  conclusion  is  that  at  least  two  rv’s  exist  which  do 

not  coincide  with  the  others.  The  net  result  is  that  the  expectation  in  (3.6)  will  always  be 
zero,  establishing  the  approximate  uncorrelatedness  of  Nii(xki)  and  Nii(xk2).  The  same 
exact  argument  establishes  that^{NQQ(xk, )  Nii(xk2)}  *  0,  while  it  is  even  easier  to 


show  that  NQQ(]T|j23|a  s  0.  Thus  all  of  the  above  can 

be  combined  in  (J.4)  to  illustrate  that  different  samples  (^i  of  the  same  process  N*  we 

approximately  uncorrelated.  It  can  similarly  be  demonstrated  that  the  same  is  true  for 
Nq(t),  defined  in  (5.45c). 

eq 

Finally,  we  are  concerned  with  the  variance  of  N  j  (Tjj)  and  NqCx^).  Clearly,  since 
the  respective  means  are  zero,  those  variances  coincide  with  the  second  moments. 
Furthermore, 


(J.7) 


But, 


r  P  ‘^!"iCtO''iOr-'iO"iCt2>iC^2-TO[dt^dt2 

■^k 


(J.8) 

where  a  simple  change  of  variables  has  been  performed.  We  shall  now  use  a  familiar 

4 

property  of  four  jointly  Gaussian  rv's  {xjj.i,  namely. 


<?|xjX2}  ^1x3X4}  +^{xiX3},r|x2X4}  +<?{xiX4|^{x2X3| 


(J.9) 


Applying  (J.9)  into  (J.8)  results  in  (for  any  t): 


H 

*{"iC‘2*0"iC‘2)}  <“i"2 

Th-t 

0 

‘^\p 


T„-t 

4.  /*  /*  d2  _4-** 


•  «nn(-)CVO  W-t'i)  dt^dt^ 

*  jj  'n,£»r‘2-0 

0  (J.IO) 

where  RnnW  =  #{tij(t)nj(t+T)}  is  the  correlation  function  of  nj(t). 

Equation  (J.lO)  can  be  further  simplified.  First,  we  note  that,  for  x  =  =  kB*i, 

Rfln  (“^k)  =  so  the  first  term  drops  out.  Furthermore,  one  can  employ  the  even  symmetry 
of  Rnn(^)»  i  c.,  the  fact  that  Rnn('C)  =  Rnn('‘')»  'M'der  to  reduce  the  two-dimensional 
integrals  to  one-dimensional  ones  (see  [49,  page  325]  for  details).  The  result  is  then 

0 


(J.ll) 


With  the  change  of  variables 


■'k 


ai2) 


and  defining  the  normalized  parameter  as 
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we  can  rewrite  (J.ll)  as 


/ '  o-o'  )*■  E®("  '(^H-’0*’k)] 


-°°-  (;F^(k)  t  F2(k)) 


(JM) 


where 
F, 


l(k)  =  (l-p')Sa^  -^p*  dp'  =  y*^  (l-p-)Sa^[uGCl-cOp]dp 


(J.15) 


and 

Fo 


(k)  =  ySi-P')  Saj^G(p'(;i-fO  *  *  ^kE*^**' 


(J.16) 


Combining  (J.14)  with  (J.7)  yields  (5.47a). 
Quite  similarily,  it  follows  from  (5.45c). 


Th 

ff  <^{"iC‘i)"QC‘rO"iC‘2)"qC‘2-o}‘“i‘“2 

Th  - 

yj  "nnCM-tijdtjdtj  .  yy  ii‘„cti-ta  dt^dtj 


(J.18) 

since  ni(x)  and  ngCx)  are  independent,  with  the  same  correlation  function.  The  integral  in 
(J.18)  is  identical  to  the  second  term  in  (J.IO),  which  corresponds  to  the  Fi(k)  term  in 
(J.14).  Furthermore, 


Th 

f{NiQ(T)  Nqj(t)}  •yy^{niCt0np(;ti-T)np(;t2)  np^tg-x^  dt 


T 

Th-t 


jf  """(‘i '  ‘2  *  ’)''nn(‘i  -  4  -  ')''4><4 


(J.19) 

which  is  identical  to  the  third  term  in  (J.IO)  corresponding  to  the  F2(k)  term  in  (J.14). 
Combining  the  above  with  (J.17)  (note  the  minus  sign  of  the  second  term)  results  in 


as  per  (5.47b) 


(J-20) 


4 
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APPENDIX  K 


DETECTION  OF  FREQUENCY-HOPPING  SIGNALS 

VIA 

AUTOREGRESSIVE  MODELING 

KI.  INTRODUCTION 

The  purpose  of  this  Appendix  is  to  present  a  number  of  methodologies  that  can  be 
used  for  the  detection  of  frequency-hopping  signals  in  AWGN.  These  methods  operate  in 
the  autocorrelation-domain  (ACD)  and  are  based  on  autoregressive  (AR)  power  spectrum 
estimation  concepts.  They  are  simple  to  implement  and  have  low  computational 
complexity,  i.e.,  they  have  properties  that  make  them  very  suitable  for  real-time 
implementations.  Specifically,  four  different  methods  will  be  discussed  and  their 
performance  demonstrated  and  compared  with  that  of  the  conventional  radiometer  and  the 
Correlator-Detector  [see  section  5.4]. 

The  organization  of  the  Appendix  is  as  follows.  Section  KII  establishes  the 
problem  formulation  and  its  imposed  constraints.  Section  Kill  presents  the  four  different 
detection  methods,which  are  all  based  on  AR  models.  The  discussion  emphasizes  the 
methods  mathematical  description,  performance  evaluation  and  computational  complexity. 
Concluding  remarks  and  recommendations  for  future  research  are  given  in  Section  KIV. 

KII.  PROBLEM  FORMULATION 

Assume  that  the  output  of  the  correlator  y(x)  can  be  expressed  as  in  equation  (5.44), 
where  now  we  allow  the  signal  frequency  fo=coo/27t  to  be  arbitrary  (i.e.,  different  from  the 
center  frequency  fg) 

y(T)  =  C(T)  cos  tOo'f  +  N|^(t)  cos  tOgt  +  Nq(t)  sin  tOgt  (K.  1 ) 
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where 


or 


W(k)  =  S  exp(jQk)  +  n(k) 


(H,) 


n(k)  =  n,(k)+jnQ(k) 
njC*),  nQ(»)  ~  Gaussian  distributed,  U.d 
^{ni(k))‘SS{nQ(k)}  =  0  for  k  =  1,2,...,  G 
var{nj(k)}  =  var{nQ(k)}=  (NoBp  Fi(k)  for  k  =  1,2 . G 

and 

^{n,(0)}=NoB,^{NQ(0)}  =  0 

Note,  that  the  Gaussianess  assumption  for  {n(k)}  weakens  as  k-»  G.  Therefore,  we 
will  concentrate  on  the  observation  interval  k  =  0,1,2,. ..,^G  where  k  =  0.1-0.5.  Under 
these  assumptions,  (G(l  -  X)  »1)  a  good  approximate  formula  for  Fj  (k)  can  be  obtained, 

i.e., 

1 


Problem:  Given  {Wc(k)},  (W,(k)}  k=0,l,..,  XG,  G»l,  decide  between  Hi  (signal 
plus  noise)  and  Hq  (noise). 

A  simple  autocorrelation-domain  (ACD)  algorithm  is  desired  that  can  be  of  low 
complexity  and  very  suitable  for  real-time  implementation. 

Important  Observation:  The  SNR  in  the  ACD,  is  generally  higher  than  that  in  the  time 
domain.  For  example,  take  the  case  where  Cl  =  0.  Given  the  input  SNR  per  band 


Tin  i 


«o“> 


which  typically  is  very  small  (Yjn  =  10"^  lO'^),  we  observe  that  the  local  SNR  in  the  ACD- 
Baseband  Jb  >  is  given  by 


or  equivalently 


y,(k)  .  rl^  (G-w 

k=0,U,MXG. 

if  Yin  =  10-1  (=-10dB).  G  =  103, X  =  0.1.  then 
9.54dB^YBOc)^10dB 

which  indicates  a  significant  increase  in  SNR. 

The  YbOc)  decreases  as  k  increases.  So,  the  problem  here  is  how  to  employ  the  first 

XG  samples  of  {W(k)}  so  that  significant  gains  can  be  obtained.  For  example,  the 
radiometer  uses  the  first  sample  W(0),  and  the  method  of  Kay  [34]  uses  the  first  two 
samples,  W(0),  W(l)  to  form  a  detection  parameter. 

Kill.  DETECTION  METHODS  BASED  ON  AR  MODELING 

Three  different  methods  have  been  investigated  and  their  performance  evaluated  via 
Monte-Carlo  simulations.  The  first  one  is  based  on  complex  data  ({W(k)}),whereas  the 
other  two  are  on  real  data  ({Wp(k)}).  All  methods  employ  an  autoregressive  (AR)  model 
in  which  AR  coefficients  are  used  to  form  a  detection  parameter  in  the  ACD  without 
forming  an  AR  spectrum  estimate.  An  important  result  obtained  during  our  research  is  that 
the  first  method,  which  is  based  on  complex  data,  does  not  exhibit  any  sensitivity  to  angle, 
whereas  the  other  two  methods  do. 

Km  A:  METHOD  I;  COMPLEX  DATA 
(i)  Preliminaries 

Suppose  we  are  given  the  following  set  of  data  samples,  {R(0),  R(1 ),..., 
R(m)}  that  have  the  properties  of  an  autocorrelation  sequence.  Note  that  R(x) »  R*(-t). 


An  equivalent  unique  representation  of  this  set  of  samples  can  be  obtained  in  terms  of 
reflection  coefficient  {kj,  i.e.,  {R(0),  kj,  k2,...,  k^}  where  IkJ  ^  1,  i  =  1,2,..^  [46]. 

This  is  achieved  by  fitting  an  Mth-order  AR  model  to  the  data  {R(0),....  R(M)}  in  the  MSB 
sense,  and  solving  the  resulting  linear  system  of  equations  via  the  Levinson-Durbin 
algorithm.  The  steps  are  as  follows; 


‘  *•*  ►  * 


I  •  •  ™  •  * 

■vaV- 


R(0)  R(-l)  R(-2)  ...  R(.M)1  fl  T  fo 

M 

R(U  R(0)  R("l)  ...  R(l-M)  0 


LR(M)  R(M-1)  R(M.2)...  R(0)  J  La^^iLoJ 

M 

(K.9) 

The  unknowns  are  the  AR  coefficients  tjjg  variance  Pm,  of  the  linear 

prediction  error.  Once  the  above  linear  system  of  equations  is  solved,  then  the  resulting 
AR  spectrum  of  the  data  is 


I  +  £  expj(-Wi)  ^ 

i-1 


.'Vv>:. 

a  •  •  •  »  • 

jkS»._0 


(K.10) 


i.2)  Levinson-Durbin  Algorithm  [50],  [51]. 


There  is  an  elegant,  recursive  and  fast  procedure  for  solving  the  normal  equations  in 
(K.9),  the  so-called  Levinson-Durbin  algorithm.  The  method  requires  ~M2  operations  only 
comparing  to  more  general  methods  that  require  ~\P  operations.  Note  that  the  B  matrix  in 
(K.9)  is  Toeplitz  and  Hermitian.  The  algorithm  proceeds  as  follows: 


t  I 


I _ 


0 


-[R(i)  ♦  I  »B  R(i-n)l 

m-l 

^  ’  Pi-1 

j  -  0,1  .... .1-1 


1—Pi  -  Pi-1  (1-  P) 


(K.11) 


Note  that 


0^'"  '  . " 


one 


and  ajj)  =  0  for  X>i.  So,  in  obtaining  the  solution  for  an  AR  model  of  order  M, 

actually  computes  the  solutions  for  all  AR  models  of  orders  less  than  M,  i.e.,  i-l,2,...,M'^ 
M.  The  quantities  {kj}  1  ^  i  ^  M  are  known  as  reflection  coefficients  or  partial  correlation 

coefficients.  A  necessary  and  sufficient  condition  for  the  stability  of  the  resulting  Mth- 
order  AR  model  is  that  [51],  [52]. 

Pi  >  0, 1  ^  i  ^ 

or  equivalently 

|kil<l,  l^i^M  (K.12) 

It  can  be  shown  that  iff  Kn  >  0  in  (K.9),  then  (K.12)  is  satisfied  [51].  It  can  also  be 
shown  that 


0  S  Pi  <P:  1,  i  =  1,2,...,M 


i.3)  Detection  Parameter 

The  detection  parameter  that  will  be  adopted  corresponds  to  the  geometric  over  the 
arithmetic  mean  of  the  Mth-order  AR  spectrum  given  by  (KIO).  It  is  described  by  [52] 


It  is,  therefore,  clear  the  V^  depends  completely  on  the  shape  of  the  AR  model  spectrum. 
Since  GM  <  AM,  it  follows  that 

O^Vm^I  (K.15) 


It  is  easy  to  see  that  when  Sar^^  (©)  is  flat,  i.e.,  white  noise,  then  vm=  1  (FigureK.l(a)). 
On  the  other  hand,  if  the  spectral  data  spread  is  large,  such  as  the  pattern  shown  below  in 
Figure  K.l(b), 


^AR 


JL 


4- 


(a) 


V-'O 

M 


(b) 


Figure  K.l.  Spectral  spread:  (a)  white  noise  ,  (b)  sinusoid  in  noise 


then  V,^  »  0.  So,  another  way  of  looking  at  V^,  is  in  terms  of  the  flatness  of  the  spectrum 
152]. 

It  can  be  shown  that  V|^  can  be  easily  expressed  in  terms  of  the  reflection  coefficent 
of  an  Mth-order  AR  model  as  follows  [52]. 


Vm 


n  (1-  P) 

1-1 


(K.16) 


However,  a  closer  examination  of  (K.16)  and  (K.11)  reveals  that 


which  is  the  normalized  by  the  signal  energy  variance  of  the  linear  prediction  error 
associated  with  the  Mth-order  AR  model.  Since  Pq  =  R(0),  it  is  apparent  that  from  (K.17) 

and  (K.13),  equation  (K.15)  follows. 

Hie  detection  schemes  described  in  this  section  will  be  based  on  comparison  of  Vj^ 


to  a  specified  threshold  value,  i.e., 


<  Y  (Hi) 


>  Y  (Ho) 


(K.18) 


where  7  is  a  number  between  zero  and  one. 


i.4)  M  =  1  AR  CASE;  Equivalent  Detection  Parameters 


If  we  are  given  the  two  first  samples  only  {R(0),  R(l)}  .  |R(1)  I  <  R(0),  R(0)  >  0, 
then  a  Ist-order  (M  =  1)  AR  model  can  be  obtained  with  the  detection  parameter 


Vi  =  [1  -  Iki  p] 


Vi  - 


R{0) 


(K.19) 


However,  |kjp  can  also  be  used  as  an  equivalent  detection  parameter  as  follows 


>6  (Hi) 

<6  ( hQ  ^ 


(K.20) 


This  is  the  scheme  suggested  by  Kay  [34]  (note  t»  ay'  *  VO  and  shown  to  be  equivalent 
to  a  comparison  of  a  normalized  spectral  peak  with  a  threshold  v,  i.e. 


Sar  (CtT)  > 
oax  _ 

L  Pi  I  < 


V  (Hi  ) 

V  (Hq  ) 


(K.21) 


Another  equivalent  interpretation  of  the  detection  parameter  V|  in  (K.19)  can  be  established 
using  the  determinant  of  the  matrix 


Frco)  iCt)  1 


normalized  by  R^(0).  Specifically,  the  detection  parameter  is 


detCRi  3  r  ^  ^ 

- ^ —  .  1 - 

R(0)2  |_  r2(0)  ^ 


Hi  =  [l-|ki|2]=Vi  (K.22) 

It  is  important  to  note  that  in  the  case  where  {R(0),  R(l)}  are  autocorrelation  samples  of  a 

Gaussian  weakly  stationary  random  process  then  H}  corresponds  to  a  normalized  entropy 

measure.  Under  Hg  (white  noise)  the  entropy  is  higher  than  that  Hj  (signal-t-noise).  This, 

in  the  case  of  a  Gaussian  signal,  in  Gaussian  white  noise,  the  detection  scheme  in  (K.22) 

compares  a  normalized  entropy  measure  of  the  process  to  a  threshold. 

We,  therefore,  conclude  that  Vi(i.e.,  GM/AM  of  the  spectrum  or  normalized 

CO  ^ 

variance  of  the  linear  prediction  error),  la'i  maxISAR(to)/Pi]  and  Hj  (i.e.,  normalized 
determinant  of  the  covariance  matrix)  are  all  equivalent  detection  parameters  for  a  Ist-order 
AR  model  and  complex  data  {R(0),  R(l)}.  Their  generalization  to  higher-order  AR 


i.S)  Mth-order  AR  case:  Generalizations 

If  we  are  given  the  first  M  +  1  samples  {R(0),  R(1),...,R(M)},  then  the  following 
detection  parameters  can  be  obtained  from  the  Mth-order  AR  model: 


(K.23a) 


Pu  H 

-  -  n  n-  |ki  p.  (GM/AM) 

R(0)  i-1 


*M  - 


|a(M)  p.  (Magnitude  of  the  poles) 

M 


Tm  -  I  K"’  p!  138], 
1-1 


Cm  -  I  |ki  p. 

i-1 


(K.23b) 


(K.23C) 


(K.23d) 


^yy'< 

Xv:v:-:-: 

•  •  •' 


(w*) 


PEAKh.  -  aea 

Mirr 


r”(0) 


R(0)  R(0) 


(K.23C) 


(K.23f) 

The  detection  parameters  Rm>  I'm*  And  Cf^  may  be  seen  as  generalizations  of  |kiP  shown  in 
(K.20).  However,  they  are  not  generally  equivalent  with  each  other  nor  with  (K.23e), 
(K.23f)  and  (K.23a).  For  finite  length  data  N  and  order  M  >  1,  each  one  of  the  detection 
parameters  given  by  (K.23  a  -  K.23f)  is  expected  to  have  different  performance.  Some 
asymptotic  results  (N— »<»)  related  with  (K.23a),  (K.23c)  and  (K23e)  are  given  in 
Addendum  I. 


W-  .*  *  »  *  !► 

t-.. 
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(ii)  The  Method 


The  method  employs  as  a  detection  parameter  ((K.16),  (K.17))  and  compares 
that  to  a  specific  threshold  value  as  shown  in  (K.18).  However  the  is  estimated  by 
performing  AR  modeling  in  the  autocorrelation  of  autoconelation  domain  AACD.  This  is 
because  the  the  AR  coefficients  of  an  Mth-order  model  are  computed  using  (K.9)  or  (K.1 1) 
with  "data"  {R(0),...,R(M)}  generated  as  follows; 


1  iG-x 

R(x)  -  -  I  Vl(n)  W*(n*x) 

XG  n-1 


where  W(n)  is  given  by  (K.7).  Since  W(n)  is  the  output  of  the  correlator,  {R(t)} 
represents  an  estimate  of  the  autocorrelation  of  the  autocorrelation  of  the  original  time 
sequence  (i.e.,  the  input  to  the  correlator).  Equation  (K.24)  corresponds  to  the  biased 
autocorrelation  formula,  which  has  been  shown  to  guarantee  a  stable  Mth-order  AR  model 
(i.e.  mi  <  1,  i  ==  1,2,...,M)  [51].  It  is  for  stability  reasons  that  the  biased  estimates  are 
used  here  instead  of  the  unbiased  ones. 

In  summary  the  method  proceeds  as  follows:  a)  Given  the  data  in  (K.7),  compute 
R(t),  X  =  0,1,.. .>1  from  (K.24).  b)  Use  the  Levinson-Durbin  algorithm  in  (K.ll)  to 
compute  Pm  c).  The  detection  parameter  is  =  Pm/R(0).  Note  that  the  first  sample  W(0) 
in  (K.7)  is  not  taken  into  account  when  computing  {R(t)}  in  (K.24). 


(iii)  Performance 

The  detection  performance  of  Method  I  has  been  tested  via  Monte  Carlo  simulations 
for  G  =  1000,  X  =  0. 1-0.5,  and  input  SNR  =  0  dB  down  to  -25  dB.  The  probability  of 
detection  Pq  was  computed  for  Pfa=  assume  statistical  confidence,  a  very  large 

number  of  trials  (1(X)0)  was  performed.  These  results  are  illustrated  in  Figure  K2.  Several 
features  can  be  noted  from  Figure  K2; 


I 

I 


m 


Figure  K.2  Monte  Carlo  Simulation  results  on 
Method  I  with  G=1000 


«w^  -."  <■  »-w? 


a)  For  a  constant  order  M,  the  performance  of  the  AR  detector  I  improves  as  X 
increases. 

b)  For  a  constant  X,  the  performance  also  improves  as  the  model  order  M 
increases. 

c)  For  M=2  and  X.=0.5,  Method  I  outperforms  the  radiometer  detector  ((§))  when 
Yin  ^  dB  and  the  correlator  detector  (0)  when  Yj^  S  -21.4 

d)  For  G  =  1000,  ?[>=  0.95,  Pfa=  M=1,  X  =  0.5,  the  AR  detector  I  will 
require  3.6  dB  and  8.2  dB  less  than  the  correlator  0  and  radiometer  0  detectors, 

respectively. 

e)  There  is  approximately  0.4  dB  improvement  if  M  changes  from  M=1  to  M=2. 
Furthermore,  an  additional  improvement  of  ~0.4  dB  is  achieved  by  changing  the 
order  from  M=2  to  M=7. 

Figure  K3  illustrates  similar  improvements  but  for  the  case  where  G  =.100.  For 
example,  if  M  =  1,  X  =  0.5  then  Method  I  outperforms  the  correlator-detector  (  0  )  for 
7m<  -14dB.  Furthcrmore,Table  K.I  summarizes  the  effect  of  X  (length  of  the  data),  model 
order  M  and  input  SNR  Yin  on  the  distance  d  =  |mo  -  mjl  between  the  two  means  and  on  the 
variances  Oq^  ,  of  the  detection  parameter  Vm-  From  Table  KI,  it  is  apparent  that: 

I  O  9 

a)  ifM=const,  Yin=^o*^st and XT  then  di, 

b)  if  X=const,  Yin=const  and  mT  then  dt  Oq  ♦  *  * 

c)  ifM=const,  X=constand  YuJ^endi,  o^-const, 

The  AR  detector  is  therefore  independent  of  the  noise  variance  (Yin)  for  a  given  order  M 
and  length  of  data  XG.  It  is,  thus,  unnecessary  to  estimate  the  noise  variance  in  order  to  set 
a  threshold  at  the  receiver. 


DETECTOR 


22.807x10^  3.775x10^  31.055x10^  5.801x10  32.532x10  6.394x10 


The  computational  complexity  of  the  AR  cfetector  I  is  compared  with  diat  of  the 
Periodogram  (FFT  based),  using  as  a  figure  of  merit  the  number  of  complex  multiplications 
(which  are  the  most  expensive).  Note  that  a  complex  multiplication  is  equivalent  to  4  real 
multiplications  and  2  real  additions.  The  results  are  summarized  in  Table  KH.  From  this 
table  several  important  features  can  be  noted: 

a)  If  the  length  of  the  data  is  XG  =  2* ,  then  the  Periodogram  detector  requires 
Mu  =  (a  +  1)  •  X.G,  whereas  the  AR  detector  requires  Mu  "  (M  +  1)  XG. 

b)  The  AR  detects  requires  less  multiplications  than  the  Periodogram  for  M  <  a . 

Note  that  a  rectanglar  window  on  the  data  has  been  assumed  for  the  Periodogram  in  Table  KII. 
If  any  other  window  is  used,  which  is  usually  the  case,  then  Mu  (Periodogram)  =* 

(a  +  1)  XG  complex  muiltiplications  plus  2XG  real  multiplications. 

Km.  B  METHOD  H:  REAL  DATA  IN  AR  MODEL 

This  method  employs  real  data  to  form  a  detection  parameter  using  a  second-order 
AR  model.  The  sequence  of  samples  is  given  by  (K.3)  and  (K.5),  i.e., 

Wc(k)  =  S  cos  Ok  +  n,(k)  (H,) 

k  =  0,l,...,G  ;|0|Sn 

W,(k)  =  n,(k)  (Ho)  (K.25) 

with  the  statistics  of  {n](k)}  given  in  (K.7)  and  (K.8).  The  approach  taken  here  models 
{Wc(k)}^  as  a  2nd-order  AR  process  of  which  the  parameters  are  estimated  using  the 
least-squares  (LS)  method  for  spectral  estimation  [46].  The  detection  parameter  is  the 


AR  DETECTOR 


i  ndow  on  the  data  . 


coefficient  of  the  AR  model,  which  corresponds  to  the  product  of  the  poles'  magnitude.  It 


is  given  by 


1  ■ 


1^2 

Rii 


(K.26) 


where 


AG-2  XG“2 

«n  -  I  w2(k^i)  ,  R22  -  I 

k-1  e  c 

AG-2  AG-2 

Ri2  ■  £  WqC  k)Vle(  k  +  1  )  t  ■  E  We(k  +  1)  We(k*2) 

k-1  l(.1 


AG-2 

^2  -  r  We(k)  Wc(k-^2) 
k=l 


(K.27) 


It  is  important  to  note  that  for  the  signal  only  case,  Pn  =  1  exactly.  Furthermore,  for  XG 
large  (XG  2;  100)  the  Burg  technique  [53]  and  the  Yule-Walker  (YW)  method  [54]  for 
spectrum  estimate,  both  give  estimates  very  close,  if  not  identical  to  those  of  LS[46].  This 
is  because  as  XG  — »oo  the  LS,  Burg  and  YW  methods  become  the  same. 

The  detection  scheme  compares  Pu  to  a  threshold  value  y  4  e., 

>Y  (Hi) 


<y  (Ho) 


(K.28) 
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We  have  tested  the  performance  of  this  scheme  for  yi„  =  0  dB  down  to  -20  dB, 

G  =  1000,  k  =  0.1,  0.5,  and  signal  angles  G  «  11.25*  p,  [i  =  1,2,... ,10.  The  sample 
mean  of  the  detection  parameter  Pj]  was  computed  from  1000  sample  trials  under  H 1 .  The 
key  result  obtained  from  this  group  of  simulations  was  that  Pu  is  very  sensitive  to  the 
values  of  the  signal  angle  G(i.e.,  p).  This  high  dependence  of  Pu  (Hi)  on  G  may  cause 
Pii  (Hi)  =  0  not  only  using  the  LS  method  but  also  the  Burg  or  YW  techniques.  The 


•V  ••  •  *.  -N\'. 


results  are  illustrated  in  Figure  K.4(a,b)  for  the  LS  and  Burg  techniques,  respectively.  A 
theoretical  justification  of  this  phenomenon  can  be  established  using  asymptotic  (XG-><») 
results.  As  the  detection  parameter  Py  in  (K.26)  becomes  (LS  or  Burg  or  YW 

method) 

r2(i  )-r(0)  R(2)  „ 

Ptt  "  •  jr 

r2(0)-r2(1  ) 


(K.29) 

with  R(0)  =  1/2  +  Qb,  R(1)  =  1/2  cos  Cl,  R(2)  =  cos  2Cl 

Thus 


N  =  1  -  cos  2Q  -  4  Qb  cos  2Cl,  and 
D  =  2  -  cos  2C1  +  8[Q^b  QbI. 

where  Qb  is  the  "average"  variance  of  the  additive  WGN.  Note  that  if  Qb  =  0,  then  N  =  D 
and  therefore  Pn  =  1.  However,  when  Qb  ^  0  then  N  =  0  if 

1 

cos  2Cl  =  — - -  (K.30) 

1  ♦  »«Qb 

So,  for  those  values  of  Cl  for  which  (K.30)  is  satisfied,  ^n(Hi)  =  0.  Table  Kin  illustrates 
the  cases  where  N  =  0  for  various  Yin  using  (K.30)  and  Figure  K.4.  Apparently  the 
simulation  results  (Figure  K.4)  are  in  good  agreement  with  the  asymptotic  results  (K.30). 
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Table  Kin.  Values  of  £i  *1 1.25*  for  which  Py  »  0  under  various  y„,. 


Yi  n  ( dB  ) 

Q (theory) 

fi  (theory) 

^(simul at1on,Fig.  K.4] 

0 

5x10"^ 

1.81“ 

0.16 

<v/  0 

-5 

5x10-2 

5.68“ 

0.5 

-10 

BSElHi 

16.77“ 

1.49 

^1.5 

-15 

5x10"^ 

35.26“ 

3.134 

rv  3 

-20 

5 

43.63“ 

3.87 

<v4 

Due  to  the  high  dependence  on  the  signal  angle  D,  Method  H  was  not  pursued  any 
further.  It  is  important  to  note,  however,  that  while  Methods  I  and  H  are  equivalent  from 
the  spectrum  estimation  viewpoint,  they  have  diverse  performance  as  far  as  the  detection 
scheme  in  (K.28)  is  concerned. 


III.c  METHOD  m;  REAL  DATA  IN  ARMA  MODEL 


This  method  is  a  modification  of  the  LS  approach  and  employs  the  same  real  data 
record  {W<.(k)}‘»®.  It  models  the  data  as  an  ARMA  (2,2)  process  from  which  the 

magnitude  of  the  second  AR  coefficient  serves  as  the  detection  parameter,  i.e.. 


1*2  1  -  h 


*21  -Zj  t2 


*1 1  *22"*12*21 


(K.31) 


where 


zii  - 


M-2  M-2 

I  Wc(k-M)  Wc(k*3)  ,  *22-1  Me^X)  Wc^^***) 

k-l  k-1 


M-2  M-2 

*12  -  I  WpCk)  We(k>3)  ,  221  "  I  We^*****) 

k-1 


k-1 


M-2  M-2 

I  We(k>l)  We(k*3)  ,  tg  -  I  «c(k^2)  Mo(k*a) 

k-1  k-1 
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In  the  "signal  only"  case,  |a2l  =  1  exactly.  Furthermore  under  (signal  plus  noise),  as  the 
length  of  the  data  record  M^«o  then 

Zji  =  1/2  cos  2(2,  Z22  -  ^  ^12  = 

Z2i  =  1/2  cos  3(2,  tj  =  1/2  cos  3Q  and  t2  =  1/2  cos  4£2. 

Thus, 

l/8Ci-eoa2Q]  I 

- f  t. 

t /8[l-cos2fl  I 
i.e , 

asymptotically  (M->o«>),  the  parameter  |a2l  is  free  of  the  AWGN.  However,  the  price  paid 
is  that  under  Hq  equation  (K.31)  becomes  a  very  bad  estimator  of  |a2l  which  as  M->oo 
becomes  undefined  (|a2l  =  — ^  ).  To  remedy  this  problem,  we  sacrifice  a  part  of  the 
method's  good  performance  under  by  using  as  detection  parameter 


^  *21  *1  "  *11  t2 

Pjjj  .  - 

*11  *22  "  *12  *21  *  Q 


where  Q  = 


Wc^CO) 


The  detection  scheme  is  the  same  as  that  shown  in  (K.28). 


(K.32) 


We  have  tested  the  performance  of  this  approach  for  Yu,  =  -10  dB  down  to  -20  dB, 
G  =  1000,  X=0.1,  0.5  and  signal  angles  SI  =  11.25*  p,  p  =  1,2,...,10.  The  was 
computed  for  Ppy^ »  lO-^  and  the  results  are  illustrated  in  Figure  K5.  From  this  figure,  it  is 
apparent  that  the  performance  of  Method  in  is  determined  by  the  angle  signal  £2,  having  as 
a  lower  bound  curve  the  AR  (min)  and  an  upper  bound  curve  the  AR  (max).  Thus,  while 
the  angle  dependence  of  Method  in  is  much  less  than  that  of  Method  n,  it  is  still  large 
enough  to  degrade  the  performance  of  the  detector.  We,  therefore,  did  not  pursue  any 
further  Method  m. 
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Figure  K.5  versus  in 

method  III,  with  varying  U, 

G  =  1000,  A=0.1,  0.5  and 
Pp^=10”l.  AR(max)  corresponds 
to  the  "best"  Q  and  AR(min) 
corresponds  to  the  "worst"  Q  . 
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IV.  CONNECTION  WITH  OTHER  DETECTION  SCHEMES 


Given  the  data  r(k),  k  »  0,1,2,..  in  the  time-domain  where 


r(k) 

=  exp  +  e(k) 

(Hi) 

r(k) 

=  e(k) 

(Ho) 

(K.33) 

various  detection  schemes  have  been  developed  to  decide  between  Hq  and  Hj. 

1)  Radiometer  The  detection  statistic  is 

1  L 

W(0)  -  -  X  |r(k)  |2  (K.34) 

L  k-0 


which  is  the  first  lag  of  the  autocorrelation  of  {(k)). 


2)  AR(l)  Method  of  Kav  [341:  The  detection  statistic  here  is  the  magnitude  of  the 
parameter  of  a  Ist-order  AR  model,  as  explained  in  section  III.A  (i.4).  The  method 
estimates  the  first  two  autoconelation  lags 


and  fc»ms  the  parameter 


1(1  >  p 
r(o)  I 


(K.36) 


which  is  compared  to  a  threshold. 


3)  Adaptive  AR  Detection  Scheme  (551:  Considering  the  binary  detection  problem 
r(k)  =  s(k)  +  e(k)  (Hj) 

r(k)=  e(k)  (Hq)  (K.37) 

the  method  models  the  unknown  signal  as  an  AR  process 


s(lc)  +  I  ai  s(k  -  i)  =  u(k) 
i=l 


and  uses  as  detection  parameter  the 


SNR=  JL 


(K.38) 


The  method  estimates  from  the  observations  {r(k)}  the  {aj}**,  &2  using  a  nonlinMr 


m  ot  Signs 


Considering  the  binary  detection  problem 
r(k)  =  s(k)  +  n(k)  (Hj) 
r(k)  =  n(k)  (Ho) 

where  {s(k)}  is  assumed  to  be  stationary  Gaussian  random  process  with  a  broadband 
unknown  PSD,  the  method  models  {r(k)}  as  an  AR  process  of  order  p,  i.c., 


r(k)+  I  ai  r(k  -  i)  *  u(k) 
1-1 


and  uses  as  detection  parameter 


t 
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which  actually  is  the  GM/AM  of  the  signal  spectrum. 


S)  RF  Correlator-Detector  fSection  5.31 

This  is  a  detection  scheme  that  employs  all  autocorrelation  samples  of  {r(k)}i.e., 

W(x)  =  -— r(k)r*(k  +  2) 

L  K“0 

X  =  0,1.— 

by  forming  the  following  detection  parameter. 

,  XG 

R(0)  =  _^  I  W(t)W*(t) 
x  =  1 

which  is  actually  the  first  lag  of  the  autocorrelation  of  the  autocorrelation  sequence.  Thus, 
the  Correlator  Detector  may  be  seen  as  the  radiometer  operating  in  a  different  domain 


ADDENDUM  I 


V. 

In  the  asymptotic  case  (XG  ->«>),  the  "data"  generated  by  (K.24)  becomes. 

Hq:  {R(0)=Q,  R(1)=...=R(M)=  0} 

Hi:  {R(0)=1+Q.  R(l)=e-j«  , ....  R(M)=e-iMn  }.  (1.1) 

The  AR  parameters  generated  via  the  Levinson  recursion  take  the  values:  (Mth-order 
Model) 


Ho:  a(M)  =  0,i=l,2,...J^ 

kj  =  0,  i=l,2,...,M 

Pi  =  Q,  i=0,l,...,M  a.2) 


1 


a(M)-  - 

- ^ -  e'Jox,  A-1,2 

X 

(Q  +  M) 

1 

kj- - 

e  ,2,  . 

(Q 

+  i  ) 

(Q  +  i  +  1  ) 

Pi-  ♦  Q 

(Q^l) 

a.3) 

Therefore,  the  detection  parameter  Fj^  (Ketchum,  [38])  (K.23c)  becomes 


Ho:  Fm^O 


Hi:  Fm  = 


M 


M' 


(Q+M)2 


and  thus  the  distance  ‘>'m=  irM(Hi)-FM(Ho)|  is  given  by 


a.4) 


M 

(Q*M)2 


On  the  other  hand,  Vj^  becomes 
Ho:  Vm  =  1 


Q(Q  +  M  +  1 ) 

Hi:  Vm=  - 

(Q  +  1)  (Q-^M) 


and  therefore,  the  distance  vm=|V,^(Hi)-V^i(Ho)|  takes  the  value 


M 

(Q+1)  (Q+M)  (1.5) 

Comparing  with  it  is  apparent  that  Vi4>  ^  for  M>  1  and  Vj^  =  for  M=L. 

For  the  asymptotic  case  (XG-*«>),  we  can  also  find  equivalent  detection  parameters 

to  a  normalized  AR  spectrum  peak.  If  we  define  as 

(M)  ^ 

Sar(w)  a  1 


Pm  ^(w)  P 


(1.6) 


where  D(co)  is  obtained  from  (1.3)  and  (K.IO)  as 


M  1 

1  ^ —  e-J[ti>*n]x 

X-1  (Q.^M) 

then 

1  M  ’  M 

J(W)  -  )D(6#)  p  .  1 - ^  eJ(W‘'QA)- - - -  y  eJ(W*a)4i 

(Q  ♦  M)  X-1  (Q  ♦  M)  jiti 


’  M  M 

^  -  1  1 

(Q  ♦  M)2  xty  i.t'i 

By  differentiating  and  forcing  ■■■  —  =0,  we  obtain 

9U) 


eJ«*>*a)(u-x) 


1 


H  M 


♦  M  X*1  W»1  li»l 


After  some  algebra,  we  find  that  the  solution  of  (1.7)  which  corresponds  to 
minimum  J(ca)  is 

W*  ♦  Q  -0 

and,  thus,  the  spectrum  peak  is  given  by 


I  I 

•  .*«  V- 


• 


ipicj*)  (2 


(,  -  -JLy 

\  Q  ♦  M/ 


We  conclude  that  equivalent  detection  parameters  to  (1.8)  can  be 


n  M  M  1  ,  , 

-  •  2  I  or  - -  -  rs,  or  —  -  la^«)  |. 

+  M  i-1  ^  (Q  ♦  M)2  Q  ♦  M 


«  1  i 


For  example,  when  M=2  equation  (1.7)  becomes 


♦  Q)  -  0 


r  ’  1 

2  3in2(OI  +  Q)  +1  -  - -  sin( W  ♦  Q) 

L  Q  "  2j 

.n(o>  ♦  n)ri4  cos(GI  •*•  Q)  ♦  1  -  - = -  j-  0 

L  •  Q  -  2j 


Therefore  ^ - 

|d(W*)  P 

(Q  ♦  2)2 


1  yQ  ♦  1v 

cos((i)*  ♦  Q)  - - I  I 

\Q  ♦  2/ 

where  r>2. 


('-T^y 


(Q  ♦  2)2 


jD(W*)  p  (Q  ♦  r)2 
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